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Chapter  I 
INTRODUCTION 


!•  Background  of  the  Problem. 

From  the  moment  a  projectile  leaves  th-  muzzle  of 
the  gun  which  launches  it  on  its  flight  "down  range*,  it  is 
subjected  o  a  complicated  variety  of  forces.   The  st  dy 
o*  the  effects  of  these  forces  upon  wh  re  the  projectile 
goes,  and  upon  how  long  it  takes  to  get  there,  has  for  cen- 
turies fascinated  mathematicians  and  physicists  ali;  e. 

Modern  guns  are-  rifled,  aad  hence  impart  to  the 
projectile  a  spinning  motion  which  affects  its  "stability" 
and  thus  influences  its  flight.   But,  even  if  one  is 
wiling  to  defer  the  tempting  study  of  the  gyroscopic 
properties  of  a  spinning  shell,  on  the  ground  oha*  th  y 
are  jecondary  in  importance  to  he  mo  ion  o  its  center 
of  mass,  one  finds  t^at,  to  compete  a  trajectory  o  t  e 
ma  -  center,  one  must  solve  a  di  ferential  system  involving 
a  number  of  functions  which  depe  d  upon  the  place  at  which 
the  projcti  e  happer.  to  be  as  we  1  as  upon  th*.  time  when 
it  happens  to  be  the^e,  not  to  menti  n  the  spe-    t  wh.ch 
it  happens  to  be  mo  ng. 

For  example,  the  air  resistance  whic-  a  buile 
nccun  ers  is  dependent  not  only  upo.  the  proj  ctil  's 


3?ee»  relative  to  the  air  bu  also  upon  ti.e  ntiOspheric 
d~*isit„  and  tempera1  ire,  and  th-.«-e  lsttei  quentl' ies  vary 
from  plaee  to  pie  t   and  from  ins  ant  to  inetont  aa  the 
shell  traverses  te  path.   The  alr-soeed  of  the  projectile 
la  itself  dependent  upon  the  speed  and  direction  of  the 
wind,  !>nd  v'iuda  ere  notoriously  "contrary". 

At  th  s  writing,  not  even  with  the  aid  of  high- 
speed omputing  machinery  haa  it  bee.  practicable  to  obtain 
a  complete  tabulation  of  aolutiona  of  the  differential  equa- 
tions of  motion,  eovoring  all  possible  sets  of  initial  con- 
ditions, wind-structures,  density  and  temperature  distribu* 
tions  eto    The  evice  which  balllstl  ien-  have  employed 
to  overcome  thia  difficulty  ia  that  of  integrating  the  equa- 
tions under  oertaln  relatively  simple  "standard*,  or  •nor- 
mal" conditions  (such  as  zero  wind,  a  n  n-rot*ting  earth, 
ate),  and  then  computing  a  1  lnear  approximation  to  the 
actual  effect  (on  range,  on  tlme-of-  light,  or  on  some 
Other  trajectory  "element")  due  to  the  "disturbances"  (de- 
partures from  "normal"  conditions'. 

*.  First  Differential  Effects. 

The  mathematical  theory  underlying  the  computa- 
tion of  these  linear  approximations,  or  "differential  ef- 
fects" as  they  are  usually  oalled,  was  developed  by  G.  A. 
BIJaa  while  he  waa  at  the  Aberdeen  Proving  Grcund  during 
World  War  I.   some  of  his  results  were  published  in  refer- 


ences  [l]  nd  [2]  of  the  bibliography  at  the  and  of  this 
paper;  others  appear  in  referenoes  [3],  [4],  and  [5]. 
Other  investigators  have  provided  modifications 
and  improvements,  both  in  the  theory  and  in 
its  applications.   The  most  recent  such  contribution  is 
a  paper  by  B.  J.  McShane,  entitled  *The  Differentials  of 
Certain  Funotionals  in  Exterior  Balllstlos*,  and  soon  to 
appear  in  the  Duke  Mathematical  Journal.   In  this  paper, 
which  will  hereafter  be  designated  as  [MoS]»  Mcshane  shows 
that,  for  moat  trajectories,  the  "norm"  (or  measure)  of 
the  disturbance  as  used  by  Bliss  in  his  proofs  can  be  re- 
placed by  a  simpler,  more  natural  "norm". 

In  order  to  make  clear  the  distinction  between 
the  norm  used  by  Bliss  and  that  used  by  McShane,  we  shall 
quote  below  an  excerpt  from  the  introduction  to  the  paper 
[Mcs] -  the  manuscript  of  which  Professor  McShane  very 
kindly  loaned  the  writer  in  advanoe  of  its  publication. 

"To  be  specific,  let  us  restrict  our  attention  to 
the  effect  of  wind  on  range.   The  'normal*  condition,  un- 
der which  trajectories  are  computed,  is  zero  wind.   The 
range  is  a  functional  of  the  wind.   If  the  wind  is  sup- 
posed to  be  a  continuous  function  of  altitude  alone,  we 
are  considering  a  real-valued  functional  on  a  space  of  con- 
tinuous functions .   In  this  space  we  can  introduce  a 
familiar  norm,  the  norm  of  a  continuous  function  w(.)  being 
defined  to  be  the  least  upper  bound  of  its  absolute  value. 


<: 


Then  the  classical  definition  of  a  (first)  differential  can 
be  applied;  it  is  a  linear  functional  of  w(.j  which  differs 

>m  the  change  of  range  /\I  due  tow(»!  by  a  quantity 
whose  ratio  to  norm  w(.)  approaches  zero  with  norm  w(,}. 
But  the  existence  of  a  differential  in  this  sense  has  noi 
heretofore  been,  established*   lo  the  best  of  my  knowledge.. 
lie  nly  adequate  matheraati  :.l  is  ussion  of ■  differential 
affects  s  due  to  G.  A.  Bliss.   Bliss,  hcwev3rs  does  n 
•  'he  norm  mentioned  above.   Ins-ead,  ha  defines  the 
norm  oj  w  .)  to  be  the  greater  of  the  least  upper  bound  of 
|w|  and  the  lea^t  upper  bound  of  the  absolute  value  of  its 
derivative  with  respect  to  altitude  y:  that  is,  he  uses 
norm  of  the  Banach  space  C'  instead,  of  the  norm  of  the 
vace  C.   Now  with  certain  disturbances  this  would  be 
quite  unobjectionable.   For  instance,  the  Cor lolls  force 
a  projecti  e  is  a  disturbance  whose  alur  and  rates  01 
change  with  respect  to  position  and  velocity  are  comput- 
a  le,  and  the  inclusion  of  the  rates  of  ciange  in  the  defi- 
nition of  the  norm  produces  no  annoyance,   But  he  'ind  is 
an  experimentally  determined  function  and  the  experiment 

-nishes  only  the  mean  values  between  a  finite  numb?r  of 
successive  altitudes.   The  experimental  result  may  be  ex- 
pected to  differ  from  the  exact  wind  by  an  error  whose 
norm  in  the  space  C  (i.e.,  maximum  absolute  valued  is 
small;  but  we  have  no  assurance  that  the  rate  of  crange  of 
wind  wih  altitude  does  not  at  some  places  differ  widely 
from  the  experimental  estimates    In  order  that  the 


mathematical  treatment  of  differential  effects  shall  be 
applicable  to  the  experimental  methods  involved,  it  is  of 
great  importance  that  the  norm  be  determined  by  the  magni- 
tude of  |w|  alone,  and  not  depend  on  the  derivative  w1 . 
......  We  have  seen  that  disturbances  can  be  classi- 
fied into  two  aggregates;  for  those  of  one  class  it  is 
harmless  to  admit  the  derivatives  in  the  definition  of  the 
norm;  the  others,  specifically  wind,  temperature,  and  den- 
sity, are  to  be  regarded  as  functions  of  altitude  alone, 
and  their  derivatives  are  to  be  excluded  froci  the  defini- 
tion of  the  norm.   It  is  possible  to  avoid  the  annoy  roe 
of  considering  two  separate  classes  of  disturb  nces,  and  at 
the  same  time  to  gain  slightly  in  generality,  hy   defining 
the  norm  of  a  disturbing  |f unction  to  be  the  gr  atest  of  the 
maxima  of  the  absolute  values  of  the  function  itself  an5  of 

r 

its  partial  derivatives  with  respect  %o  all! variables  ex- 
cept  the  altvt-de  y.    •    .    .     .        •        *he  norm  being   defined, 
....    ,,  with  given   fn  rmal*   conditions   and  given  dis- 
turbances,   ......  do  departures  from  'normal'    condi 

tions  produce    Mfferentiabl^  effects?       We  find  that  they 
do      ,     Is 'the  difference  between  differential  ef- 
fect and  actual   change  an  infinitesimal  of   second  ordar  in 
the  norm?       We  find  that  it   is  it   the  trajectory  is  every- 
where ascending  or  everywhere  descending;  otherwise' it  may 
not  be  as  small  as  second  order," 


3.  Need  For  A  New  Theory.  Summary  of  Results. 

For  most  purposes,  the  use  of  differential  effects 
in  place  of  actual  effects  of  disturbances  has  produced 
quite  satisf acto  ^  results.   Recently,  however,  evidence 
has  accumulated  to  support  the  belief  that  linear  approxima- 
tions to  effects  of  cer  ain  disturbances  (e.g.,  a  wind;  or 
the  combined  presence  of  wind  and  abnormal  temperature)  do 
not  give  sufficiently  acourate  estimates  of  the  actual  ef- 
fects.  Because  ij  the  need  of  improved  estimates,  and  be 
cause  the  literature  contains  no  mathematical  theory  o 
second-or  er  effe  s  if  disturbances,  Professor  McShan 
suggested  to  the  writer  that  be  undertake  the  developmen 
of  such  a  theory. 

~The  present  pap  r  consists  of  a  portion  of  the  re- 
sults of  that  nvestigation.   It  includes:  (1)  a  proof  of 
the  existence  of  a  second  differential  of  the  mapping  de- 
fined by  uhe  dif  .eren^ •  al  equations  of  mofon  of  the  mass- 
center  of  a  proj  ct'l  i  (2)  a  demonstration  hat  this  se  ond 
differential  is  expressiblo  in  terms  of  the  second  variations 
of  the  coordin^t;-   f  the  center  of  mass,  (3)  a  derivitlon 
pf  the  system  of  differential  equations  satisfied  by  the 
second  variat  on-  ,  a  system  which  ;as  the  same  hompgei  eous 
terms  as  d  es  the  classical  system  of  first  variat:. on  equa- 
tions; A)   a  formula  for  expressing,  a  solution  of  toe  dif- 
ferential system  or  second  ariat.ons  in  , erms  of  double 
Stieltjes  i,  terra's,  (5)  a  simple  formula  for  oa..oula,ing 


the  second  variations  of  trajectory  elements  from  the  second 
variations  of  t;ie  coordinates  of  tne  >a88a-center ;  and  (6)  a 
theorem  on  second  differentials  enalowus  to  ;..c;hene's  main 
theorem,  (11.7 )   of  [l.cj]  ,  on  first  differentials. 

The  nor-  of  disturbance  employed  In  obtaining  the 
first  :ive  resultr  above  is  defined  in  terms  of  the  absolute 
values  of  the  disturbance  functions  and  their  partial  de- 
rivatives of  the  first  and  second  orders,  with  respect  to 
all  the  variables  except  the  independent  one,   A  somewhat 
simpler  norm  ir  used  in  establishing  (o). 

4.   Fields  for  further  Investigation. 

It  har.  Ion  been  known  that  (first)  differential 
effects  are  expressible  in  terms  of  (single)  Stieltjes  in- 
tegrals.   (By  "otieltjes  integral"— both  the  single  ^nd  the 
double  varietie.  — v.e  mean  the  one  oft.^n  referred  to  as  the 
"Aiemann-Stieltjes"  integral,  not  the  "Lebesgue-Stielt jes" 
integral.)   This  fact  has  made  possible  the  use  of  a  "bal- 
listic wind" — e.   fictitious  constant  wind  which  has  the  same 
differential  effect  on  range  as  does  the  actual  renge-wind 
wx(y)«   (A  "ballistic  density"  and  a  "ballistic  tempera- 
ture" are  similarly  defined.)   The  ballistic  wind  is  com- 
puted in_  the  field  from  observed  wind-data  by  the  use  of 
so-called  "wind  weighting- factor  functions"  which  have  been 
pre-computed  at  the  ballistics  laboratory.   iach  wind 
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wei ght ing-ftfc tor  function  is,  by  definition,  the  rtitio  of 
the  differential  effect,  g(k)„  0  «  k  =  1,  of  a  range-wind 
waich  is  zero  at  all  levels  above  kY  (Y  being  the  altitude 
of  the  summit  of  the  trajectory)  and  +1  at  all  level*  be- 
low :-Y,  to  the  differential  effect  g(l)  of  a  range-wind 
whi'.h  is  one  length-unit  per  t'nas-unit  at  all  levels.   The 
graph  o;.  a  weighting-factor  function  is  known  as  a 
"weight 'ng-f act or  curve",  and  i>  plotted  at  tLe  laboratory. 
Thus  the  fact  that  (first J  diff ■  ential  effects  a>.e  ex- 
pressible as  (single)  Stieltjes  integrals  has  rendered 
possible  a  method  whereby  practically  all  the  computations 
for  differential  effects  can  be  performed  by  a  group  of 
people  who  can  be  situated  far  from  the  combat  theater,  and 
the  results  of  their  calculations  can  be  recorded  in  such  a 
manner  that  applications  may  be  made  in  the  field  quickly 
and  by  means  of  very  simple  computations.   (For  a  more  de- 
tailed account  of  weighting  factors  and  their  applications, 
see  References  [5]  ,  [6]  ,  and  [7]  . ) 

Now  the  results  of  the  present  investigation, 
which  enable  us  to  express  second  differential  effects  in 
terms  of  doubl e  stieltjes  integral*,  lead  us  naturally  to 
conjecture  that  it  is  possible  to  devise  a  method  of 
weighting-faotor  functions  for  second  differential  effects 
analogous  to  those  for  first  differential  effects  described 
above.   But  no  attempr  has  been  m*de  to  do  this  in  the 
present  paper.   In  fact,  such  an  effort  does  not  appear 
feasible  until  the  formulas  developed  herein  shall  first 


have  besn  applied  to  the  ballistic  equ?.:ions    (aifferential 
equations  of  motion  of  the  mass-oenter)  themselves.       This 
latter  task  is  Itself  one  of  some  magnitude,   as  will  become 
increasingly  evident  in  later  chapters         It  has  been  post- 
poned as   constituting  an  appropriate  objective  for  a  '^.ter 
investigati on, 

5o  Remarks  on  Notation. 

After  experimenting  with  several  notations*  the 
writer  concluded  that,  for  the  purpose  of  tnis  paper,  the 
most  satisfactory  system  to  use  would  be  a  slight  modifies- 
tion  of  the  notation  employed  in  [mcs] •   Since  the  latter 
paper  is  not  yet  available  in  published  form,  it  has  seemed 
advisable  (for  the  convenience  of  the  reader)  to  quote  sev- 
eral passages  from  the  manuscript  form  of  fucs] .   These 
will  be  clearly  indicated  by  quotation  marks,  and  grateful 
acknowledgment  is  hereby  made  for  Professor  McShane's  per 
mission  to  make  these  quotations* 

We  shall  use  the  summation  convention  with  regard 
to  repeated  indices  throughout  the  paper,  except  in  certain 
parts  of  Chapter  VI  where  the  departure  from  this  notation 
is  clearly  indicated. 
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Chapter  II 

THE  EQUATIONS  OF  MOTION.   A  MORS  GENERAL 
DIFFERENTIAL  SYSTEM. 


60  The  Equations  of  Motion. 

In  this  section  we  shall  develop  the  differential 
equations  of  motion  of  the  mass«eenter  of  a  projectile,  sub- 
ject to  certain  simplifying  assumptions,  in  the  form  in 
which  they  are  usually  encountered  in  ballistic  literature. 
For  the  notational  reasons  indicated  in  $  5,  we  quote  the 
following  passage  from  [Mcs] • 

"Let  us  ohoose  a  rectangular  coordinate  system  in 
which  the  direotion  of  the  earth»s  gravity  field  at  (0,0,0) 
is  along  the  negative  y-axis.   The  system  is  supposed  fixed 
with  respect  to  the  earth.   In  order  to  avoid  inessential 
annoyances,  we  shall  select  some  number  B  and  deoide  to  con- 
fine our  attention  to  trajectories  along  which  the  in- 
equalities 

(6*1)   |ac|<B,  |y|<B,  |z|<B,  |vx|<B,  |vy|<  B,  |vz|<B, 

are  satisfied  and  whose  time  of  transversal  is  less  than  B» 
The  set  of  all  points  satisfying  {6.1)  will  be  named  R. 
This  restriction  is  not  unrealistic;  at  any  particular 
epoch  there  is  an  upper  bound  to  attainable  velocities,  and 
if  trajectories  depart  more  than  a  few  hundred  miles  from 
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the  earth's  surface  problems  of  ballistics  blend  into  those 
of  orbit  computation.   The  components  of  -velocity  of  the 
center  of  mass  of  the  projectile  (which  for  brevity  we 
shall  call  the  components  of  velocity  of  the  projectile) 
will  be  denoted  by  (vx,  Vy,  vz)  ,  and  the  length  of  the  ve- 
locity vector  by  v.   The  components  of  wind-velocity  will 
be  denoted  by  (wx,  wy>  wz) «   It  is  convenient,  though  not 
absolutely  essential,  to  assume  that  Wy  is  identically  zero. 
The  velocity  of  the  projectile  with  respect  to  the  air  then 
has  components 

(6.2)    UX  =  VX-WX,     Uy  »  Vy-Wy  =  Vy,     Uz  -  VZ-W2  . 

The  length  of  the  vector  (ux,uy,uz)  we  denote  by  u.   We 
assume  that  the  components  of  wlnd-veloc ity,  the  air- 
density  P  ,  and  the  absolute  temperature  ©  ere  func  .:. ens 
of  position  only,  independent  of  time.   We  shall  assume 
thet  for  a  projectile  of  given  shape,  size,  and  mess,  mov- 
ing through  air  of  given  chemical  composition,  the  princi- 
pal aerodynamic  force  has  direction  opposite  to  the  vector 
(ux,Uy,uz)  and  magnitude  muE,  where  m  is  the  mass  of  the 
projectile  and  E  is  a  function  of  u„  P  ,  and  *),  and  is 
continuous  together  with  its  partial  derivatives  of  first 
and  second  order.*   If  the  gravitational  acce       a   were 
of  constant  magnitude  and  direction  throughout  .;,  it  would 


For  a  more  detailed  discussion  of  the  nature  of  the  func- 
tion E,  see  references  [5],  [o]  ,  and  [7]. 
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produce  a  force  (0,  -iug,  0)  on  the  projectile.  Thus  the 
force  on  the  projectile  had  components  approximately  equal 
-mEujc,  -nifiuy-mg,  -mSuz,   All  the  other  forces,  whether 
rodynamic  or  due  to  the  .*    e  i  i   ravity-field  or 
oriolis  forces,  we  lump  lct     i:    r    and  designate  by  the  sym- 
ol  (max, may,  az)  •        FO]  pn  purposes  it  is  immaterial 

how  these  forces  can  be  calculi  ec  ;  .  ,  is  enough  that  they 
can  be  found,  and  that  they     ordii-jrily  considerably 
smaller  than  'drag'  mEu  aid   •       .-.-.     We  can  now  write 
the  equations  of  motion  in  the  form 


(6.3) 


dt 

-  Euy   -   g 

dt 

-    iSux    +     &y 

..      -Tn         J-     1 

-  ■  vy. 


. 


dt        z    Zl  * 


"Henceforth  we   shall  assume   that    there  is   a  posi 
vc  number    *l~   such  that   everywhere  in   R  the   inequality 

(6.4)  ay   -      g  -    ^ 

Then  at  every  value  of  t  such  that    uy  ■   vy  =   0, 
we   have 


dtg 
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This   implies   that  whenever     Vy     vanishes     y     has   a  proper 
maximum.        Then     y      can   have  no   relative  minimum.      Between 
two  proper  maxima  a  minimum  must   occur;   therefore  there   can- 
not  be  two  times  at  which     Vy     vanishes.        If  on  a  trajec- 
tory there  is  a  time      ts»    at  which      vy  »  0,    the  point 
(x(ts),    y(ts),    zUs))    is  called  the    summit.        For  t   <    ts 
we  have     vy   >    0;   the   corresponding  arc  of   trajectory  is  the 
ascending  branch.       For  t    >  ts  we  have  vy  <;   0;  the    corres- 
ponding arc  is   the  descending  branch,.       Artillery  trajecto- 
ries,   for  fire   at   ground  targets,    ordinarily  have   summits. 
Bomb  trajectories  in  level  or  dive  bombing  have  the  descend- 
ing branch  only;    anti-aircraft  trajectories  usually  consist 
of  an  ascending  branch  only. 

"The  assumption    (6.4)    is  a  very  mild  one  for  con- 
ventional trajectories.       However,    by  adopting  it  we  exclude 
projectiles   •   .    •  which  are   capable  of  sustained  level 
flight." 

7 .     The  Equations  of  Motion  as  a  Special   Case  of  a  More 
General  Differential  System. 

The   system   (6.3)    is  only  one   of  several  ways  in 
which  the  equations  of  motion  of    the  mass-center   can  be 
written.       Other  useful   forms  are  obtained  by  using  alti- 
tude or  slope  as  independent  variable,    instead   of  time. 
In  each  of  these  forms,   however,    the   right  members  of  the 
equations   are   rather   complicated   expressions   involving  the 


compound  function  E.        But  all  these  forms  of  the  equations 
'  motion  have  in  common   the  property  of  being  a    special 
i  of  a  differential  system  which  can  be  written  quite 
simply  in  the  notation  introduced  in    [Mcs] .        Our  investi- 
gation,   t     refoie,   will   take   the  form  of  a   study  of  this 
:   general  differential   system,   not   only  because   of   the 
greater  ecopc  r    of   notation  but  bIso  because  the  theory 
which  we   shall  c.evelop  will   then  be  applicable  to  several 
forms  of  tbe    equations   of  motion   instead  of  being  confined 

aly  to  the  system   (6.3) .        For  this  reason,  we  again   quote 
from   [KcSJ ,    this    time   from      f  4  of  that   paper. 

Tre   variable  v/hi  ch  is   chosen  as  independent  var- 
)le  in  the  equations   of  motion  shall   be    called     s    ,    and 
11  be  assu     d   to   vary  over  some  interval    (si,sg)» 

e,    this   corresponds   to  assuming  a  bound  on  the   times 
flight  of  all   projectile      considered;    if,  say  y  or  x  is 
?n  as   independent      a     a  )ie ,  it  follows  from  the    assump- 

lat  the   trajectori    o  lie   in  a  bouided  part  of  spa 
depende.i .     rar      bles  we    shall        q1,.    .    .  ,qn   ,      and  the 
(q   »•        .    q11)   we   shall  abbreviate   to      q     when  we 
Sup    r.c      pts   aff ■ x   d   directly  to  letters  are 

ve ;  expone   ts,  when  we  have   occasion  to   i 
be  af   lxed   to  parentheses.        Thus    (q3)2     wou! 
square    of   the   third  component   of  the  n-uple      q 
dent  v   ria   les   ra    ge    over  some  open  set   Ro   in  n- 
3n   the     ?et      (si,      2)  )<  R0   ,    consisting  of  all 
^les    (s,q)   with  si  *    s    «  sg  and  q  in  R0   ,    certain 


3 


0 
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functions     y>  ■*■    (s,qj,    .    .    .    ,     lp  m    (s,q)   are  defined;   these 
will   always  be  assumed   continuous  except  perhaps  at  a  finite 
number  of  values  of      s,    at  which  they  may  have   simple   Jump 
discontinuities.        We    shall  also  assume   that   all  the      ^ 
considered  are   less  in   absolute  value   than  some  uniform 
bound  B;   this   corresponds   to  deciding  in  advance  that  we 
shall   only  consider  winds,    velocities,    etc.,  which  do   not 
exceed  a   chosen  bound.        On  the  set   in  m-space  defined  by 
Ir^J    ■  B    (  **  ■  l,...,m)    there  are  n  continuous  functions 
pl(r ),...,   Fn(r).       We  seek  solutions  of  the  equations 

(7.1)  dqVds   -  F1!  <f>  (s,q(s))) 

-  F*(  if1(s,qV...,qn),...,ipni(gtql q11)), 

(i»l,...,n) 

with  initial  values 

(7.2)  q1(si)   -  qo  (i«l,...,n) 

lying  in  Ro  and  with  q(s)   [i.e.,   (q1(s) , . . . ,qn(s) )  ] 
lying  in  R0  for  si  -  s  ■  sg. 

"The  equations  (7,1)  and  the  initial  conditions 

(7.2)  are  determined  when  we  specify  the  system  of  numbers 
and  functions 

(7.3)  (q£.. ..,<£,   f1 f  m,F1,...,Fn). 

"By  well-known  theorems      [see,    e.g.,   E.   J. 
McShane :      Integration,    (Princeton  Univ.   Press),  p.   344.]    on 
differential  equations,  for  each  system   (7.3)   with  the 
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properties  previous    specified,    equations    (7.1),    (7.2)    have 
a    solution  defined   on  some   sub-interval    (si,    S3)      of 
(sl>s2)      an^   lyin8   i-n  &o»        But  we    cannot   guarantee   that 
the   solution    is  unique.        Systems  lacking   unique  solutions 
are   of  no    interest  to   us,    so  frora  the   aggregate   of  all   sets 

(7.3)  we  select  the  subset  (J)   consisting  of  all 
(q0,  w   ,F)  which  furnish  unique  solutions. * 

Since  on  the  set  (J)  the  solutions  may  not  be 
continuable  to  s  =  S2»  we  select  from  (j)  a  subset  (Jo  de- 
fined as  follows. 

(7.4)  Definition:  A  "point"   U0»  f  »  F)  belongs  to  <t>0  if 
it  belongs  to  (ji   and  the  solution  q(s)  of  (7.1), 
(7.2)  is  defined  and  in  R0  for  S]_  *  s  =  33. 

That  is,  the  equations  (7.1),  (7.2)  map  <j)0  uniquely  into 
the  space  Cn  of  n-uples  of  real-valued  functions  defined 
and  continuous  on  the  interval   (s]_,S2),  the  "image"  of 
(qo»  f   >  *")   being   q(.)  S  {  q(s)  |  si  -  s  =  S2  )   .   With 
reasonable  norms  in   (f>0  and  Cn,  it  can  be  shown  that  this 
mappin-  is  continuous  on  <|)0-   In  fact,  this  conclusion 
is  established  in  Theorem  (4.9)  of.  [,..cs]  for  the  norms  N 
and  N0  defined  below. 

(7»5)   Definition:   Usint;  the  summation  convention,  we 
write   |q|  =  (qV)1/2. 
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(7.6)  Definition;   In  the  space  Cn  the  norm  will  be  the 
customary  one,  N(q(.)),  given  by 

N(q(.))  -  sup  {  |q(s)|   |  si  «  s  =  s£  } 

(7.7)  Definition:   In  $Q   the  norm,  which  we  shall 
designate  N0(q0,  if  ,  F)  ,  will  be  the  greatest 
of  the  three  numbers  :  |q0l» 

sup  {  |  if  (s,q)  |  J  si  =  s  «  s2,  q  in  R0}  ;   and 
sup  {  |F(r)|  |  Ir^l  <  B,  <*  »  l,...,m]  .    [By 
|  tf  |  we  mean  ( if "  cf *"  )1/2;   similarly, 
|P|  »  (F1?1)1/2.  ] 

McShane's  continuity  theoremj  on  the  mapping  of  (Jj0  into 
Cn,  is  as  follows: 

(7.8)  Theorem:  When  the  spaces  ty     and  Cn  are  metrized 
in  accordance  with  the  respective  norms  No  and  N, 
the  set  <J>0  is  an  open  subset  of  <j),  and  the 
mapping  of  <{)0  into  Cn  defined  by  the  equations 

(7.1),  (7.8)  is  continuous. 

In  the  opening  paragraph  of  this  section,  it  was 
asserted  that  the  system  (6.3)  is  a  special  case  of  a  more 
general  differential  system.   That  more  general  system 
consists,  as  the  reader  may  already  have  inferred,  of  the 
equations  (7.1)  and  (7.2).   To  see  this,  we  identify  the 
independent  variable  s  with  t  and  the  dependent  variables 

q*  with  x,  y,  z,  vx,  Vy,  vz though  not  always  in  the  order 

named.  The  functions   ^    we  tal?e  to  be:   t,  x,  y,  z,  vx, 
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vyi  ?z>  wx,  wz,  ^  ,  ©,  ax,  ay,  az.   Thus,  for  the  system 
(6.3),  we  have  n  -  6   and  in  »  14.   Given  a  law  of  air  re- 
sistance (i.e.,  a  particular  "drag  function") ,  the  function 
B  can  be  calculated  from  the  values  of  the  14  funotions 
f  **  ,  and  hence  the  right  members  of  (6.3)  are  the  func- 
tions F*  with  the  r1  replaced  by  the  14  functions   ij?*  , 
which  in  turn  are  functions  of  the  6  dependent  variables  q* 
and  the  independent  variable  t. 

Theoretically  there  is  a  one-to-one  correspondence 
between  projectile-types  and  dreg  functions.   But  the  de- 
termination and  tabulation  of  a  drag  function  is  an  expen- 
sive process;  not  only  are  costly  experiments  (firing  pro- 
grams) necessary,  but  extensive  calculations  are  required 
to  reduce  the  data  obtained  from  the  experiments.  For- 
tunately, the  function  E  has  as  a  factor  not  only  the  drag 
function  but  also  the  reciprocal  of  an  empirical  quantity 
called  the  ballistic  coefficient ;  by  altering  the  value  of 
the  ballistic  coefficient  when  changing  from  one  projectile- 
type  to  another,  sufficiently  similar  to  the  first,  ballis- 
ticians  have  found  it  possible  to  use  the  same  drag  function 
for  more  than  one  projectile-type.   The  functions  F*  are 
always  the  same  unless  the  drag  function  is  changed;  since 
it  is  rarely  desired  to  calculate  even  first  differential 
effects  due  to  changes  in  drag  function  alone,  we  have  not 
felt  it  worthwhile  to  complicate  the  notation  in  our  theory 
of  second  differential  effects  by  including  terms  which 
arise  from  changing  the  functions  F*. 
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Accordingly,  the  only  "disturbances'*  which  we 
shall  consider  in  the  remaining  chapters  are  those  which 
consist  of  changes  in  initial  conditions,  or  of  changes  in 
the  functions  t^**  ,  or  of  both.   A  disturbance  which  con- 
sists of  changes  both  in  the  initial  conditions  and  in  the 
functions  if  we  shall  denote  by  the  symbol  (  A   qo*  A*p  ). 
In  accordance  with  this  notation,  (  A  qo»  0)  will  denote  a 
disturbance  whioh  consists  of  changes  in  initial  conditions 
only;  (  0,  A«f  )  will  mean  a  disturbance  which  consists  only 
of  changes  in  the  functions  ip   •   Therefore  we  shall 
suppress  the  symbols  F  and  A  F  from  such  notations  as 
(qo,  4?  ,  F)  or  N0(  A  qp,  A«p .,  A  F). 


20 


Chapter  III 

THE  (FIRST)  VARIATION  EqUATIONS.   CONTINUITY 
AND  DIFFERENTIABILITY  THEOREMS  OF  BLISS . 

8.  The  (First)  Variation  Equations.' 

At  the  beginning  of  §   7  we  made  certain  assump- 
tions regarding  the  continuity  of  the  functions 

(«"l,...,m)  and  F1  (i  m   l,...,n)«,  For  our  discus- 
sion of  first  and  second  differentials  of  the  mapping  de- 
fined by  (7.1)  and  (7.2),  we  shall  need  additional  assump- 
tions concerning  the  differentiability  of  the  functions 

and  F*.   These  additional  properties  are  described 
in  the  following 

(8.1)  Definition:  Let  (^a8  be  that  subset  of  (j)0  which 
has  the  properties 

(8.1i)  the  functions  F*  are  twice  continuously  differ- 
entiate on  the  set  of  r  with 
|r  **  j  »  b,  <*  *  l,...,iij  and 

(8.1ii)  for  all  s  in  lsi,S2)  and  all  q  in  R0,  the  functions 
tf '   have  bounded  partial  derivatives  of  the  first 
and  second  orders  with  respect  to  the  qi,  continuous 
except  perhaps  at  a  finite  number  of  values  of  s. 

-  should    pointed  out  here  that  the  assumptions  con- 
cerning sec  one;  derivatives  contained  in  (8.1i)  and  (8.1ii) 
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are  unnecessarily  restrictive  for  the  theory  of  first  dif- 
ferential effects.   In  fact,  McShane  uses  in  [mcs]  the 
subset  <t>3  of  4>0  defined  by  altering  definition  (8.1)  in 
such  a  way  that  the  word  •twice*  is  omitted  from  property 
(8.1i)  and  the  words  "first  and  second  orders"  in  property 
(8.1iil  are  replaced  by  the  words  "first  order".   Since 
evidently  ($s8  c  $St   the  results  in  [mcs]  are  valid  on 
the  set  <frss.  ] 

We  wish  to  emphasize  that  in  the  definition  (8,1), 
as  in  the  analogous  definition  for  ()s  in  [Mcs]  ,  the  var- 
iable s  chosen  as  independent  variable  is  distinguished 
from  all  the  others,  i.e.,  from  the  variables  q*.   For  in 
(8.1ii)  the  partial  derivatives  of  the  functions  <f  **  with 
respect  to  b   are  not  mentioned.   This  is  the  reason  for 
using  the  letter  s  in  the  subscripts  on  the  letter  <). 
For  example,  the  set  0tt  (i.e.,  $8s  for  s  -  t)  is  not  the 
same  set  as  <)yy  (the  set  4ss  when  y  is  chosen  as  inde- 
pendent variable ) • 

Now,  let  (qo,  <f  )  and  (q0,  if  )  be  points  of 
4SS,  and  let 

*1  1o        *      \~    X 
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}for  each  real  number  f  ^  0  •     f    *  1»   let 


H 


(S,   T    )    I       (q1     (S,«),     .     .      .     ,    ^Q     (S,f  J     ) 


be  the  solution  of  the  corresponding  differential  system 


J!  S 


(8.3) 


?V*i>*>  5  t. 


f  A 


i 


i.e.,  we  assume  that,  for  each  f  ,  0  ■   ^  ■  1, 


:-'  9 


(8.4) 


Let 


<rc*t,r)  s  V  +  *  Ai^ 


tS)        "2 


a?  cs>°) 


(.s.sr) 


v 


C*> 


and,    for  brevity,    let 
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^M=  F^  (t?(s,Q<s,*))   +  *   A±<f  Cs,  Q<*rJ)J 


£.6) 


o( 


with  similar  abbreviations  for  F  «.  a    »   *P  :  ^.  k 
A^  cp  ,  A^  tp  j   ,  etc.   In  this  notation,  the  result 


of  differentiating  each  member  of  (4)  with  respect  to  X   » 
and  then  interchanging  the  order  of  differentiation,  may  be 
written 


f  TW  9rfc*,rj  ^f  1tH*Axy  jW  Q^  c*.*) 


Lt 


f*.7) 


Thus,  the  functions 

(8.8)      iftt,  S  i?tc*'i^#,.«f;A1^..M^si?rf*'#; 


satisfy 


i3l  =  fMoJ 


(8.9) 


ds 


1r 


»|£  <"*«.»    5    A±  V 


s 
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which  is  the  same  as  the  system  (6.1i),  io.lii)  of  [tics] 
with  A  replaced  by  A  i  and  with  A?  »  d  [a?  »  0 
because  we  ;  ,'re  U3ed  the  seme  functions  P*  in  both  point 

(<lo>  <P  »V)*    do*   "(?$$)•       In  th9   closing  paragraphs  of 
|   ?  we  explained  why  this  was  done.] 

?he  equations    (3.9)   have   come   to  be  known  in  the 
theory  of  differential  equations  as  the   "eolations  of  varla- 
tioa*«*       In  this  paper  we  shall  refer  to   -Ghent  as  the  "firsx. 
variation  equations",    to  distinguish  them  from  another  sys- 
tem (which  we  a hall  call  the   "second  variation  equations") — 
to  be   deducec;  in  Chapter  IV.        The  functions  I?  Ms)  which 
constitute  the  solution  of    (8.9)   are  cal-ec  the  "variation*, 
or,    as  we  shall  refer  to    them,   the   "firs t  variations".** 
In  fact,   \,Le    .:'irst)    differential  effects  on  trajectory 
elements  &rs   certain  linear  combinations  of  the    (first) 
variations;;  4     the   first  variations  also  tre  prominent  in 
the  theory  of  second  differential  effects,   as  we  shall  see 
in  the  next  chapter.        In  applying  the  theory  of    (firs*,, 
differential  effects,   it  is  seldom,   if  ever,  necessai- 
solve  the   equations  of  variation    (8.9).       A.  simpler  pro- 
cedure   (also  due  to   Bliss)    is  to   solve   a  related  homogeneous 
differential   system,    said   to  be   "adjoint"  to    (3.9),   and 


'  See    [5],   p.   59;    [6],    Chapter  VIII;  and    [t], Chapter  VIII. 

**See  also    [9] ,   p.   165. 

■  See   [5],   chapter  V;    [e],   Chapter  VIII;  and  [?], 
Chapter  VIII. 
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then  Dtain  the  de  differential  effect    (linear  com- 

bination of  the   function        f£  *)    by  means   of   a  fairly  simple 

[2]    between  the    solutions   of  the 
system    (8.9)    and   those   of    its  adjoint   system. 

9.     Linearity   and   Continuity  of  the   First  Variations. 

i  we    shall  state   three  results   from 
>]   which  we  shall  need   in   the   next   chapter. 

As   observed  in     §   6  of    [lies] ,    the   solutions   of 
(8.9)    are   obviously  linear  in  the   following  sense: 

(9.1)  If     Tp  1   corresponds  to    (  A  1qQt  Aj_tf  ) 
and     1^2  corresponds    to    (  Ag  q0>   A  gl^  )(    and   if  a]_ 

a2  are  real  numbers,    then     ai  *?  1  ♦   a2    ^?  2 
corresponds   to    (ai    A  iq0  +   a2  &2lo» 

al  ^lf     +   a2   ^2^)0 

Moreover,    the    Y~l       are  continuous  functions  of 
the   disturbance    (  A  q0,   A  tp  )    in  -the  sense  of  the   following 

(9.2)  Theorem:      There  exists  a  constant  K  such  that 
N(  »?(.))    -     m0(  Aq0>  Af ). 

The  proof  of  this  latter  theorem,  as  given  in 
[lies] ,  depends  upon  the  following  lemma  on  differential 
inequalities. 
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(3.1)    of    [Mcsl  J    :     Let    A(s)    be  a  bounded 
Riemann  integrable    (or  Lebesgue  summable)    function, 
on  an  interval    (31,32),  and  let   B(s)    be  bounded  and 
integrable    (or  Lebesgue   summable).        If  f(s)    is  a 
function  Which  satisfies   the   inequality 

(9.4)  f'(s)      -     A(s)f(s)      ♦      |B(s)| 

at   almost   all  values  of  s   such   that      f (s)   >  0,    it  also 
satisfies 


(9.5)     $cs>    =    ex 


f,$*AU)b{  l-f Ch)  l+  /,B  C*;l l**tS/'A(t))lr]lt}  , 


\ 


We  shall  have  occasion  to  use  this  lemma  in  our 
proof  of  the  principal  theorem  in  Chapter  VII.   Mcshane 
also  employed  it  in  his  proofs  of  the  two  theorems  of  Bliss 
which  are  referred  to  in  the  next  two  sections. 

10.   Continuity  Theorem  of  Bliss. 

An  important  theorem  of  Bliss  (  [1],  p.  85),  con- 
cerning vihat  he  describes  as  "continuity  of  Type  I*,  fur- 
nishes us  with  an  estimate  of  the  magnitude  of  the  actual 

Tect  of  a  disturbance.  We  state  it  below  in  the  form  in 
which  it  is  proved  in  [Mcs] .    (See  (5.8)  of  that  paper.) 

Theorem  (Bliss):  Let  (q0,  <f  )  be  a  point  of  (}s 
and  q(.)  the  corresponding  solution  of  (7.1),  (7.2). 
Then  there  is  a  constant  Kl  such  that  for  each 
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(<lo»     *f  )    in     ^o     the  corresponding  solution  q(.) 
of    (7.1),    (7.2)      satisfies  the   inequality 
N(    q(.)    -   q(.))      ^KiN0(q0  -   qot     «f  -  *f  )  • 

11.      Definition  of  First  Differential.        A  New  Norm. 
Bliss'    Theorem  on   the   First  Differential. 

As   intimated  in  the   last  paragraph  of      §    8,    one 
does  not  have    to   solve    the  first  variation  equations    (8.9) 
in  order  to   obtain  the    first  differential  effects   on  tra- 
jectory elements.        The    system  actually  to  be   solved  is 
either  the   system  adjoint  to    (8.9)   or    (see   Chapter  VII) 
the  homogeneous   system  obtained  from    (8.9)    by  discarding 
its  non-homogeneous   terms.        Nevertheless,  given  the  equa- 
tions of  motion    (6.3),    one  finds    it  necessary  at   least  to 
write   the   corresponding  first  variation  equations  before 
either  the  adjoint   system  or  the  corresponding  homogeneous 
system  can  be  written.        This   can  obviously  be  done  by 
calculating  the   partial  derivatives  F^      and      4*    j     and 
then  substituting  the   results,  together  with  the   disturb- 
ance functions      A±  *f      in  the  right  members  of  the  equa- 
tions   (8.9). 

But  there   is   still  another  way  in  which  the  first 
variation  equations   can  be   deduced  from  the   system    (6.3). 
Let     x  *  i(t),    y  *  y(t),   z  ■   z(t)      be   the   parametric  equa- 
tions of  the  normal  trajectory  defined  by    (6.3)   with  normal 
initial   conditions   and  standard  functions     ^>  ,   ©,   wx  -  0, 
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wz  -  0,  ax  -  0,  ay  ■  0,  az  «  0.   Let  x  -  x(t)  *     A  x(t)  , 
y  »  y(t)  +  Ay(t),  z  -  z(t)  «■  A  z(t)   be  the  parametric 
equations  of  a  disturbed  trajectory  defined  by  (6.3)  with 
non-standard  initial  conditions  and/or  abnormal  functions 

^  +  L  ^    ,  9  ♦  A  ©,  wx,  wz,  ax,  ay,  az.   Since  these  lat- 
ter functions  and  E  are,  by  the  hypotheses  stated  in  $  6, 
functions  of  x,  y,  z,  vx,  vy,  vz,  so  are  the  right  members 
of  (6.3).   Let  us  expand  the  right  members  of  that  system 
as  power  series  in  A  x(t) ,   A  y(t) ,   A  z(t) ,   A  vx(t) , 

A.    Vy(t),  A  vz(t),  assumed  to  be  convergent  for  suffi- 
ciently small  values  of  these  latter  differences.   If  we 
retain  only  the  linear  terms. of  these  series,  and  if  from 
the  results  of  this  procedure  we  subtract  the  differential 
equations  of  the  normal  trajectory,  we  obtain  the  first 
variation  equations  with  /\x(t),  .  .  .  ,  A  vz(t)  in  the 
roles  of  the  variables   1? fc  which  appear  in  (8.9).   (See 

[10],  Chapter  IV;  [ll] ,  Chapter  II;  and  [7],  Chapter  VIII.) 
This  is  the  method  usually  employed  in  elementary  exposi- 
tions of  the  subject  of  approximations  to  effects  of  dis- 
turbances. It  has  certain  undeniable  pedagogical  advan- 
tages, and  undoubtedly  the  solutions  of  the  first  variation 
equations  are  "approximations"—- of  some  kind— to  the  dif- 
ferences  Ax(t),  .  .  .,   Avz.(t), 

Neither  of  the  two  procedures  described  above  for 
deriving  the  variation  equations,  however,  throws  any  light 
upon  the  important  question  of  the  accuracy  of  the  approxi- 
mations.  McShane,  Kelley,  and  Reno  (  [6],  Chapter  VIII) 
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call  attention  to  the  analogy  between  the  expansion  of  our 
equations,  to  terms  linear  in  the  disturbances,  and  the 
directional  derivatives  of,  say,  a  function  of  two  real 
variables:   it  might  be  thought  that  if  such  a  function  has 
a  directional  derivative  at  the  origin  which  is  zero,  re- 
gardless of  the  direction,  then  that  function  must  be 
closely  approximated  by  a  constant  on  some  region  about  the 
origin.   But,  as  is  well  known,  there  exist  functions  with 
such  directional  derivatives  which  are  discontinuous  at  the 
origin.   Bliss  appears  to  have  been  the  first  to  recognize 
the  need  for  an  estimate  for  the  error  involved  in  using 
the  variations  instead  of  the  effects;  i.e.,  in  the 
s  -  q  _  tf    notation,  an  estimate  of  the  magnitude  of  the 
quantities  |  Aq^s)  -   J?*(s)  |.   He  obtained  an  upper 
bound  on  the  error  by  showing  [see  (11.3)  below]  that  the 
mapping  of  the  space  4>o  into  the  space  Cn»  defined  by 
(7.1)  and  (7.2),  has  a  first  differential  in  a  sense  analo- 
gous to  that  defined  by  Frechet  in  the  theory  of  functions 
of  lines. 

In  order  to  clarify  the  meaning  of  Bliss1  theo- 
rem, we  state  next  a  definition  of  a  differential  of  a  map- 
ping such  as  that  just  referred  to.   This  definition  para- 
phrases, in  the  notation  of  the  present  paper,  one  given  by 
Hildebrandt  and  Graves  [8],  who,  in  turn,  assign  credit  to 
Stolz,  Young,  and  Frechet. 
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(11.1)      Definition:      Let     3~  ($)    denote   the    mapping  of      (J>0 
into   C     defined  by    (7.1),    (7.2)>let     lj[  =    (q0,  <f   ) 
be   a  point  of      ()>o,    and  let  N'(lj0    be   e   norm  defined 
in  <j)0.        Then    &    is  said  to  have  a  differential 
relative    to  N'    at    the  point  Tjr  if   there  exists   a 
function  dj?-(  A  TJT)    on  (J»   into  C      linear   and   continu- 
ous   [in  the    sense   of    (9.1)    and    (9.2)]    on   4>»    such 
that   the   function  R  on  ^o  to  Cn  defined   for 
N»  (  A  Tjj)    t  0  by 

(11. li)       |  9  (ijT  +    A  Tfr)    -    3-  (T|r)    -  d   *  (  A  TJF)  I 

•R(t^t)    N»(  A  ifr) 
satisfies 

(ll.lii)  lim  R(Tjr  +  *  t)   =0. 

H»  (  A  ifr)  — ^0 

11.     A  New  Norm.      Bliss*   Theorem  on  the    First   Differential. 

In  order  to   state   the    theorem  of   Bliss  on   the 
first   differential  of  the  mapping  defined   by    (7.1),    (7.2), 

shall  need  to  define   a  new  norm.      This   new  norm  is   the 
one   denoted  by  the   symbol  Ng(q0,  if  ,   F)    in    [mcs] ,    but  we 
shall  designate    it  as  Nis(q0,  if  )    in  order  to  distinguish 
Lt  from  still  another  norm. which  we   shall   introduce  In  the 

chapter.       Moreover,    when  there   is   no  necessity  of 
emphasizing  its   arguments,    we  shall   refer  to    the  new  norm 
as  Ns  or  as  Nis.        It  is  defined  as  follows,    (see     #  6 
of    [mcs].) 
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(11.2)   Definition:  Let  N]_s(q0,  <f  )  denote  the  greatest  of 
the  three  numbers: 

(ii.2i)  N0(q0»  ¥  ); 

(11.211)  sup  |  <fqj  (s,q)|   (j»l,...,n;  si  -  s  =  s2;  q  in  R0) ; 
(11.8111)  sup  |F^a  (r)|     (  «*-l,...fiB;   |r*|  -  b)  . 

It  should  be  emphasized  that  the  definition  of  the 

norm  N15  does  not  involve  the  partial  derivatives  of  the 

functions  if      with  respect  to  the  independent  variable  s; 

the  derivatives  of  the  £P   with  respect  to  the  dependent 

variables  qJ  are  the  ones  which  appear  in  (11.2ii),   It 

is  the  omission  of  the  partial  derivatives  *t   in  "the 

s 

definition  (11.2)  which  causes  the  norm  N]_y(Nig  when 
s  =  y)  to  reduce  to  the  more  reasonable  norm  No  in  the  case 
of  a  one-branched  trajectory  (on  which  the  choice  of  y  as 
independent  variable  is  possible)  along  which  the  disturb- 
ances being  considered  are  assumed  to  be  dependent  on  alti- 
tude alone.   Such  disturbances  include  the  important  ones 
due  to  winds,  abnormal  density,  and  non-standard  tempera- 
ture.  (See  §   2).   For  two-branched  trajectories  we  cannot 
use  y  as  independent  variable.   Therefore  it  is  by  no  means 
obvious  that  on  such  trajectories  the  differential  exists 
when  Niy  is  used  as  norm.   This,  however,1  is  established 
in  the  main  theorem    [(11.7)3  of  [Mcsl   • 
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In   terras   of    t  is   norm  Nis  we   now    state   Bliss' 
theorem   on   the   first  differential    of    the   mapping   defined 
by    (7.1),    ( 7 .  £ )  .        Thi  s    is  t  he  th  eorera    (0.5)    of    [[..  cs]  . 

(11.3)      Thee  rem    (Blirs)  ;      Let    (q0,  <f  )    be    en   eletrert    of  (j>£ 
such    that  the    i«'i      have   continuous  partial    deriva- 
tives   of   first   and    second   order  with  respect    bo   the 
r  *,    and    the     ^       have   bounded  first   and    second 
order  partial    derivatives  with   respect  to   the      qJ, 
continuous   except  perhaps   at  a    finite    number  of 
values   of    s.        Then   there   exists   a    constant  Kg   such 
that  whenever    (qo   +  &  qo»    if   +    &■  (f  )    is  in  <|>s  end 
|0   ( . )   =     {  Y)   s)    |    si  ■  s   *  sg  }      is  the  function 
defined   by    (8.9),    the   inequality 

(11.31)    N(   17   (.)    -   Aq(.)     )i     K2K0NIS   -  Kg[Nls(  AqD^f)]2 
is    satisfied. 

With  respect   to    the  norm  N]«i    the    first   inequality 
(11, 3i)    is   a  stronger   conclusion    than  the    statement   that 
Y^    (.)    is   the   differential   of    o(.).        The   second  inequality 
in    (?.l„3i)    is  obviously  e  consequence  of    (11„2),    sine?'   from 
the",    definition  we   immediately  have     N0  »  N]». 


Chapter   IV 

ON   THE   EXISTING ii  Or    flffi   SjSWwD  DIFFERENTIAL, 


THE   oECO..l>    VA^IATIUK    EQUATIONS 


1 2 ,      Yeaning   of   the  Terra    "Second  Differential". 
A  New  Nomi. 

In   the    present    chapter,    it  will  be   our   objective 
to    show    the   existence    of   the    second    differential   of      q(o)P 
i.e.,    of  the    (first)    differential  of     f}   (.),    at    (q3,  f   ). 
But    this   latter   ter^  must  he   more    clearly  defined   before 
we    can  proceed   further.        Let   us    first  assume    that 
(qo>    *f  )    belongs  to     <$ss  ,    and  that    (q0   +  Aiq.0»  f  +  ^VP  ) 
belongs   to      <p      .     Let   q(s}    be   the    solution   of 


(12.1) 

S     is 

Fl 

(  \'t*i) 

and  let 

(12.2)       •(    AzfCf,<p  =     If  c*.  V    -    fC  V 
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Finally,  let  V?  . )  be  the  differential  of  q(.)  at 
(qo,  $   )  relative  to  the  norm  Nis(A1  q0l  A  i^P  )    i.e., 

(IE. .-5)   V?  I.)  =   ^  »2  (s)  I  si  »  a  -  s2  }    , 

in  which 

(12.4)  7}   (3)  -  r?  Is;  q0,  «f  5  Al  <lo»  Al  'f  ) 
satisfies  the  system 


(12.5) 


d 


.*  ~ 


a  -  -  FJL  ( f C5'  $w))  ?  J  c*,y^?J+ AJ^v^ 


l  v 


f?l<«4>'*   4.V 


obtained  from   (8.9)    by  replacing     (^  (s,    q    (s)    )      with 

^   (s,   q(s)    )«       It  is  our  desire  to  find  a  function  which 

approximates  to  the   difference 

(12*6)  A   ^  (s)    a     ^    (s)    -     17   (s) 

in  a  manner  (to  be  more  carefully  described  later)  such 
that  it  may  appropriately  be  called  the  differential  of 
1^  (.),  or  the  second  differential  of  q(.).   But,  as  we 
have  seen  in  our  review  of  the  theory  of  the  first  differ- 
ential, the  definition  of  a  differential  requires  the  prior 
specif icetion  of  the  norm  of  the  disturbance.   As  might  be 
expected,  we  find  it  necessary  to  introduce  a  new  norm— in 
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general  larger  than  Ni3— in   or<Inr  to  establish  the    exist- 
ence of  a  second  differential    of  the  mapping  defined  by 
(7.1) ,(7.2).        This   new  norm  we    shall  denote   by  Nga    (q0,  <f   ) 
and  define  as  follows: 

(13.7)  N23    (q.o»  H  J    "  tfte   greater  of  the  numbers: 

(i)      Nia    (qot   f  )J 
(it)    ^p  {    |<fqi    q*    (s,q)     1}  , 
(J,k  -  l,...,n;   si  -   s  «  32   ; 
q  In  Rq). 

Y«e  observe  at   once  that  Ngg  share;   a;i  important 
property  with  the  norm  M^g  -—on   trajectories   alop.2;  which  it 
is  possible  to  use   altitude  y  as   independent  variable  of  dis- 
turbances due  to  wind,    abnormal  temperature,  and         lortHal 
density,    the   norm  N2s  reduces  to  the   simpler,  more  natural 
norm  N0.       And  now  we   can  formal:,.,  e    our  definition  of  the 
second  differential  of  q(.),   relative   to  the  norm  N2S* 

(12.8)  Definition:      For  fixed  Ai  tf,   the   first  differential 
d  9-  (ijr,  A  1  XJT >   depends  on     f     and,    relative  to  some 
norm  N"(l)f)    defined  in  <t>0i  may  have  a  differential 
which  is  a  function  of     TJT     and     A  2  $»   linear  and 
continuous    [as  in    (9.1),    (9.2)   ]    in     A  2  $•        *f  so> 
we  denote   it  by  the   symbol     d2  -9-  (V»^l  $»  ^  2  10 
and  call  it  the   second  differential  of    d"  at     if 
relative   to  the  norm  N*. 
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From  the  de  fin  it  ion  (11.1)  of  a  firs  6  differen- 
tial, the:  function  R»(t,  A  i  $.  ^  2  fS  defined  by 

|d?(?  +A2$,Aitf)   -  d5rt,/l|t]    -   d2^  (IJF,  Ait,A2<r)| 

-   Rt(ijr;  A^,  A2f)   N«(A2f) 

tends  to  zero  as  N*(  ^2$)    does.         «re  ahall  show  more  than 
this.       It  turns  out  that  d2  9~     is  bilinear  and  symmetric 
in  A  i  |,    ^  2  1*»   ttnd  satisfies 

l^ttjr  ♦  A  *)  -   -sUf)  -  d  2  Ofr,  A    .) 

-  £  d2  3-(#,  A  *,  4  *)|    -  R«tftA*)   N2aUt)   Ng8Ul» 
wherein  Rn{$»  A  ifr)   tends  to  zero  with    A  $• 

13,     X  yirst  Approximation  to  the  Second  Differential, 

Since  our  objective  is  to  show  the  existence  of 
the  second  differential  of  the  mapping   (7.1),    (7.2)   rela- 
tive to  the  norm  N28»  w©  begin  in  this  section  our  search 
for  the  function      <L*  3~     referred  tc    in  the  definition 
(12.8).       Our  quest         is  not  immediately  successful;  we 
find,    instead,    certain  functions  which  -tie  designate  T  *  ( 
and  which  are  almost  the  ones  we   desire.       They  lack,   how- 
ever,   the  property  of  bi linearity  in  the  disturbances, 
the  other  hand,   this  observation  enables  us  to  alter  cer- 
tain steps  in  the  argument  of  the  present  section  in  such 
a  manner  that,    in  the   section  whioh  follows,   we  are  able   to 
obtain  the   functions     y^s)   which  we  seek. 
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In  order  to  make  use  of  results  in  Jkes],  we 
shall  find  it  convenient  to  write  (8,9)  in  the  form  which 
(8.4)  assumes  when  t  -  0;  i.e.,  in  the  form 


(13,1) 


1}   f  l£  -*1  <<f  (.,  q)  ) 


<4  («i)  -  qj  , 


the  system  satisfied  by     q(s).       With  this  purpose  in  mini 
we  define,   for  i  -  1,   2,    ,  n,     the  functions 

(13.fi)     O1   (p1,...,?*  )   e  p1   , 

(13.3)  ^V^i^CT^^fc^V^)^*   \<?*W>]    , 

(13.4)  ift&I   ^/^V^>)[fev«)^4   A<«p"c**«o]  f 


>5  (s,i?)  -  (\51(b,»?) t9n(s,  7  )   )  , 

(13.6)  $   (8,7)  -  (  \3  1(str^  ),  ....      £n(e,v?  )   ), 

(13.7)  Atf  (s.ijJb    9  (■•  17  )  -  19  (•»  17  )• 

Then   (8.9)   end   (12.5)  become,   respectively, 

U3.3,        <      ^     -     "»".»)» 
H1    (81)       SA^ 
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and 


If  we  write 

N0(  4*  )   »  No(0,  A  If  )     % 
N8(    **)    «  Ni8(0,4  *    I  , 

end  let 

be  the  function  which  satisfies  the   system 

then,  by   (11.8)   —since  the  functions  G*,   T&       clearly 
satisfy  the  hypotheses  of  that  theorem— it  follows  that 

(13.11)     N(  A»?(.)-  f  (.))-  0  (N0(  A*  )  NB(  **)    ) 

But  the   "disturbance"    A1?   depends  upon  the 
"disturbances"      (Ax  qo,   A  ±Lp  ),    (  A2q0,  A  2  If  ) , 
Hence  it  is   to  the  norms  of  these  latter  "disturbances" 
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that  we  wish  to  compere  the  left  member  or  (13.il).   Thla 
requires  that  we  express  A^  in  terms  of  functions  whose 
magnitudes,  relative  to  these  latter  norms,  we  can  estimate 
by  previously  established  results.   The  first  step  in  thia 
process  is  to  observe  that 

(13.12)  A*  t*,i>J  =  V\^Jf('4c*>>)$*,(sAiti} 

To  evaluate  the  right  member  of  (13.12 • ,  we  shall 
find  it  convenif)nt  to  use  the  abbreviation*) 


and  similar  ones  for  the  partial  derivatives  F  ^  /$     » 
({>"  »    &  *  W    •      t  A*  <f  .  '   etc.   In  this  notation, 
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(13.U)  Vrt»,q(e))  ■   Atf[o]  ♦  {^f'ni-  \f[o}} 


'•  Y 


<13»16)   F^(f(a,${8)))   *  ?;  [0]   ♦   {  T\  [1]    -  F^    [0]} 

♦  A^  (e,q(o))J      f'   „Wat 
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and 


(13.16/ 


'*#-  -/ 


,<*,     ~, 


I*,-*. 


dt 


Substituting   (13.14),    (13.15),  and   (13.16)   In    (13.18) 
yields 


4£ 


;ia.i7)  4eN.,.?).tfJd,qW>]/<f«    w  «* 

*  F\, [0]  A  <f "  [0] 

+  A.i  (•)*■;« mj  ^ .   w  a* 

*  a,qK<  •)  \«(£  M  f  *  r^ „  M  (f£  W  a  * 

♦  ( *tq1  a) )  •/  F^  W  <$M  a  */^„  W  a* 

+  ^«fflro]fJ  [o]  /  »*.  [r]  ar 

♦A  ffcU  q(»)  J  *i  *W**fAf  ?  Wdt 
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♦<MK<*)<f"ro]  /  £„[•]««/  aa  »  I  Hat 

•  Off  *        )| 

♦A^qW.f  ?M  J  *'  ,  Hd*J  ^f  *W  d* 

♦  fl,  q*i  »)f;«  lo]  ^  d1<f  "   W  c-r 
-dxqK(»)  Atf"M  J^rC  M  f^K  W  a* 
♦Vfc^qW*:, >Mf^C4d«J  d4f*H  a* 

o  •  •  T 

♦  ^f  *[o]  A,f  >]  J  f  i  .„  K  a-r 

♦  4xq  <•)*..  f  "[o]  J  f1    ^  W  a*/  *4  v  i  W  a  * 
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6fc  qK(s)  A4(fM  J'V^^Hdr  J\  f'>]d 


Since  wo  assume  that     A  F  ■  0,  and  alnce  we  ahall  naed  to 
distinguish  between     N^gfA  i  qoi^l^)      *Qd 
1*18^2^01  ^2  H    )•  we   sna*!  tint  it   helpful  to  abbreviate 
these   last  named  noras   by     $\&[  A  l)   and    Nia (  A  g),   re- 
spectively,   as  well  as   to   employ  similar  abbreviations  Tor 
other  norm: 

There  ere   twenty-two   terms   In  the  ritfht  member 
(]     •  I?  , »        To  avoid  writing    each  one   o$aln,    let  us  da- 
note   J  -   after  replacing     V}  *  by     ^  *    (s)   —  by 

Ti,    •    .    •    ,  respectively,   the   subscript  indicating 

tho  position  of      <o    term  in   the   order  of   its  occurrence 
in  the    ri    ht   m<         r  of    (13,17). 

recoil:      f i r r- 1 .    thet  by    (10.1}   we  have 
(13.16)      &  2   q*    (S)    .   Q   (i;0    (  *  g)    )    ; 

secona,    that  by  the  continuity  of  the    functions    19^  (s) 
[see    (9.2)],   we  have 


(13.19)    9  i  (s>  =  ^  (N°  (  A  i)  )  ; 

and,    third,    that    the    first   anc.    second   partial    derivatives 
of  aid     F1   are,    by    (3.1),    bounded.         .'cnce  v/e   conclur]e 

that 


(1C .20.1)  Ti        =  0 


113.20. 2)  T2        =   0 


1 15.20.3)  T3 


(13.20.4)  T4        -   0 


(13.20  .5)  T5        =   L 


,13.20.6)  T6        =   <J 


(13.20.7)  T7 


(13.20.8)  'C8 


(15    -J. 9)  To 


(13.20.10)         T10      ■    ° 


(13.20.11)         T11      =    ° 


(13.2J.12)         Tig      »   0 


(13.20.13)         T13      -   0 


(13.20.14)         T14      *    0 


(13.20.15)         Ti5      -   0' 
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N0    ( 
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N0    ( 
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A   2)    K28    <   A  2)    ) 


^  2)    ') 


A2)2    ) 


A  2)   Nl8(  A  ,  i 
4   2)2  N2s(^2)    ) 

A2)    ) 


a  ?.r  ) 


*  2)  Nl8i  ^  2)  •• 

^   2)2   N2s(A  2)    ) 
A  g)      Nls(A2)    , 

A   2)2  Nls(  A2)    ) 
A  2)      Nls(A?J£    ) 
^2)2  Ni8(A2)    N2si  *-2;     ) 
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(13.P0.16)    T 


16 


0(Nf 


(13.20„17;    T17   =  0(N, 


(13.20.21) 


(13.20,22.) 


'o 


(13.20.13)    T18   =  0(No 
(13.20.19)    T19   -  0(N0 


(13.20.20)    T20   =  0(No 


T22 


=  0(No 


=  0(No 


Consequently,  from  (1Z.17 
(13.20.22),  inclusive,  we 

N0  (A  ,5  )  *  £  (Nig 


A2)  Nls(2l1)  ) 
A2)  N0  (^x)  ) 

A2)2  Nl8U!)  ) 

*l)  No  (42)  ) 

A2)2  Nls(Ai)  ) 

A2)  N0  (  AX)  Nls(^2)  ) 

^2)2  Nl8-(^l)  -ls(A2)  ) 

and  the  equations  (13.80.1)  to 
see  that 

Al>  N2s(A2)  ) 


and 

tfl8 (  A?)    *  Ns(  A*  )    =    0    (N2s(A2)    ). 
Hence    (13.11)    becomes 

(  A  9  (.)    -  J  (.)     )    «  fl   (Nls(Ai)    N2s(^2)2    )• 

v;e  restrict  our  attention  to  the  case  in 

which  A  £f  =S/A^l(f»  and  the  ra'tio  /^ 
is  independent  of   q^   and  also  satisfies 

^2  1o  =A  Al  q0) 


e's  we  shall  later  wish  to  do,  then 
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so 


N0(  A  l9  )    ~  ^   0  (Nis(  AX)    M2b(  Al)     )• 
^ls(  A*  )    =/u^(N2s(^l)    ), 


and 


(13.22)      N    (A!?(.)    -    r  (•))    "  /^2^(Ni8(  Ax)    N2s(  Ax)2 


Now,   because   of    (13021)    or    (13.22),   we  might  be 
led  to    conjecture  that   the  function    3"  ( . )    defined  by- 
is  the   function  we   are   looking;  for,    viz.,    the   differentia 
of    Yp  ( . )  .        But   a   moment's   reflection   is   sufficient   to   cor 
vince   us   that  this   is  not   the    case.        In  fact,    J  (.)      ack 
one   of   the   essential   properties   of   a   differential :      it 
fails    to  be   linear  in  the   disturbance    (A  2  Q.o»^2^) 
which   the   difference      A  1?  ( • )    is    attributable.        To    see 
this  we  need  only  to   observe, from  equations    (13.20.1) 
through    (13.20.22)    inclusive,    that    Al?(s,    9(3))      is 
linear  in  ,(  A  2  <lo>  ^  2  *f  ) »   an(i»    consequently,    that 
is   also  not   linear   in  this    disturbance.        On  the  other 
hand,      A  l5  (s,   r?  (s))    is^  linear  in  the    disturbance 
(Al   qQ»    A  1  cp   ).        L.oreover,    since  both     Yf  ( . )    and 

) 
are    (first)    differentials   of   q(.)    due   to   the  disturbanc 

(^1   Q.o»    ^  1  ^f   )    a^d    hence   are  linear  in  that   disturbs- 
it    follows   that     A  \T}  ( . )    =    f?    (.)    -    »?    (.)    is   at  most    . 
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in    (Ai    q0,   ^l^'   )•        These   considerations   suggest   that, 
in  our   search  for  the    (first)   differential   of  •?  ( . ) ,   we    re- 
place  the    term       Ai9     (s,  *?    (s))    in    (13.10)    by  an  approxima- 
tion  thereto  which  is   bilinear   in  the    two   disturbances 

(^1  ciD^l  4*    )    end    ^2   9-o>  ^  2  *P  )  °        Such   an  approxima- 
tion is    furnished   by  the    sum  of  those   terms   in 

A  19  ^(s,    i?  (s))    which   pre    linear   in    (  A  2  Qo»    A  £  ^  )  • 
From  equation    (13.17)    we    see    th«t  th^se   terms   are    T^_,       -» 
T4,    T3,    TiQt    T17»    ana   ™19«        Their   sum,    however,    is    still 
not    quite   the   approximation  v/e  desire,    for  two  reasons: 
first,  we  have   no  means  whereby  to  evaluate   the   integrals 
in  these   terms;  and,    second,   ultimately  v/e   shall  need    co 
consider  the   case   for  which  Agqc  «    L  -^  q0  and    A  ^  =    - 1 4 
under  which  circumstances  the   terms   T^,    T4,    Tig»   and   r;x7 
contain  the   very  functions   A 1  q^   (s)   whose   computation  we 
wish  to  avoid: 

14.      The   Gecond  Differential. 


But  the  difficulties  mentioned  in  the   last  para- 
graph of    §•  13  are  not   insurmountable.       We  may  overcome   th 
first,    in  fact,    by  the   simple   device  of  replacing  each   in- 
tegral  in  Tii    T4,    T8,   Tig,   Ti7,   T19     by  its   integrand 
evaluated  at     X  -  0.        For,    by  the  mean   velue  theorem  of 
the    integral   crlculus,    each  such  integral  is   equal  to   i 
integrand   evaluated  at     V  *  0  plus   a   term  or   terms   of   at 
least  the   first   decree   in    the    disturbance    'Ao   q0,   A  g  <f    j# 
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The    second    difficulty,    too,    is  easily  circumvented.        ,/e 
need   only  define 

V?  2    (s)    2.    ^    (s;    q0,  if   ;  A2   q0,    &  2  Cf>    ) 

as  the  solution  of  the  system 

j*  L  * 

so  that  V"p2(s)   is  identical  with  IQ  (s)   when  ^2  2  &1» 
and,  by  the  continuity  of  the  functions  Vj)  2(s)> 


(14.2;   r)  ]     (s)  =  O  (N0(A2)  ). 


2 

Then  we  may  replace  A  2  q(s)  by  *?2  (s)  in  Tl»  r^4»  '^16 » 
T17,  without  disturbing  the  estimates  contained  in  the 
right  members  of  the  corresponding  equations  (13.20.1), 
(13.20.4),  (13.20.16),  and  (13.20.17). 

Let  T'i,  T2,  T4,  T8,  Tig,  T17,  Tig  denote  the  re- 
sults of  revising  the  corresponding  terms  T]_,  .  .  .  , 
in  the  manner  described  in  the  preceding  paragraph. 
the  corresponding  equations  in  the  sequence  (13.20.1), 
.  .  .  ,  (13.20.19)  are  still  valid  when  Ti,  T2,  T4,  TQ, 
T16»  T17»  T19»  are  replaced,  respectively, by  t"i,  T2,  T~4, 
T8,  T16,  T17,  T19.   Let 

d  15  i(s,  9  (a);  A  X(f  9  A£cf»  ) 
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denote   the   sura 

Tl     +   T2     +   T4     +      8     *   T3.6  *   ?17  *    T19    J 
1  *e  •  9 
114.3)      J^Cr5,  r?«->i     A±cp,A^cf; 

+  f; («pc^^;  a  to0*.  fe*"0  i?J  co 

+F "   „Gf  ^V'%,"  (* V*°  *P- k  *'  >"°  ^  ^ 


F  i       (cpC^-0  A4Cf  V.  V">  **¥  W"' 


r'r^ 


Now  let     ^  (s;   A  1  <f   ,  A  2  ^  )  be  the  function 
defined  by  the   system 
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(14.4)\    ds  t 


j£^c._  G-lJ&Cw*>>)  &.&,7**>)54+-l'& -(w*>i*t%K<l) 


L? 


(  y:Ok;  VP,*»<f^«  ?c?±>-i?«i>  s*     , 


and   let   us    compare   the  function      7(s;    &  £.iX>  ,   A  g  cp  )   with 
the    function     f"(s;    ^tf,    A2  P)    defined  by    (13.10).         <* 
,»e   can   estimate  their  difference   by  Theorem  (10,1).        To 
this   end,    let 


H1    (u1,    . 


,    u11)    =   u1, 


yl(Sj  s)  =  Q^(d  <Ts  v*)))  jV/*  V"^rJ+  A  ^  i?«>.  ^<fV\fJ 


end 


—  «. 


4>    v  '  I    ? 


Then    (13.10)    and    (14.4)    become,    respectively,. 


d  5* 

~d*s— 


H1       (-Cp  (s,   j   )    ) 


(14.5) 


X1    (si)      =   0, 


and 


ds 


i c  (  y  (s,  r  )  ) 


(14.6) 


31  (sx)     =  o, 


and    by    (10.1),    we   hove 


But 


(14.8)   \pu<r^i>-  v"^^ 


and,  since  from  (11.3)  it  follows  that 


N(?2(.)  -  A  2q(.))  -  O   (N0(A2)  Ni8(Ag)  ), 


we  have  from  (13.17)  and  (14 , 3)  the  relations 


TX  -  \   -  C   (N0(AX)  N0(  A2J  Nls(A8)  ), 
T4  -  T4  -  Q   (N0(AX)  N0(  A2)  Nl8(Ag)  ), 
*16  -   Tie  =0  (N0(  A2)  Nls(Ag)  Ng8(A1)  ), 


and 


T17  -  T17  =D  (NoiAi)  W0  (Afi)  Nls(A2)  )• 


reover,    a    comparison    of    (13,17)    ond    (14.3)    also   yields 
the   followin-  relations: 

TQ   -  ?8      m  O  (N0  (  ^  1 )    N0(  A2)    N0(  A8J     )  9 
?19   "  T19   -£?(N0(Ai)    N0U2)    N0(A2)    >• 

Finally,    from    (14.8),    the    seven  relations    above,    and   the 
equations    (13.20.1),    .    .    .,    (13.20.22),   we    conclude  that 

(14.9)  N0(\P(s,  r  )-4*(»t.T  ))-    0(N2s(Al)No(  A2)N2s(  A2)    )• 
Consequently    (14.7)    may  be   replaced  by 

(14.10)  N  (f  C-i  VfWpJ-jOi  &JP,  ^xf)) 

Hence  we   conclude,   from   (13.21)   and    (14,10), 
that 

(14.11)  /V  (a*?6J-  r  C-;  Vfl.'Vf  » 
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uf  course,  In  the  C3se  to  which  the  remork  following  (13.51) 
applies,  i.e.,  when 

and 

with  /^   independent  of    q* ,    (14.10)   may  "be  replaced  by 

(was)  N(T.(-iw,/*\t)-  r<s  Vf«/*M'0 

so  that   from   (13.22)    and    (14.12)   we  have 
(14.13)      N  (L^C-i  -  S  O  I  \<? , /*■  *t.'f  >) 

instead   of    (14.11).        L'oreover,    from   (13.19),    (14.2),    and 
(14.3) ,    it  follows   that      d  19  i(S;    r?  (s)  ;   A  1  Kp   ,  A  2  i^>  )    is 

bilinear  in  the  disturbances    (A^  q0,    ^  l  ^f"  }    and 

(^2  <lo>   ^  2  'y  )•        Consequently,   we   conclude   from    (14.4) 

that    ?  (s;  A  i  t     ,   &  2  ^'   )    aiso  has   this  property. 
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15 •        An   Improved    Approximation    for  the  Actual  effect  i 
of   a   jisturbance  on   the    Functions    q    (s). 

are   now   in   n  position   to   turn   our   attention 
to   the   problem  or   approximating     Ai    qMs)    by  an  expression 
of  the    fern 


lT?i(s)+|-ri(B;A1(pfA1y>) 


2      3 

and  of  obtaining  an   esti^nate   of    the   error  of  such  an  ap- 
proximation . 

Let  us   suppose    that  v/e   are   £iven   a   disturbance 
(^  coi  ^  '"f  ) ,    and  let   us  continue   to   suppose   that    (q0,  f   ) 
and    (q0  +    £  q0,    (_p  +   a  ^   )    aTe   in     4>ss.        Then,    as   a  con- 
venience   in  what   follows,    and   better  to  emphasize    the   de- 
pendence  upon   qo,   <<f       of    the    solution  of 

(15,1)        / 

q      (si)    s   qQ    , 

let  us  denote  that  solution  by  the  symbol 

q(s;  qojf  )» 

Because     A  r   =  0  we  do  not include  F  in  the  symbol. 
Now  let  •»>     be  any  positive  integer.   It  is  easily  veri- 
fied that,  for  j  -  1,  2,  .  .  .  ,   \>   -  1,  the  points 
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(flo  ♦  -i-    A   flo,         «p      ♦  -i-   Acp   ) 


are    in     4>ss.  rom   (8.8)   artf      11,3)    it  follows  th?t,      or 

each  i   ■  1,    .    .    .    ,    n,    and    J  ■  1      .        .    ,    v> 


But   from    (14.13)   we    conclude   that 

J    L 


■15  3    1?^.^.+  ii*^,   f+  Vr;  ^n.>t"   f^ 


ore  .  ver 


(15.4)        NQ(~-  Aqo,~~    Af)    =  ~"     N0Uq0.  *f)      , 
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and 


Since      j    -   1    =     ^   ,    we  also   have 


(15.6)         N2S(J~     ^qof^^f       -N2s(Aq0,^f)     . 

Hence,    by   comb   niif    (15   c)   with    (15.3),    utilizing    (15.4), 
(1    ,5),    and    (15.6),    and  b     recalling  the   linearity  of  the 
VO  s      a   ^      '  e   bilinearity  of  the      "jx  i      w©  obtain 


Hepla    inr?      ti      in    (15.7)    by  1,    2,    .    .    .    ,   v>   ,    consecutively, 
and  adding  the   results   yields 
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i 


A_    ,f    ^/V,   f  A  ^   Af      /V        CA%i   A  rfj^ 


upon   taking  liiuits    as       v>  -*-oo   ^    Ve   find  tha 

2. 


+  c  r^fH./^c^/^^^]^ . 


Finally     when  we  write 

qi(s)    2   qi(s;    q0,  Cf  )    , 

(s)    «   qHs;    q0    +Aqol   if  +  *  ^  )    -    q^s,    qQ,  <f    )  , 
*l    its)    =^(3;    q0,    (f   ;    Aq0,A<f) 
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and 

It    is   re8di.lv   seen   from    (15.9)    that   there   exists  a    constant 
K3  such  that 

(15.10)  a/(&i<-> -•?<•■>  -%,Z«>) 


i   K     ^^..^Ks^V*^ 


Thus  we  have  not  only  shown  that,  as 
^2s  (  A  Q-o »  A  ^  )  tends  to  zero,  the  quotient 


/Vf*<fCO  -  ^O  -  A-  yro) 


l>*^nOJ  Af ;  3 


also  approaches  zero;  (15.10)  is  a  somewhat  stronger  con- 
clusion than  the  statement  that  the  function 

J    (.)     =         J     ^(S)      I     81    «     S     *"  S2    ^ 

defined  by    (14.4)    with   A^   qQ=    A  g  qQ  s   A  q0   , 

A  2tf  *  A  xCf>  5A^,    is  the   second  differential  of   q(.) 

with  respect    to  the  norm     Ngs(A  q0,  A<f) . 
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16,     The  second  variations;   The   .jecond  Variation  aquations* 

We  may  now  obrerve    that    (14.4)    can  be   expressed 

Ln  terras  of  the    given   functions     /*,   tf°* ,  A  if  *"  ,   and  initial 

i  i 

<lo   »    A  q.3    •        For  if  we    compute     19  **j      from    (13.3), 

put     A  q0  in  place  of    A  ^q0     and   A  2q0   ,   put    A<f  in  piece 

Agf     In    equations    (3.9),    (14.1), and    (14.3), 


and  if  we  note   that      G 


•  0  for  »  *   i,   G" 


r 


a  1,  we 


fine  that  the  system  (14,4)  becomes 


(16.1) 


'lJ^Bfq(s)))%i(Bfq{8))JJ 


GS 


*  P^rfi(f(s»^8)))H)^e»<l(s))f^(s,q(s))^(s)^K(s) 

*  2Fr«r(|(^8»^8)))^y-<s»ti(8))A<f<(s,q(s))l?J(s) 
+  F?1,r^(«f(Bfq(«)))Af,(s,q(t))Af(f(8,q(8)) 

*  2F^(«f(efq(i)))A<f£  (s,q(s))>?j(s) 

+  ^(^(8,^8)))^^   (8,q(«))t?J'(s)9K(s) 


3fc(^)  - 


[n  fact,  (16.1)  is  the  system  which  one  obtains 
differentiating  (8.7)  with  respect  to  f  and  then  set- 
r  ■  0  in  the  result.    Thus  the  functions   Ti(s) 
(s,0)  are  identical.    .,e   shall  call  them  the 
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"second  variations"  of  the  funotione     q   (a).        (see    [fij, 
p.   165) •       The   equations    (16,1)  which  determine  the   second 
variations  we  shall   call  the   "second  variation  equations fT. 

If  we  oompare  (14.4)  with  (8,9),  we  see  that 
Theorem  (9.2)  applied  to  the  function  y  (.)  defined  by 
(14*4)   yields   the   result   that 

N    (  I  (.))   -  K  No(0,   d  tf  ) 

From   (14*3)   we  have,  when  A  ^  *    A  2, 

No(0,    dtf  )    -  6  lN0(Aq0,A(f>)   Ni8(Aq0»  ^  <p  )    )  . 

Hence  there  exists  a  constant  K*  such  that 

(16.2)       N(?(.))±  K*  N0(Aq0,Af)   Ni8(Aqo,A<f  ). 
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Chapter  V 

S3Q0NP  DOTBBfflTTXAl  mfBCTS  Of 
TRAJECTORY  SLBQCWTS 


17.     Effects  on  Trajectory  Elements  Compered  With  Effeote 
on  coordinates  of  the  Ma as- Center. 

Let  x(t),  y(t)»  z(t)  ha  the  funotlona  whloh,  to- 
gether with  their  derivativee  x(t)  •  vx(t),  y(t)  ■  Ty(t), 
s(t)  -  vB(t),  constitute  the  aolutlon  of  the  system  (6.3) 
aubjeot  to  the  normal  initial  conditions,  whloh  inoluda 
z(0)  ■  0  »  y(0)  -  z(0)»       That  la,  tha  ourra 


(17.1> 


la  tha  "normal  trajectory**   If  at  eome  time  t  •  T  >  0  wa 
have  y(T)  -  0,  then  wa  call  T  tha  time-  of  -flight*  and 
X-  x(T)  the  "range".   If  Ax(t),  A  y(t),  A  a(t)  denote 
the  aotual  effeota  on  the  coordinates  x,  y,  z,  respectively, 
due  to  a  given  disturbance,  then  the  ourve 

x  -  x(t)  ♦  Ax(t) 

(17,2)    \       y  -  y(t),  ♦  A  y(t) 

z  -  z(t)  ♦  A  z(t) 
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is   the    "disturbed  trajectory"*        If  T'    is  a  positive  value 
of     t     such  that     y(T»)   4     Ay(T')   -  0,      then 
X1   -  x(T')   ♦  Ax(T')    is  the  range   on  the   disturbed  trajec- 
tory and 

(17.3)  Ai-r  -  X 

-  x{V)   ♦     £  x(T»)    -  x(T) 

is  the  "range-effect"  due  to  the  given  disturbance. 

We  note  that  the  range-effect  A  X  is  defined  as 
the  difference  between  the  x-coordinates  of  a  point  R  on 
the  disturbed  trajectory  (17.2)  and  apoint  P  on  the  normal 
trajectory  (17.1),  the  points  P  and  R  having  been  "matched" 
(made  to  correspond)  by  the  requirement  that  they  have 
equal  altitudes  (y  ■  0  at  both  P  and  R) •   Since  we  can 
compute  A  X  by  (17.3)  only  if  we  know  both   Ax(t)  and  T', 
which  knowledge  can  be  acquired  only  by  computing  the  dis- 
turbed trajeotory  (17.2),  and  since  it  is  this  additional 
computation  which  the  ballistician  wishes  to  avoid,  we  need 
a  method  of  approximating  AX  with  reasonable  accuracy. 
Our  theory  of  differential  effects  should  furnish  such  a 
method. 

But  most  trajectories  are  computed  with  time  t  as 
independent  variable,  so  when  we  attempt  to  apply  our  theory 
of  differential  effects  we  find  that,  while  we  can  with  a 
fair  amount  of  ease  obtain  a  first-order  approximation  fon® 
of  the  functions  »pi(t)  defined  by  (8.9)  for  s  «  t  ]  to 
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A  x(t),  we  obtain  no  immediate  information  about  the  value 
of  T*.   The  fundamental  difficulty  is  that  the  variation 
equation?  enable  us  to  find  approximations  to  the  actual 
differences  in  the  coordinates  of  two  points,  P  on  the 
normal  trajectory  and  4  on  the  disturbed  trajectory,  which 
correspond  by  virtue  of  having  equal  values  of  time  t; 
whereas  we  require  an  estimate  of  the  difference  between 
the  x-coordinates  of  p  end  a  point  R  (generally  different 
from  Q)  which  match  because  they  have  equal  values  of  an- 
other variable,  the  altitude  y.   Similar  difficulties  are 
experienced  in  attempting  to  approximate  the  effects  on 
time-of-f li^ht,  on  deflection  (departure  from  the  plane  of 
fire),  and  on  other  elements  of  the  trajectory. 

If  we  are  to  be  able  to  apply  our  theory  of  dif- 
ferentials, then  it  is  clear  that  we  must  be  able  to  ex- 
press differential  effects  computed  at  equal  value8  of  one 
variable  in  terms  of  differential  effects  computed  at 
equal  values  of  some  other  variable. 

18 •   Relation  Between  First  Differential  effects  on  Trajec- 
tory Elements  and  First  Differential  Effects  on 
Coordinates  of  Mass-Center. 

The  objective  outlined  in  the  last  paragraph  of 
the  preceding  section  was  attained  for  first  differential 
effects  by  Bliss  in  reference  [2] ,  through  the  use  of  cer- 
tain theorems  on  implicit  functions  contained  in  his  earlier 
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paper    [l].        In  a   simpler,  though  perhaps    less  general 
manner,   1  cShane,   Kelley,   and  Reno  develop  a  similar  rela- 
tion  in  their  book    [e] .        It   la   this  latter  relation  whieh 
we  shall    state   in  the  present  section,   since  we   shall   find 
it  useful   in  deducing  an  analogous   formula  for  second  dif- 
ferential effects   in     #  19. 

Let  P  be  any  point  of  a  normal  trajeotory.       Let 
A  and  B  be   any  two  variables  whioh   can  be   used  as  indepen- 
dent variables  on  an  arc  of   trajectory  inoluding  P.       Let 
Z  •   (AqofAf7  )   be  any  disturbance.       Let  C  be  either   (i) 
the  value,   at  a   specific  point  of  a  trajectory,  of  one  of 
the   fundamental   quantities      q*(s)      defined  by    (7.1),    (7.2) 
or    (ii)    a   twice   continuously  differentieble  function  of 
the   q*.       Let  Aq,  Be,  C0     be  the  values  of  JL,    B,    C,    respec- 
tively,  at   the  point  P  of  the  normal  trajectory.       We 
shall  aasume   that   the  norm     N2a  ■  Ngfl(A  qQ,  A  <p   )    of  *ne 
disturbance     A     is   sufficiently  small  to  inaure  the   exist- 
ence of  points  Q  and  R  on  the   dlaturbed  trajectory  auch 
that    (i)    at  ^  the  variable  A  has  the   value  Ao,   and    (11)   at 
R  the  variable  B  has  the  value  B0.       Let     AB(A  |  A  -  A0) 
denote   the  difference  between  the  value  of  B  at  Q  and  the 
value   of  B  at  P;   similarly,    let      A  c(  X  |   A  -  Ao)    he  the 
difference  between  the  value  of  C  at   Q  and  the  - value  of  C 
at  P.        That  is,    at  £  we  have 
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(18a)   A  -  Ao,   B  -  Bo  ♦  AB(S|A-io)i 

e«0o*A0(£|A-Ao)* 

Likewise,  we  denote  the  values  of  A,  B,  end  0  at  g  by 

(18.8)   A  -  Ao  ♦  AA(SlB-Bo).  B  -  Bo* 

0  -  Co  ♦   ACi  &  |  B  -  Bo)» 

• 

reepeetively.   Finally,  if  we  let  d0(  A  |  A  *  Ao)  denote 
the  f irst  differential  effect  of  £  on  0  at  A  -  Ao,  end 
if  we  assign  similar  meanings  to  the  symbols  d0(£|  B-Bo) 
and  dC(  A  |  B  -  Bo)  end  dB(  &   |  A  -  Ao),  then  the  rela- 
tion eatabliehed  by  MoShane,.  Kelley,  and  Bono  (  [d],  chap- 
ter VII,  f  1)  la 

(18.3)  dO(Z|  B  -  Bo)-  A0{  A |  A  -  Ao)-  0*  (B0)oB(2  U  -  Ao), 

in  which  the  symbol  0v(Bo)  meene  the  value  at  B  •  Bo  of  the 
derivative  C (B)   of  c  with  r aspect  to  B  on  the  normal  tra- 
jectory* * 

Tor  example,   to  obtain  the  first  differential 
effect  on  range  attributable  to  a  disturbance   A    (the  oaae 
discussed  in    §  17),  we  apply   (18.3)   with  C  ■  z,  B  ■  y, 
Bo  -  0,  A  -  t,  Ao  -  T,   and  x»  (y)  -  i(t)/y(t).       If  we  de- 
note the  angle-of-fall  by  *>    ,  the  application  of   (18*3) 
with  the  substitutions  just  indicated  yields 
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(18.4)     6x(  2  |  7  -  0)   •  <bc(  2  |   t    <*  T)   ♦  cot  £  dy(  K  \  t  «  1) 

*  dx(2  |  t  •  T)  -  $Hi    6y(A  |t  •  t) 

7(T) 


as  the  first  differential  sffsot  on  rango  dus  to  ths  dis- 
turbsnos  S  •     If  i  •  t  and  ths  variables  q1,    ,   .    •   ,  <$  in 
(7*1)  srs  x,  y,  g,  vx,  vy>  v1#  respectively,  then  ths  dif- 
fsrsntlsls     dx(2|  t  •  T)  and  dy(JT|  t  »  T)  in  (18.4)  &rs 
ths  values  st  T  of  ths  functions  $  *>(t)   and  >?~(t),     ra- 
spsotivsly,  dsflnsd  by  (6*9)  with  &  •  t. 

19.     A  Rslstlon  Bstwssn  Sscond  Differential  r.ffaots 

Coaputed  With  Two  Different  Indsay- '■  dent  vgrlablas* 

If  ws  ettsmpt  to  iaproTS  ths  sstimats  furnished 
by  (18*4)  for  ths  range  sffsot    A  1  by  v,tili7,iuig  ths  second 
variations     5  *•  dsfined  by  (16*1)  whsn  &  •  t,   ws  fsos  ths 
same  difficulty  whi oh  confrontsd  ua  in    J  17:     ws  rsquirs 
an   (improved)  sstiasts  of  ths  difference  bstwsen  ths  z» 
coordinate t>  of  ths  points  R  and  P— -matched  bsi.suss  thsy 
hsrs  equal  altitudes—whereas  ths  variation  aquations   (both 
first  and  sscond)   furnish  us  only  ths  tenia  for  estimating 
the  differences  between  the  coordinates    (end  their  time- 
derlvstives)  of  the  points  Q  and  P,  which  correspond  by 
virtus  of  having  equal  vsluss  of  time.       That  is  to* say,   if 
ws  wish  to  smploy  sscond  differentials  to  improve  our  ap- 
proximation to  ths  range-ef fact    A  I,  we   shall  have  need 
for  a  formula— analogous  to   (18.4)—  which  enable  a  ua  to  da- 


66 


duoe  second  differential  ef facte  computed  at  equal  values 
of  one  variable  from  eeooad  and/o?  first  differential  ef- 
feote  computed  at  equal  values  of  some  other  variable.       It 
ia  our  purpose  in  thia  eeotion  to  develop  auoh  a  relation* 

We  shall  uae  the  notation  introduoed  in    f  18, 
and.  in  addition,  we  ahall  denote  by  dsc( £ |  A  -  Aq)  the 
second  differential  effeot  on  0  it  A  •  io  attributable  to 
the  diaturbanee  E  ;  i*e.,   the  second  differential  defined 
in  (lfi.8)  for  the  case  in  which    &iH>    lA.lp   i    S, 
Be  aaorlbe  similar  msaninga  to  d^cC&I  B  •  B^)  and  to 
d*B(&  |  A  -  Ao)»  and  mi  let  G*(B)  denote  the  seoond  deriva- 
tive of  C  with  respect  to  B  on  the  normal  trajectory* 
Then  the  relation  which  we  ahall  ®slabliah  ia  that 

(19.1)      dSC(A  |  B  -  B0)   -  A8C<£  |  A  -  Ao) 

-  C*(B0)   d«B(£|   A  «  JLo) 

♦  0»flM      [dB(I|  A  -  Ao)]8 

-  SdC'C*  I  A  -  Ao)*dB(Z|  A  -  Ac). 

Toward  thia  end  we  note  first  ,that,  if  B  ia  the 
Independent  variable,  then  from  point  %  to  point  R  on  the 
disturbed  trajectory  the   change  in  B  is     -^B(&|  A  •  Aq) 
while  that  in  C  is     A0(&|  B  «  B0)   -  A0{S  |  A  -  Ao)  • 
Hence  by  Taylor's  Theorem  there  exists  a  number  Bi,  be- 
tween    B0  and  B0  ♦  AB(£|  A  -  Aq)  ,    such  that 
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(19.2)  AC(  A  1    3  -  Do)    -    AC(  A  |   A  -  A0) 

-   -   A  B(ff  |    A   -  A0)    £•  (B0   ♦   ABlA  |    A  -   A0)    ) 

»   |   C-    A  BIX  |    A  -   AQ)      ]2  ^»(B1)     § 

Wherein   the   symbol    C '  (Bo  ♦AB(5|    A  •   Ao))   means   the    value 
of  the   derivative   of   C  with  respect   to  B  at   the   point   on 
the   disturbed  trajectory  at  which  B  ■   Bo  ♦  AB(  £  |   A  ■  Ao), 
with  a   similar   interpretation   for  C""(Bi).        From   (10.1)   we 
see   that      AC     will  not   exceed   a  finite   multiple  of 
No  "  N0  (A).        This   applies   equally  well    to  A  B,    AC*,    end 
A  C*.        In  particular,    the   values  of  CO)   at  the  points 
of   normal  and  disturbed   trajectories  with  B  •  Bl  differ  by 
at  most  a  bounded  multiple  of  N0.       koreover,  Bl  differs 
frors  Bo  by  at  most   a   bounded  multiple   of  No»    so  the   values 
of  C*(B)   at  the   points  of  the   normal   trajectory  corres- 
ponding to  3  -  Bi   and  to  3  »  Bo  also  differ  by  at  most  a 
bounded  multiple    of  N0.      Hence 

(19.3)  7*(3i)    -   C*(B0)      ♦  €  i  Nol 

in  which  €  ±   is  bounded.   3y  hypothesis, C (B)  is  a  twice 
continuously  dif ferentiable  function  K(q)  of  the  q*,  so 
the  equation 

C»(B)  -  K  j.(q(E))  q'^B) 
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expresses  Cf(B)  as  a  continuously  differentiable  function 
of  the  q1.   If  B  can  be  used  as  independent  variable  along 
the  whole  trajectory  (and  not  Just  alone:  an  arc  which  in- 
cludes F)  ,  then  by  (11*3)  the  q*  are  differentiable  func- 
tional® of  the  disturbance  A   ,  and,  by  a  slight  adaptation 
of  a  proof  to  be  found  in  most  advanced  calculus  texts 
(e.g.,  [is] ,  p.  77;  or  [13] ,  p.  72),  it  follows  that  the 
differential  dc * (  2t  |  3  -  B0)  exists.   On  the  other  hand, 
if  B  cannot  be  employed  as  independent  variable  along  the 
whole  trajectory  (as  is  the  case,  for  example,  when  B  ■  y) , 
then  theorem  (11.3)  does  not  apply.   But  for  the  important 

case  B  ■  y,  I..cShane  has  shown  (  [Wcs] ,  theorem  (11.7)  ) 

tk<n 
that  even  the  q^-  are  differentiable  functionala  of  the 

disturbance.   His  method  of  proof  depends  upon  the  intro- 
duction, through  a  transformation  of  coordinates,  of  a  new 
variable   s  capable  of  use  as  independent  variable  along 
the  whole  trajectory  and  uniquely  determined  by  y  when 

rt.  r\u»fccr       L«  *r      Him*,     bit  m.r  Attr-Arl  lu       *AjCr» 

y  does  not  exceed  the  maximum  ordinate  of  the  normal  tra- 
jectory.   In  Chapter  VII  we  shall  describe  in  more  detail 
the  coordinate  system  used  by  McShane;  it  suffices  here  to 
say  that  the  method  is  quite  general  end  can  be  applied  to 
any  variable  B  which  can  be  used  as  independent  variable 
on  an  arc  of  trajectory  including  P.   Thus,  in  either 
case,  the  existence  of  dC'(&|B  ■  Bo)  is  guaranteed.   Of 
course,  a  similar  conclusion  holds  concerning 
dC*(A|  A  -  Ao).   By  (11.31),  or  by  (11.9)  of  [kes]   in 
case  B  »  y,  we  have 
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2 
(10.4)     C'(D0)  *  OMPo)  *  duf  (  a  I  3  -  Bo)  ♦  c  2  "Is  , 


in  which  Nis  -  Wls( a  )  and   £  g  is  bounded. 
Hut  fro-n  (15.10)  it  follows  that 

(19.5)   AC  (a  |  A  »  A0)  -  dC(S|A  -  A0)  ♦  (l'2)d2C(A|  A  »  A0) 

+   €  3  N2s3 

and 

.  )  AB(a|a  «  A0)  =  dB(A|;4  -  A0)  ♦  (l'8)d23U|  A  -  A0) 

♦   e  4  N883, 

which  £  3  and  £  4  are  bounded  and  *:g3  •  N2g  (£). 
>rder  to  simplify  the  not&tion  somewhat,  throughout  the 
remainder  of  the  present  argument  v:e  3hall  write  Ab,  Ac, 

1  etc.,  for  A3(A  I  A  •  A0),  AC(£  |  A  ■  Ao), 

;  A  I  A  ■  A0) ,  etc. 

Using  this  abbreviated  notation,  we  observe  th 

Taylor's  Theorem 

(19.7)   Cf(30  UB)  «  C»  (B0)  +  A3-C"(B0)  +  €5   N02 

in  which  €  5  is  bounded;  and  also  that,  by  (18.3)  we  have 

119.3)    dC»  (  £  I  B  -  E0)  -  dC»  -  C"(B0)  dB. 
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We  are  now  ready  to  combine  equations   (19.2) 
through    (19.8),   inclusive*       The  process  will  be  eesier 
to  follow  if  we  perform  it   in  sevoral  steps. 

First,   we  transpose  the  second  term  of  the  left 
member  of    (19.2)   to  obtain 

9.9      ACU|S  -  Bo)    -    AC   -AB'?1  (B#  ♦  AB}*   jIaBJ^CBi)  . 

We  then  substitute   (19.7)  and  (19.3)  in  (19.9). 
Thie  yields 

4C(£|8-Bo)-AC  -AB  CMB0)-    [j&B]2[c«  (Bo)-  |c«(B0)] 

-   *5N0*AB*  |61^02   [ab]2. 

But  since  neither  C*(Bo)   -  C*(Bo)   nor  AB     exoeeds  a 
altiple  of  N0,   we  conclude  from  (19.10)   that 

(19.11)    AC(&|E-B0)    "AC  -  AB-C1  (B0)-  |  kEj*-C»(Bo.H  *  6NQ3, 
in  which  €  §  is  bounded. 

Next,   we   substitute    (19.8)    in   (19.4)   and  put  the 
result  Id  place  of  <J*  (B0)   in  the  right  member  of   (19.11). 
V/e  thus  obtain 


(19.18)     AC(a|B-B0)!    AC  -ABtcMBo)    ♦   <*C*    -   C*(3j    dB  ] 


-  |  Tab]2  c"(Bo)  *€7  Ni83  . 
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with  t  7   bounded,   since  No  ■  Nifl  and   A  B  -  0(No).        But 
from    (11.31)   we  have 

(19.13)  AB  -  dB  ♦  €8  Nls2 

in  which  £  9  is  bounded;  end,  since  from  (9.2)  it  follows 
that 

(19.14)  dg  -  ;-:0  hqi 

we  see   from   (19.13)   and    (19.14)    that 

(19.15)  [>b]2    -    (dB)2  *€g   Nls3 

with      €  g  bounded. 

Now   in   the    first   two   terras  of    the  right  member   of 
(19.12}   we   replace    AC,     AB  by  the    ri  ght  members  of    (19„5), 
(19.6),    respectively;   and   in   the   third   terra  of  the  right 
member   of    (19.12)    we   replace      [ab]2     by  the   right  member 
of    (19.15).        Wc  find   in  this  way  that 
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1    ,2.  ^     „      3 


(19. lo)      ac(a|b-b0)   e  dc  t  -d  c  ♦  e  3M?. 


2 


:g 


-  C'(B0)    dB   -   I  C'(B0)    d2P-    G» (B0)  ^4>T283 

-  dC'dr.    -    ^   dCM2B   -    dC     €4    r:2g3 
♦    Ctt(Bj     [dB]2*    i  C'w(Bo)    dBd^B 


♦    C^-^o)    ^  4   ^28^   ^ 


-  j^"^c>    M2"  ^6^is3  c-(ae)>  c7:.l33. 


But    fro::     (lo.£)    we    set:    t  .at      d23      does   not    exceed   a    bounded 
multiple    3f        N'ls    1    30,    rince   fro:n    (9.2)    v-e    have    do'    -   KqIIq, 
it    follov.s    fro:,    11  j.  10)    that 


(19.1' )        4  :(  a  |:.'--0)    -    0.    -    C  IB0)    dB 

♦  ^{d2C  -   CMP0)dan  -  2    du»dBtc"(b0)[d::]2} 

*  C10   -2s3» 

in  v/bich    £  lo    is   bounded. 

Finally,   since      do  -    C1  (Bo)    dB   is  by    (13.3)    the 
first    differential      dc(  a|    B   -  Bo)  »    it    follows    that   the 
expression    in    braces    {*••]    in   the   ri'ht   member   of    (19.17) 
i3      d2C(£|    '^   ■  Bo),   the    special    form  assumed   by  the   second 
differential — defined    in    (12.^) — for    the    case 
&1{V   =    &2Vm~K     .        -phis  establi-s-  es    (19.1). 
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As   an   illustration  of    the    applic    Lion   of    (Ij.1), 
let   up    find    the   second    differential    effect      d2x(  Z  \    y  *   0) 
on  range    due    to  the   disturbance  A  .        Just   05   in  deducing 
the    first    differential   ranpc    effect    (1  >. 4),   v/e    let   0    -   x, 
B   *   v,    Bo   =»   0,    A  ■   t ,   Ao   =    T.      Hence 

CME)    ■   x' (y)    =   vx(t)yvyit).  en,    with   the   help   of   the 

nornol    differential    equations  of  motion 

vx   =   -   .5vx 

Vy    =    -ivy   .-■"«.• 

we   find    that 

(19.15)         C»{d0)    -   X"(Q) 

vy'.f)  Vx(T)    -  Vy,(:i')vy(T) 

~Xvy(T)] 

_      g   Vx(T) 


Lvy(T)]3 


.  ore  over 


(19.19)       dC'Ul    A=A0)    -    [vydvx   -    vxdvy]    ,     [v^] 2 


in  which  vx,    Vy  are  evaluated  at      t  =  T  and 

dvr=   dvxU"|t   =   T)  ,    dvy  =   dvy(S  |t   =   T)  .  Substituting   the 

appropriate   values  for  C,   B,    .    .    *    ,    do'  in    (19.1) 

yields 
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119.80)      d2x{2|y  »   u)-   d^x(a|t   -  T)-   ^~~-  d2yU  U   -   T) 


dy  (  A  1 1  «   T) 


-    {  vy{T)dvK(  Z  |  W) 


Cvy(T)J2 

-  vx(T)dvy(  A|t-T)  J 

[vy(T)J3 

In   case   the    striking  velocity 

vO    "      (    f^x(T)]2  *    [vy(T)]2}    1/Z 

and    the    first   differential   effect  d  <3    on  tvie  arigle-of-f all 
&   have    eiready   been   calculated  for  the    trajectory  to  which 
one   wishes    to    apply  the    formula    (19.20),    tie    latter  equa- 
tion  can    be    written    in  a    somewhat  iiore  conveniently  applied 
form.         /or 

(19.21)  *3     *     TC     -  0(t) 

in  which     9{T)    is  the    velue   at   t   »   T  of   the   an^le   &{t)    of 
inclination   of   the    tangent  to    the    trajectory.        Conse- 
quently,  s  ince 

(19.22)  C»(3)    -   i»(y)    *   vx(t)/vy(t)    »   cot   e(t), 

we  have     vx  \/vy(T)  *  -  cot  ^>   ,  and  (19.18)  may  be 
written  ir  the  alternative  form 
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{iv. i  c"(n0)  =  x"(o) 

=»   -P5   cot.    £  /[vy(T)]2 

-    -7   cot   ■£  csc2    <u)  /v^  2. 

Moreover,    fro-     (19.2;)    we   have 
(19.84)      dC  (  a  |A   -   Ao)58    -    esc2   6(T)    d©  (  2  |t    »   T)  . 

But    the    differential    effect    dO    on    the    en^le-of- 
fall  *>   is   the   differential 

(1.9.8  d^>    =   -deiJ|y-0), 

not    -de(A|t  =  T) .    / ortunately,  the  relation  (18.3) 
enables  us  to  express  this  latter  differential  in  ter:.is  of 
d  <*)  .    :iov;ever,  to  use  (18.3)  we  first  need  to  .calculate 
the  value  of  d0/dy  at  y  =  0  on  the  normal  trajectory.    To 
this  end,  v/e  observe  that 

(19.86;     x"(y)  =  d(cot©)/dy 

=  -  esc2  9  de/dy, 

from  which  it  follows  that 

(19.87)    (de/dy) |y  ,  0   =  -  x"(0)  sin2  cD     . 

Now   setting  c   =  G,    E  =  y,    3o  «  0,  A  =  t,  Ao  -  T  in    (13 

and   usin;?    (19.25),     (19.27),   we   obtain 

(19.28)         -   dcD   =   de(A"|t=?)    +    x"(0)    sin2  cD  dy(  a     t-3 

It  follows  from  this  last  equation  that  (19.24)  may  be 
written  as 
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(19.29;   dw1  (  A  U  "  A)-  cs°C  ^   d^  +  x*(u)  dy("2  |t  -  T), 

e 

in  which,  we  recall  fro;::  (19.22),  C  is  /.'(y). 

jou.<=,e'.;uo!itX.Vi  whe  i  we  ret  0  ■  x,  ^  ■  y,  Ho  "  os 
A  »  t.  Ao  »  I",  and  utilize  (19.21),  (19.22),  (19.23;, 
(19.29,  L  i  i  [  1 9 , 1 ) ,  v.* e  o bts  in 

(19.30    \-:y.[A\v   -  0)-d^x(A"|t  *  T)  ♦  cot^d2y(A"|t  -  T) 

♦    (p'v^2)    cot  *>csc2<3  [dy(£  |t«T)]2 
-    2    crc?'£d<2dy    (  5  |t   -    T)  , 

in  which   the   differential      d  <^  13   that  defined   by    (19.25). 
Thus    (19.20 1    13   en   alternative    for...    of    (19,20), 

our   improved   estimate   of  tie   ran^e   effect      A  A 
caused    by  e    disturbance   A    is   now 

(19.31)       A,.    «    dx(Z|y   -    o)    +    (1/8)    d2x(£|y   «    0)  , 

in  which   the    first   term  of    the   ri.^ht  lueiuber    is    given  by 
(18.4)    and    the   second   ter.u  is   given  either  by    (19.20) 
or  by    (19.30;.        ..oreover,   we    have    the    assurance   from 
(19.17)    that    the    esrvor  of    tiie    approximation    (19,31)    is 
at  rent    a   bounded  multiple   of      [h'2s    (a)33# 
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Chapter   VI 

HSDUCTICH   Of,  Ttiiii  J3 ECO : PD  _V  A . { L V r  1 Oj£3  _TO 

DOUBLE   STJELTJZo    I  ii'uGic'iLo 


2d •      The    Second   Variation   Zquotion.'.   and    the    .associated 

ilomogeneous   jyst.em. 

If      q^   *   qMs)      is    the    standard  trajectory  de- 
fined  by  the   normal    equations    (13.1),    and  if    A(f  is   any 
set   of   disturbance    functions,    then    in   the  conventional 
ballistic   notation,    let 

(20.1)        ei(s)    =   F^    (if  (s,    q(s)     )    A«fMs,    q(s))  , 

so   that   the    first   variation  equatior.3    (8,9)    become,    upon 
our  writing    A    inrtead    of    A\, 


<|*_  «   F*       («f(s,q(s)))<f>?(s,q(s))  ?«1    +     Ci(s) 
ds  r*  4  J 

tl    (si)    z   A  q0    »    U   "   If  ...»n    ;  *  =   l,...,m). 


In  order  to  apply  to  the  second  variations  "J  (s)s  <3rr(s,o) 
as  much  as  possible  of  the  existing  theory  concerning  first 
variations,  let  us  begin  by  using  the  symbol  £^(3) 
note  the  sum  of  the  last  five  terms  on  the  right  member  of 
the  first  of  equation* (16 .1 ) .    That  is, 
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4  £yi  (  f  l  •  ,  n  (B  )  ))  f  *    n  ,  ;    r  ))  Ay^  S|    ;  (  B  ))  ^  (fl  ) 


♦  ?* 


r-r 


*    27^     (fin.qis)  ))  Af    -(p. lis))    17  J  ( s ; . 


"•en    (lu.l;    :ztiy   he  v.rltt»-n   98 


(20.4) 


^i    lai)    e  0. 


It    Is   at  once  apparent,    fro-r.  a    comparison  of 
(20.7.)    end     (20.4),    tnet    these    systems   have  the    same    corr^s- 
ponding   homogeneous   system 


(20.5) 


^  "*«  (s)*J- 


in   which 


Aijlsi    a  r'1  ^  (<f(s,q(s)))  Cf  *  (s,ql?j), 


hs   w«\l    r„    t.N»    pg^    Md  Joint    gynton 


81 


(  20  .  7  )         ^-A_ 


Aj1     ip)    X  ^    , 


.     the   matrix      ||    ;vJi    (gJ     j|       , .M  n  ,   the    trenBpo^  of    ^ 
matrix       |i    Aij     (S)     ||. 


■>  ^    v"  J»     v       x »    •    •     •    ,    n    ,    denote    n 
pearly,  independent    solution?   of    (20.5 J,    and    let 


rzr  i 


p    (si    denote    the    co-factor    of    $  1(8)    in   the    determinant 


20.8)         d(s)    = 


5  i     (a)        ?   p(s) 


•  a 


;  f  i»> 


i 

?!  «•> 


•       • 


I  i  i  <•> 


*>>  I 


D(S!)         exp    f    |An(s) 


•    +Ann(s)}   ds 


Then  D(s)  4   0,  and  it  Is  well  known  that,  fb  r  each  >>  , 
the  functions 


(P.O.  9) 


i  (S)  S 


(s) 


D(b) 


,  (1=1,  . . .  ,  n) 
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constitute    s    solution    of    th<;   ndjolnt   oyster;'    (CO. 7);    and, 
moreover,    that    the    functions    IrjMs)      and    T"     (s)    may   be 
expressed    In    tertas   of   their   respective    initial    values 
f?i^sl)    »      *5     (fll)    5  0    ,   and  the    functions      A^    (s)    by 
means   of   the    formulas 


<r*  s 

(20.10)   r?i(e)s  5*U){  ^(.njr/ui)*  J        ^(<r)er(r)d<r| 

and 

TrS 

(20.11)     ^(sJe  jMs)[  V^S].)  3*(«l)*    f         ^(tU(t)dt] 


r*sL 


=  !>/ 


*'(*)     £+(r)   dr. 


^•yi 


In  order  to  simplify  the  notation  somewhat,  we 
shall  define  functions  K  (s),  M.    (s)  by  writing 

(20.12)  K  (s)a  ?\^(si)hr(si)  ♦  f      >  *  (<r)  er(^-)  4<r, 

11     r         Jr.sL 

(r-l,  .  .  .  ,  n) 

and 

(20.13)  *.  (s)s     ,X  *  (t)  €*(v)    dr,  (p«l,...,n). 
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,Ve  ere  then  able  to  write  (20.10)  end  (20.11)  in  the  con- 
densed forms 

(20.14)  ?i(s)  -   5^1  (8j  K.p(s), 
and 

(20.15)  ji(8,  -   fj(a)  ,*,<.), 

respectively.   The  functions  ^  (3)   are  well  Iciown,  and 
have  been  employed  extensively  in  the  calculation  of 
(first)  differential  effects  of  abnormal  forces  and  non- 
standard initial  conditions.    It  is  our  purpose  to  ex- 
amine, in  the  pages  that  follow,  the  analogous   (but  con- 
siderably more  complicated)  functions  ju.    (s). 

21.   The  Functions  >^-u(s)  as  the  Sum  of  Functions  of 
Four  Distinct  Tyoe s . 

V/e  note  first  of  811  that  the  functions  y^  ^  Cs  j 
are  themselves  dependent  upon  the  functions  K   (s).    in 
fact,  when  (20.14)  is  substituted  in  (20.3)  and  the  result 
used  in  (20.13),  we  find  that  we  can  write 

(81.1) /*,<■)  .AlM   ♦A*I(S)  *^i.)   */*.?[,)    , 

in  which  the  terms  of  the  right  member  are  defined  as  fol- 
lows : 
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(21.8) /4.,1")  s    /  J^fF*      0? °*.  Cf>^  +  F*   <f".      h\    tV<U 


(  V,  5  -  1,    .. .,   n) 


x*5        x  3       u-         ** 


»      ^        t    fc     »    1^     i    ^r     being  evaluate 

*p~    and   ?P  „ 


for  the    aixur.vr.t,    t  ,    the   functions   Fp  _,     and   Fp  _     /»     for 
prcu':.ent   LP  {  t ,    q('t)),    and   the    functions    a?*  ,    if  , 

VV 


and     l^t   •  for  the    arguments  Y,    q(t) ; 


*».f 


with    A^     evaluated   for  the    arguments      t  ,    q(t),    and   with 
the   remaining  functions    in   the    integrand    of    (21,3)    evaluated 
as   described  above    for    (21.2); 

with    Af    evaluated  for  the  arguments      t  ,    q(t),   and   the 
other  functions  in  the  integrand  of    (21.4)    evaluated   in  the 
manner  already  described;   and,    finally, 


»L> 


in    v/iich    &  (T    i       la   ovnluntod    for    the    flrpinpnts      f  ,    q("r), 
and    the    remaining    funcMonq    a  ro    evaluated    as    orevlouslv 


descMb^d . 


Thus    our    study   of    the   function?    /a.^  (s)    falls 
Rurally    Into   four   pert?;    viz.,    an  examination   of   *>ach   of 
he   functions  /+ $  (a)  , /t1^  ia)  §/U.   XJX  (8  )  t/^  XV{  s  )  ,    aepe- 
rately.        This  wo   proceed    to  do. 


22.      The   runctions^A  J*    is). 

in    this    section,    we    shall    find    it    convenient    to 
use    the    following   abbreviated   rotation: 


the    functions    in    the    ri^ht   aeaber  beln$  evaluated   as 
described   above    for    (21.2).        7/e    can   thus  write    (21.2) 
what   more    ruccinctly  as 


in  which  each  of  the   functions  appearing  in  the   in- 
tegrand is  evaluated  for  the  argument      x    . 
now  use    (20.1?)    to   substitute   for     *^(t)   and    ^    («r) 
in   (22,2),    the   result  is 
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tr>   «-  J  \Js.  ir)f^(t)   |f(t)  fc^V  *,c*i>  <*■* 


t.f  <r,x 


of  the  throe  terms  in  the  right  member  of  (22.3), 

the  first  gives  ^&~  (a)  when  the  disturbance  consists  only 

altered  initial  conditions;   the  third  furnishes/**-  *(s) 

when  the  Initial  conditions  are  standard  but  are  accompanied 

or 

by_  changes  in  the  f unotlona  if    (a,  q),  8uch  as  would  be  the 

case  if  abnormal  forces  were  acting  upon  the  projectile; 
while  all  three  terms  must  be  evaluated  when  one  wishes  the 
combined  effects  both  of  non-stenderd  Initial  conditions 
and  of  abnormal  forces.   For  later  reference,  we  shall 
designate  the  three  terms  in  the  ri^ht  member  of  (22.3)  as 
/*■»* II  is),/*     ^   (fi)»  and  A?3  (s),  respectively.  Since 
it  is  the  third  of  these  terras  which  is  the  most  compli- 
cated, and  since  it  is  also  the  most  Important  in  applica~ 
tions,  we  shall  turn  our  ettention  first  to 
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(22. 4 1 l^'c* J 5  LIt )  tjr)  J;t )|  J X|Vr; C^r)  i rj  tf(f) e\ ? ;  J V J <* 


Now,  the  right  member  of  (22.4)  In  the  sum  of 
many  Integral*,  each  of  which  —  by  virtue  of  (20.1)  — 
has   the   form 

t*x  cr.ir  r^" T  (3 


In  order  to  exhibit  th*  f  unctions  ^Jf  ir)  ,  (£  (<r) , 
0  (^)  explicitly,  90  well  as  for  the  purpose  of  discussing 
the  integral  l(ai,s)  and  similar  *xpre8Fions  which  will 
occur  in  leter  sections,  we  shall  find  it  convenient  oc- 
casionally to  abandon,  temporarily,  the  summation  conven- 
tion concerning  repeated  indices.   For  the  aeke  of  clarity, 
we  shall  adopt  the  device  indicoted  in  the  following 

\ 
Remark  concerning  Notation: 

When  it  is  necessary  to  repeat  an  index 

upon  which  summation  is  not  Intended,  we 

shall  inclose  that  index  In  parentheses, 

When  an  index  is  repeated  but  not  inclosed 
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in   parentheses,    the   usual   summation 
convention  will   be   implied, 

Thus,    Tor  example ,    because   the   index   <x  range a  over  the 
feet  of    Integers   1,    .    .     .    ,        ,    the   symbol 

F^Cytir.i  qW  )  )  Af^  (r,  %(r)  ) 

will    .lean  the    suni 

s,i1(f(^q(o-)))A,.(p  1(r,q(f)W---*?i|ll  (flMlrlJlAf^lrjqlr)] 
as  heretofore;   whereas ,    on  the  other  hand,    the   symbol 


a 
r 


1rc«>iLf>{  r»Yr)))  Atf     (crr  o(<r)  * 


will  indicate  any  One  term  in  the  aforementioned  sum. 

It  is  our  purpose  to  show  that  each  term  of  the 
sun  in  the  right  member  of  (22.4)  can  be  expressed  as  a 
double  Stieltjes  Integral,  end  I(si»$)  is  a  typical  such 
term.   In  the  notation  described  on  the  preceding  page, 
it  follows  from  (22,4)  end  (£0.1)  that  the  functions 
^(t),  $(*-)»&(f)  in  (22,5)  are 
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(J) 


M 


*w  ■^»)wwij;;  <*>  s ,  w 


tes.6)     \  $   (<r)   z  ^{*?    t<r)      F(JJ}  (cf(<rt    q{<r)    )    ) 
(©(?)-    *(«,)     (f)      l^»(f  (f  »    q(f)     )    ) 


13 


Cur  first.  step  to\    id  express' i.'"*   y^- *>       (s)     as 
t    double   otieltjoe   integral  will   be   to   show  that    (r.2,5) 
cen  be  written  as 

188.7)   lUx,.)-/     J     $^if,rJ$Cr>A?^«OWr>*f*&»<Hf 


in  which 


t-.  * 


rs  > 

r  • 


/    -^Ct)  4  T  ^t  t   ?  <  ^"J 


That    is,  v?e   shall  snow   thet 


c*  s  p«r 


(2£.9)    I(slfe) 


J 


k  K 


c  r*:S  ^  flu  ^ 


•w       * 


f=*    r=r 


<«> 


<w 


00 

That  this    is  indeed  true  may  perhap*  mo«t  readily  be   seen 
by  reference  to  Figure  1  on  page     88,     In  feet,   is  follows 
at   once  from    (22.5)    tho{; 

(22.10)    I{sl9z)' J ff  W^^^^^^^t^^^^ 

and   it   is   clear  from    (22.10)    that  the   three-dimensional 
region  of  integration   is   the   "inverted"  pyremid  0-ABCD 
sketched  in  figure   1,    the   bounding  planes  thereof  being 

OCD:      (2  »  si,  OAB:     9  -  f  , 

oad  :     <r  ■  s1 1  ob«5  ;    r ■■■.-  *  , 

and 

ABC©:     r  »  s^. 

Now,    a  glance   at  the  sketch  reveals   that  this 
region  is  the   sum  of  the   two  triangular  pyramids  0-A3D 
and  0-BCD,   which  have  face   OBD  (    ^  ■  4"  )   in   common.        But 
the   integrals  of  the   function 

•$-(*)$(*-)*¥  (<r,    q(r)    )    G  (f)A(p        (f,    q(f)) 

over  the    regions  0-BCD  and  0-ABD  are,  respectively,   the 
first  and    second   integrals  occurring  in  the  right  member 
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41  f 


Mane   !  CAB:  f  »'■  J?  j 
Plane      OBCi  t  »  <r- 
Plane      GBP:  ^  »  <?~ 


Figure  1 . 
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of  122*9) *   Consequently  (£2,5)  nr^,    (22.7)  or»  equivalent, 
which  is  what  wo  wished  to  show. 

in   the  second  ploce,  we  assert  that  (2?.. 7)  con 
be  written  as  the  double  jtleltje*  Integral 


22.11)    I(sx,»)«    JJif  (•"•qCrJUf      (f.q(f))6e(si,s;f,^), 


(2 


in  which  the  integrator  function  is 

(28.12)  Rtsi.fjf.r)*   /   /   ^r  (s;u,v)${v)@(u)  du  dv 

end   the  interval  of  integration     J       is  the   square 

J:     81  •   f    •         ■•         *l  *    r  "     *   • 

To  establish  this  assertion,  let  the  square  J  be  sub- 
divided into  hf  rectangular  sub-intervals  Ijj  by  the 
lines  f -fi  ,  <7"-*"i,  (i  -  1,  .  .  .  ,  h;  J  »  1,  .  .  .,  f). 
Then,  employing  the  customary  notation  of  Stieltjes  inte- 
gral*, and  us  in*?  the  abbreviation 

F(s;u,v)£  y   U;u,v)  <f  (v)  <9  (u)f 


we  have 
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<^.*> 


(W 


C  ■*  s  ) 


<stst)  f 


Figure  2 


?. 


mAgIl3*j    j      F(ss>u,v)du  dv  ~j   f 


F(s ,u,v)du  dv 


Si      *L 


5-i  „?, 


4-i 


*±    Ji 


cr 


■t-i 


1       F(s;u,v)dudv»l      j       F(s;u,v)    du  dv 
4.        ^  «i    *4± 


*>fc 


O-i     fM 


*tp-(s;u,v)  $(v),  ©  (u)   du   dv, 
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But,  by  the  mean  value  theorem  of  the  integral  calculus, 
there  exist  numbers  ^i,  I?  j  with  ^i-l  ■  ^1  ■  Pi, 
^i-1  -  *?J   -  *j   such  that 

(22.14)  Ag  uj  -f  (85ri,l?j)$(»7j)<9(lt)Lrr^-l][Plfi.l]  .• 

Consequently  from   (22*14),  the  definition  of   a  double 
Stielt.les  integral,   and    (22,7),  we  have 


(22.15)         SS       ACf  (ffjqlr))    A  <f  f(W(f  ,q(f  )  )dg(sx,0  Jf  ,*") 


T 


■"*     ^     ;S4fp(»?j,q(9j))Z«p     (ri,q(5i))Ag  Iij 


h  -><* 

f -^oo 


/^J^f.,^ 


-   I    (si,    s) , 

which  completes  tte   proof  of   (22.11).        Thus  we  have 
established  the 


9? 


(22,16)   Theorem:   Given  y>  ,  for  each  ordered  pair  of 
indices  [<X,(3)    on  disturbance  functions,  and 
for  each  set  (j,  k,  V,  <$"  ,  i,  m)  of  the  other 
§.   indices  in  (22.4) ,  using  the  notation  of 
( 22 . 6 ) ,  there  is  e_  corresponding  term 

I  oc  /3  ik  v  £  i  m  ( si » s  )   in  the  right  member 
of  (22.4)  which  is  expressible  as  a   double 
Stieltjes  integral  of  the  form  (22.11). 
Furthermore ,  all  these  integrals  [   corres- 
ponding to  different  sets  (j,  k,  n/  9  S  ,    i,  m) 
of  the  latter  indices  ]  have  the  same  inte- 
grand, viz., 

A<f°V,  qM)     4  Cf>  ^   (?,  q  (p))  , 

but  each  has  its  own  integrator  function 
depending  upon  the  set  (j,  k,  */,  5",  i,  m) 

and  given  by  (22.12),  (22.8),  and  (22.6). 

By  a  well  known  theorem  on  Stieltjes  integrals, 
however,  the  sum  of  the  integrals  just  referred  to  is 
equal  to  the  integral  of  the  function 

A  ^P  (<r,  q(<r))  A  cp  {ft    q(<^>))  with  respect  to  the  sum  of 
the  separate  integrator  functions.   That  is,  returning 
now  to  the  summation  convention  regarding  the  repeated  indices 
j,  k,  V,  <*  ,  i,  m,tf  ,y,   if  we  define  %*}    (si,  s ;  p  ,  <r  ) 
by  the  equation 
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VsT       u- 


(22,17)    5  yslfs;f,<r}»     i  /    ^js;  u,v)  $Jv)  ©^  (u)   «u^v, 

wherein    fcf.    t32,l<n),    (22. 8j,     (22. o),    arid    (25.1)    ] 

Z=  wax  Cu^v;  ' 

[  with  the   functions   in  the   integrand   evaluated   as   in 


(22.19,      $    (v)    t     J\V(v)    F^tC^y,    q(v)    )    j 


oC  t  T^ 


and 


(22.20;     g>-       (U)     5      J\        (U)     ^-(^(u,     q(u)      )      ) 

r  m  ^  r 


then,    for  each  ordered  pair  of  indices:   <*  f  ft     On  distur'o- 
ance   functions  the   sum   of  the   corresponding  terras  in  the 
ri  p-ht    member  of    (22„4)  ,    i^ 

(22.21)   J  J  4f      (<r,q(<r)Af    (f  ,q(f)  )^tfi         (si;?;?,^). 

Finally,    since   the   ri<?ht  member   of    (22 e4)    cc»P(»J.sts  of   the 
sum  of   all   such  integrals  as    (22.21)   for  all  possible   or- 
dered peirs   °c  ,  (1   ,    [«c,(3t  1,    .    .    .    ,   m  ]    ,   Jt  follows 
that,   extending  the   summation   convention  to  the   indices 
**>  t  P   »   the   equation    (22.4)    may  be  written  as 
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J-  ' 

in  which  the  interval  of  integration  T   i  s  the  square 

T  :  (?i .  =  «^  =  s;    si  =  <r   *  s) 


it   bein.s  understood   that,    the   right   •    ml  sr  of    (22.22) 
the   gun  of     i   c      double    Stieltjes    ini      ;rcls.         e    for  e 
ordered  pair    (<*  ,  /3  )    of   the    indices    «<    and     'c     „ 

Logically,    perhaps,   we  should   at    this   tine 
cuss    the    ots:r  two   tern?    in   the    ri  rTht  )er   o^    (22,3). 

But,    because   of   their   form,   these   terras   are  more   easily 

TI  IIT 

handled   after  the    functions  A^(s)i>^    ^(s)    have    beer 

treated.        For  the  methods    by  which  these   latter  func 
are   reduced   to  do   ble  Jtieltjes   integrals  may  be   app- 
ro ->dily  to   the  remaining  terms    in    (22.3).      This  will 
do;:?   in    f  26. 


23,      The  Functions   A*-  ^   (s)  . 


Upon  using    (20.18)    to  substitute   for    Kit) 
(21.3),    end   then  utilizing    (20.1),   we  find  that 


(SSoD/JJfs)    5   2    M^M    *    2   al  "*    (s) 


m  which 


*9 


/OS 
fc^U],)    L      (*)4f  (*,    q(r))    dr, 


and 


(23.3)/^  ^     (^3/        /         Lv    (f)M^  (^)Af  (*,q(^)A«f  (<r,q(<r)   dtfdt, 


wherein,  in  turn, 


(23.4)  i/  WffX'trjF^.   (<f>(r.q(r)))f  *,  (*,q(-e))  ?•*(*)  1 

and  [  Cf .  (22.19),  (22.30)  ] 


(23.5)   M^(cr)  -  AV(<r)  F1  '  (  «f  (<r,  q(r-)  )  ). 

Ill 
Now  the   integrals  /-c  ^    (s)      are  formally  like 

III 
those  designated  >u-    .    (s)      in    (21.4),    in  that  only  a   single 

quadrature   is   involved;    consequently  we   shall  defer  treat- 
ing them  until  after  the   functions  yu.   ^    (s)      have   been 
discussed. 


nre  turn  our  attention,    then,   to    (23.3)   and  note 
that,   as   in  the   case  of    (22.4),    its  right  member  is   the 
sum  of    a  multitude   of  integrals.        3ach  such  term  in   this 
sum,    nowever,    has   the  form 


f      L(T);£(«-)Af"(rfqhrj)A<f  (<r,q(<r;  )d<rdr( 
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In  fact,    ;Tiven  >>   ,   for  each  ordered  pair    {<*  ,  A  )    Oi    the 
indices  on  disturbance  functions,    and  for  every   set    (p, 
£»    ji  V,    ij   of  values  of   the  olher  5    indices,    the   corres- 
ponding functions  L(*r)    and  k(<r)    are  riven  by    (23.4)    and 
\23.5)4respectively.        Our  next  task   is    to   show  that  the 
typical  integral    ZT (si,s)    defined   in    (23.<>)    can  be  written 
as  the   double   ..tieltjes   integral 


(33.7}  Jf  A  <f*(r,  qU))A<f   (<r,c(<r))   dg   (ei;r,t), 

in  which  the   integrator  function     g(si;<7",^)   is 
(23,3)      g(ai;^,r)    s|        J         L(v)  M{u)   KE(u,v)    du   dv, 
the   function       K^((T,r),   defined  by 


(23.9)        Kg    (  <T,  t  ) 


0,     s^  *"  T   <  r   ; 


*e 


lf      <r    a    t    -      s    , 

being  the  characteristic  function  of  the  set 

2  *  {  (  *•  t  r  )  |   si^*o-«r=   s} 

which  is   the   closed  triangle   GAB  in  the  figure   on  page 
97. 

To  this  end,   let    €    be   any  positive   number.     Let 
No     denote   the  norm  NQ(  A  qol  A  if    )    and   let    ty      be  an 
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A(*,s) 


(W 


Figure   3. 

upper  bound  for  the   absolute  values   of  the  product 
L(v)L(u)Ke(u,v)   on  the   square   3:        By  first   choosing  an 
integer     h     such   that 

h  >  2tyN02(s  -   si)2   /€    , 

we  next  partition  S  into  h2  squares  1^ j  by  means  of 
the  lines 


wherein 


s 


10X 
/'  a   _    «    \ 


j     "    81       ♦     j 


/    r    -    s] 


for  i    a   1  i    .    .    •    »    h,    and    ,1   =   1,    .    .    .    ,    h.        Then  thw   sum 
of   the   areas  of   the      h      squares   Ijj     (i  ■  1,    •    .    ,    ,    h)    is 
evidently 


,  o 


(g  -  sir-    <        * 

Nov  the    double   Stieltjes   integral    (2.5„7)    is  the 
licit,    as  h  ->•  oo    ,    of    the  ti.ouble   sum 

h  h  ...  (3 

(23.10)     ^      -Y    A«f   '•>?;.<•   »?j))A«f   fcqdi))    Ag     lij, 


in  which    i^i,f  j)      is  any  point— *to  be   specified  more 
carefully  later— of   the    s<      -"r    ^ij"        Bu"   tbi-s  double    sus 
mey   clearly  be  regerdec"    es    the   sum   of    two  pert;.    ■*-  j  a     *■"''-  g 
defined   as   follows:    2^    ;    co    -    ■.    r    r"    "..  ;ose   terms    for  which 
i   =   J  i    *•"■■-  2   consists    -      all      h<     remaining  terms.        It   is 
easy  to    ,-ee    t'.  ot  :'        <    £    ,  cr    'rora    (23.8)    we   have 


lua 


(2o.  11)    ^giij    *    glSiJfri.Tj)    -   gU^,  r  j_i) 


Vit;        -a.-z.ir- 


*    \  f  L(7)Ku)KB(u,v)    du   dv; 


ana ,    by  a   v. oil   knov.M   integral  estimate*,    the    iterated   in- 
tegral   iii  the?    ri-.u  member  of    (23,11)    is   not    greater  than 
the    product  of  $*}   by  the  are*   or   the    square   Iji«        Since 
the    terns    in  <  i   are    those  which   correspond   to   the    squares 
I; it   we   need    only   compare    each   such   term  with   the   product 
of   N0    ty    by   the   ere  a    of    the    corresponding   square   I**    to 
show  that  £  2      is    rlot   greater   than    the   product   of     N0     <P? 
by  the   sun  of   the   areas   of   tbo   squares  In.        But   since 
this   latter   sun   is    less   than   €  /S^Nq^    it   follows   that 

(23.12)  Z.  x   <   £     , 

Now  let    €-*-()•        Then  from  the   definition   of    h 
it   follows   that  K  ->- ®Q .      i.oreover  from    (23.12)    we   have 

lire     S  x  -  o. 

But  then,  since  the  e  ■   -  Jtieltjes  integral  in  (23.7) 
is  the  limit,  as  h  ~>oo  ,  of  the  double  sum  (23.10}  and 


*  Courant :  Differential  ano"  Integral  Calculus,  Vol.  II, 
p.  232. 


lO? 

since   this  cloubl«  »u»n  is  Z\  ■>  ~<c   p    f    it  follows  th«t  the 
double  ."tieltjap   inte^rnl   S,n    lF-J#v)    is  *qu*l  to 

h-**» 

The   tenr,s  in  <r 2  b«in*  those  term?   in    (£3,10)   fpr  which 
I  t   $,    we    conclude    i'ro  i    (.';''. 9).     £23,10),    nnd    (53,11.)    th*t 
©11    t';ti    tor^is    Ln    2T  •    yc       ■  i  with,  the   exception   of    fcnose 
which   correspond   to    i  1^   lyinp;  entirely  within  the 

trian.??,e   8  piaiured   in   Firurs    ?<•        On  T.  the   integrand 
L(v)K(u}i     (v.,vi      ie    '.on   .  ■  i  o  is ;   so  the  theorem  of  the  zaesn 
assures  us   the   existence  oi    ,?    peir  of  numbers    ?  i,    r?  j 

sues:    ./»  t  CTi<-;L  i   |  L   «  (T-i,    Y,   x  2   ^  j   =   -C  end  for 

which    *'■■  -.a  equation 

(23.13)   ^.:i-L(?j)t;[i[i)K3(^i>^;ffi.q..J    lrr?i-X  1 

Ls    valid      lien     I^j      lies  entirolj    witnin    the   region  a. 
Const  qvieutiy   trie  double   Stieltjes    integral   in    (23.7)    is 
equal    to 

(25a4)    h^L  i?i  4  AfVL  >'  SaPi-ri-J   [rr*j-i  1  . 

with   A^     evaluated  for  the   arguments    i»?j,  ^C?j)    )»  *  f 
for    (li,    q(5i)    ),      L   for    (  *?  j ) ,   fc  for    (5l>»    a»a     Kjg  for 
(  €  i»  *?  j )  •        But    cue    limit    (23.14)    is    the    iterates   integral 
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in    (23.6),  which  is  therefor©  identical  witb  the  St;eltjta 
integral    (23,7),   as  was   to  he  proved. 

We  are  now  ready  to   combine  all   the  terms  in  tha 
ritfht  member  of    (23.3),    by  an  argument   quite   analogous  to 
that  wh,ich  was  employed   in  deducing    (22.2.8).        We  begin  by 
defining   for  a  given    i>    ,   and  a  given  ordered  pair   ( «  » /3  ) 
of  indices  on  disturbance   functions,   the  function 

(23.15)  g**X    (ai;rfi)-f    J      l/y(v)M*  (u)Xfe(u,v)    du  dv, 

in  which  L  ,(v),  M^  (u),  and  Kj;(u,v)   are  obtained  from 
(23.4),  (23.5),  and  (23.9),  respectively,  and  in  which  we 
again  return  to  the  summation  convention  for  the  repeated 
indioes  p,£,  j,  V,  and  i.   Evidently  (23.15)  is  the  sum 
of  ell  those  integrator  functions  of  the  type  (23.8)  de- 
termined by  the  aggregate  of  indices  p,  $  ,    J,  V,  i.   By 
the  theorem  on  addition  of  integrator  functions  for 
Stieltjes  integrals,  referred  to  in  the  remark  following 
(22.16),  it  follows  that,  for  each  ^  and  each  ordered 
pair  of  values  of  or  and  {$    ,  the  sum  of  the  corresponding 
terms  in  the  ri^ht  member  of  (23.3)  is  the  integral 

(23.16)  J  j  A(f    (?,q»)Af  (^,q(«r))dgvJ-J-   (si;<r,.r). 

But  the  right  member  of  (23.3)  is  the  sum  of  all  such  in- 
tegrals as  (23.16),  corresponding  to  all  possible  ordered 


XOl 


pairs    (<*  ,  i  )    for    °C-   1,    .    .    .    ,    j   ,     f*  *  1,    •    .    .    ,   n, 
'lence,    again  extending    the   summation   convention    to  the 
repented  indices   °<  ,  /3  ,   w*j  may  now  v;rite    (23.3)   as 


(£3.17)    nI][2(s)»  „  J  /:'U",qH)/lcf  (criq(cr))(jgv    I(si;r,T) 


In  this  last  equation,  the  region   S  is  the  square 
(si  *  *"  ■   s;  si  »  r  -   s)  and  the  integrator  .function 
<'*\'     (si;^,  t)  is  ;i        (23;lo).   A3  in  (22,22),  the 
ria:nt  member  or  (23.17)  3   the   i  of  m2  integrals,  one 
for  each  ord-r-  d  pair   s,  => }  . 


The  '   p  is  in   >pen,  of  course,  to  combine  those 
terms  occurrin  *  in  (23.22)  ana  (23.17)  which  correspond 
the  same  order  d  j  ir,(  <* ,  3 }  of  indices  on  disturbance 
functions.   3ut  thi   step  we  shall  postpone  until  after 
our  discu  ion  of  the  functions  M;   ~~"  (s)  in,  order  that 
we  iay  then  combine  all  beri   occurring  in  the  right  m©a- 
jr  of  (21.1)  which  correspond  to  the  ru-ix   pair  (ocy  /S). 


24.   The  Functions  /*-  p   (s). 


It  would  appear  from  (21.4)  that  the  functions 

TXT  T  %" 

/A.~t      (s)   are  considerably  simpler  than  either  6k.   ^(c) 

Is 
or  is-    °(s;,  in  tnet  only  a  single  quadrature  each  is 

necer;P3ry  for  their  evaluation.    In  vver,  to 

combine  those  terms  in  (21.4)  determinec  ;  i    p  rticular 

peir  (  <*• ,  /3  ) ,  of  indices  on  disturbance     stic  is  with 


10© 


other   similarly  determined   terms   fVom    (r,X%'£)    e  id    (/!).,.,, 
it   will    be    desirable    to   hove    the    function:'-     /*-    ^    (r>,    ex-» 
pressed    alp  >  double   Stieltjes    IntQUpals    ovf.r   t .:»•-    same 

interval  of   into  :ration  an   that  employed   for    (21.;;.    and 
(21.3).         phis    c-'iu    be  done,    though  sOMjuwhai;   ^roster   care 
is  necessary   in    the  process. 

In  f."Ct,    if  we  let 

(24.1)       *^(*j   =  Ap   It)  ^•r^(^(t»  *(t)    )    i, 

then  (21.4)  ce.i  be  written  as 

(24. S,  uxll(s)2E      N  "(rUf  (t,q(t)Uf  U>q(-r)  )  <*«• 


y 


lut,  -ive  v  ,  to  each  ordered  pair  («*f/3)  o"  indices 
on  distur  -  ice  functions  there  correspond'  i:.  tie  rir;: 
r-ie.nber  of   24. ^^  a  term  of  the  form 

(£4,3;  a  si,s)=      II  ^   (r)Af  lt,q{«r)JA<f  ""(*,. q. (*))  d*. 

Therefore  our  next  objective  will  be  to  show  that  K(si,s) 
is  expres  ;ibie  as  the  double  stieltjes  integral 


«>  (0)  ^Lii 

(24.4}  J,   d*  (t,qW)Al/;   (*-,q(*-))  dg       (n^i?) 

1  C*;C/j; 
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in  which  th^  integrator  function  i§   ,t 


*  * 


V  =  rnln  (Cf  tr)  >  ■   : 

4.5)        g  *  "ll    (si;<r,^)   m     /  N     "    (v)    flv. 


To   this    end,   we   observe   first  that   for   each   6  >0 
we   can  fine    8   partition  oj1   the   square   3  such  that,   for 
every  re-partition,    the   corresponding  finite   sum  is  within 
£      of  the   integral    (24,4    .        Let     o<     <  <*; .  <•••  <   °C  . 

°  a.  i\ 

be  £kll  the  numbers  whi  sh  >ccur  either  as  ^"-  or  ■  a  Z- 
coord  in  ates  of  vertices  of  s-ub^iatervals  in  a  partition 
with  this  property.        Consider  the    intervals 

(24.6)        IijS{    \<r,*)    l^-^^i^H^^^} 

■'<-■  * 
for  i  ■  1,  .  .  .  ,  h,  and  j  ■  1,  .  .  .  ,  h.   For  each  one 
of  these  intervals,  either  at  most  one  cf  its  vertices  is 
on  the  diagonal  '^  ~   si  ■  <T  -  si,  or  else  itstdiagonal  lies 
on  the  diagonal  of  the  square  3.   In  fact ,  * if  i ■  y-   J,  the 
diagonal  f '-  si  *  c?"-  s\     has  at  most  one  point  in  '£■   ;; 
^  whereas j  for  each  i,  In  is  a  sq.up.re  with  its  diagonal  on 
the  diagonal  of  j. 

Nov;,  in  order  to  simplify  the  notation  throughout 
one  remainder  of  this  section,  we  shall  employ  the  abbre- 
viations 


w 


vi:. 


(84. ?>   < 


HeoeUinr  that,  by  dofinlHon, 

(24. S)    ^  J    ^^f  <r|Q(r))A'-f    p,q-(«-)-)   A^{.54;  '•-,'-/ 


v/herein   (.f  iff?j)    is  any  point  of  the   interval  I^j,   sind 
xhst 

(24,9)  A-gl^i   "  «(SI5  n»vj)    -  8!8l5^i»rj-li 

♦  sUi;ri-l»r.5-l}  "  g<Bl>*i-l  »rvn   » 

it  is  ea^y  to  see  that,  using  (24.5) ,  the  right  member  •-' 
(24#9)  it.  identically  zero  for  every  interval  Ij  4  for 
which  i  ^  J.   For  if  Iij  belongs  to  the  triangle 
S\' {     (ft  v )  I  si  ■  F*  ■"  "£  ■  s  }  9  then  the  smeller  co- 
ordinate at  each  of  its  vertice:  is  the  <*.orr?Bpan(Un.«$  #*» 
so  that  the  terms  in   the  right  member  of  (24,9)  cancel 
out  by  p^irs.   similarly,  if  X j j  belongs  to  the  triangle 
32*.  {  (?*•?•>  |  3i  ■  r  t  r  ■  a  j  ,  then  the  anally/  co- 
ordinate at  each  of  its  vertioes  ic  the  corresponding  -"  » 
an4  again  x.he  terms  on  the  right  in  (24. 9)  cancel  in  pairs. 


1Q9 


Hew :  "•<*  may!  confine  our  attention  only  to  those  tcris  of 
the  double  s;r.  In  (£4:,S)  for  Which  i  »  j , 

ff.j't  from  (£4.5),  (24.7),  and  (24.9)  we  have 
(24.10)  ^giii  =:  ;   ^(v)  dv  -  /     N{v)  dv 


N'(v)  dv  •  /      Nlv)  dy 


"i  *i 


N(t)  dv. 


Moreover,    since    .     .      is   continuous,    the  mean   value   theore  . 
of  the  integral   6a~>i  cuius  guarantees  the  existence  of  a 
nuaiber   ft  i  with  **  j  „__  "  /^  i  •  ■■*  i     end  such  that 

(24tll)        /  rf'v;    dv   -  N(/Si)     [«Ci-«<i.i]    . 

"*•       A 
fc-± 

If  we  now  choose  booh    ^  j    and      l7i  tc  he  ft  i,    it  follov/s 
from    (24.3),    (24.10),    and   (24.11),   and  from  our   conclu- 
sions concerning  those  terms   of  the  right  member  o:?    (24.8) 
for  which  i  4  j,   that 

r  r  ex)  f0J 

(24.12)       J  J      Af    ltr,q(t))Af      (r|Q»)    d3    (si;^,"c) 

h  <«*;  ((3) 

-  Urn      5l   A<p    ^i»q(pi))Af  (fli,    qlft)  )N(/H>  fcr«i..i] 


no 


But  the  xext  member  of  (24.12)  is  the  double,  stieltjts 
integral  (24.4);  while  the  right  member  of  (34.12),  be- 
cause °<  1^1   »  ft  £  »  ^  i  end  si  ■  *  o  ^  ••«  <of  h  ■  s,  is 
the  integral  £lsi,s)  of  (24,3) •   Thus  (24.12)  establishes 
the  conclusion  we  set  out  to  prove  conoerning  that  inte- 
gral in  the  right  member  of  (24 .IS)  which  corresponds  to  a 
particular  ordered  pair  of  indices  &  ,   />   on  disturbance 
functions  and  a  given  index  ^  .  If  we  extend  the  summation 
convention  t,o  the  repeated  indicos  of  9   (3    in  (24.4),  we 
obtain   [  Cf. (24.2)  ] 

(H4.13)/iv^I(s)-//^f7^q(r))Af  (^.q(»-))dg^II(.1;<r,r)  ^ 

in  which  the  ri^iit  member  is  the   sum  of     nfi     integrals 
like    (24.4),   one  for  eech  ordered  pair   (  <*>(3)    » 
(<*»  !,..♦,   in;     ^3  »   1,    .    .    .    ,   m). 


25.   The  Functions  ^-^(s). 


From  (21.5),  (20.12)  and  (20.1)  we  have 


(25.1)  /^J'U)  «  2/AIJ1(s)  +  2/**V2(s), 


in  which 

(25.2)  At  ^(sjs  /    JC  (siiP^K)^^^  (r,q(-fc))  dr, 


Ill 


eju 


(25.3)/i,%f   I  1^*)/<^4f 4i(*,V«>)Af  «W»JWr. 


?-  5    r:  S 


wherein,   in  turn 


(25.4)       Pv.     (?)   3    Xl,   [t)  Fp    (r,o;{r))   jf    (or) 


and      (or,  (23.5)  1 


(25.5.        M        (<r)    9     v\i    («")    )■'.'  ,   (  <r9    oj*-}    ). 


If  we  now   c omp ere    (25.3)     l&th    (23.3),     v     se 
that  each  integral  in  the  right  member  of   the  former  eo.'1-*'- 
tion  hes  the  seine  for;1:  as  do  those    In  the  latter:    camels 
that  exhibited   Ln    (23,6) — provided,   of  course,    that  we   re~ 
place    A  <f  (t,q(r)]    in    (23a6j    by    /,  'i'r    •   ■  -r    c(tr)    }      f  r  each 
j    (J   ■  1,    .    .    .    ,   a)  .        r'1*"    bl  an   the   argj        a    employed  in 
establishing    (23„16)    aid    (23,17)    leads  ua  to  conclude  that, 
gi  '  ,   for  eac.i   ordered  pair    (  -'"  ,  fJ  )    of  iacli   es   on   dis- 

turbs..:.a   functions,    and   for  eech   j?    the   -    ~   3f  the   corres- 
ponding tera^  in  the  right  member  of    (£5.  0,    is   the  double 
Stielt.^es  integral 


US 


re      (U)  ((*}  yiv 

(B5,u)J,J    A<f  fj,(-r»q(^))  A<f     (er,q(0  )dg  K  (»i;rfr) 


wherein  the   Integrator  function     a  ^^    (si;**,  *  )      Is  de- 
fined    C  Cf.  (S3. 5)  ]    as 


Vtr 


u.»( 


\?iv  f   r      •>>       v 

(25.7)    gv  *    (8i;tr,r)s|  P    .     (v)h      (uifc^u.vjdu  dv   . 


V*5     «**  5i 


The   functions   in  the  integrand  of    (25.7)   ere   those  defined 
in    (25.4),    (25.5),    and    (23.9).        To  obtain  the  sum  of   those 
terms   in  the  right  member  of    (£5.3)   which  correspond  to  a 
given   v>    and   an  ordered  pair  of   indices  <*  ,  /S    on  disturb- 
ance functions,      we   sum  on  the    index  j   in  the  expression 
(25.6),    thus   obtaining,    say, 

(25.8)u      cs>')J    Af     .(z*<*>)Acp    C^ICrO  d%         Cs^r,*) 


Unfortunately  it  does  not  appear  possible,  in 

general,  to  combine  these  n  integrals  in  (25.8)  into  a 

single  term,  at  least  not  by  the  device  of  adding  inte- 

oc 
grator  functions.   7or  the  functions  A^nj  (  t  ,  q(-*)) 

are  conceivably  different  for  different  values  of  the  in- 
dex j.   jlxtendin g  the  summation  convention  to  the  indices 
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c*  ,  (3   ,   however,   it  follows  from   (25.8)   that  we  can  write 
(25.3)   as 

(25.9)>cIV2(s)./f4f<*(*,q(T))Af  {*,*(»)  )*&**!   (si;<r,*), 

In  which  the  right  member  is  the  sum  of  nra2  double 
Stieltjes  integrals,  n  such  integrals  for  each  ordered 
pair  of  indices  « ,  ft  on  disturbance  functions. 

On  the  other  hand,  for  trajectories  on  which  alt. 
tude  y  can  be  used  as  independent  variable,  and  if   the 
disturbance  Atf   depends  only  on  altitude,  then  A^fq^  3  0 
and  We  do  not  have  to  contend  with  tne  functions  (25,1). 

The  functions^  *V1  (a)  defined  in  (25.2)  involve 
a  single  quadrature  for  each  term  of  the  right  member  of 
that  equation,  and  so  are  formally  lilce  the  functions 
14.  ^**(s)  which  have  already  been  discussed.   They  will 
be  treated  in  the  next  section. 


£6 


?..     The  Functions  Ayl(s),   Aw2(s),   /*  IJLl(s),   and 


A,Ivl(a>. 


In  this   section  we   shall  attempt  to  show  how  the 
four  types  of  functions, /^(s) ,    A*2(s),  /*-^(s)   and 
//cIv^(s)    can   be  reduced  to   double  atieltjes  integrals  by 
the   sorae  methods  previously  developed  for  the  functions 


/c  ^(s)  and  A  IJI    (s). 
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We  shall  consider  first  the  functions   "  Cf.  the 
first  term  of  the  right  member  of  (22.3)  2 

(26.1)  >?■(•)■   I    ty.V'(r)^(r)SK(*c)  K    (si)Kr(si)    dtr  , 

•J 

in  which  the  abbreviation    If.    (*■)      is  defined  in    (22.1). 

Let 


9 


*. 


(26.X)       q,  \{r)  «    Y".  <*>*,   M  J    "<*> 

so  that    (26.1)    can  be  written  as 

(26.3)    Ac    ^(s)    2  /  KT    Ul)  K  (si)    fcf     (r)    d-r  . 


T*  SL 


Then,    if  we  define     g^    ^si;*-,*)      by    £  Cf .    (24.5)  ] 

\f^  win  C^V** ) 

(26.4)  g^/tsiJC.t)    i/  <4*M    dr» 

an  argument   similar  to  that  which  established  the  equiva- 
lence of    (£4.3)    and    (24.4)    yields  the    conclusion  that,    for 
each  ordered  pair  of    velues  for   the    indices     V  ,    2    ,   the 

sum  of  the   corresponding  terms  in   the   right  member  of 
(26.1)    is 

(26.5)  K     (si)    K   (si)    JS    dg'*11        (si;«r,  <c)     . 

<4>  CSi  s  CYJCS) 
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Vflben  we  extend  the  summation  convention  to    Vf    S    ,   it 
follewa  from  (26,5)     that   (£6*3)  can  pe  written  as 

(33.0)/*-  J1(s)»)C(si)H:^i)//d^(ei;<r,r),    [V,/i«;,  .  ,,,nj. 

Secondly,  we   observe  '/net   the    seccad   ten.j  of  the 
right  member  of   (2£.3)    can  be  written  as 

(26.7)  1*  ?(•)*£  J     /     V  *V;f' (*-)?*(*)  It  (3i)i   (rU'Vi    drdtr. 


r*\  r**k 


JJow  let 


(86.8)  y£m  ■  y&mt^m  r^w  » 

and    [  recalling   (23.f>)   ] 

(26.9)  }.*    (<r)    e   )/(<r)   F\  (u>(<r,    q»))    . 

Then  (26.7)  can  be  written,  using  (20.1),  (26. 8),  and 
(26.9) ,  in  the  form 

(86.10)/*,  J8(s)*2  f  J       tyjsMX*(r)KJsi)ft<f   (<r,q(r)  )d<rdir. 

Comparison  of    (26,10)   with    (23.3)    shows   that,    formally, 

(26.10)  is  a   special    oase  of    (23.3)   with  L     .  ( t  )    renlnoed 
by    "^     (  t)   and  Af  (v,    q(*"))    replaced  by      K"    (si).       Hence 
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by  the   argument  which  established    (23.16),   we   conclude 
that   for  each    v>  and  each  ordered  pair  of    indices    (-/,/?) 
the   corresponding  term  in  the   right  member  of    (26.10} 
is   the   double   jtieltjes   integral 

(86.11)  K^(S1)  //  4(f    ('»q(r))d8Cy;I<%     tsu^r) 

S'  ' 

in  which  the  integrator  function  is 

(25.12)  g*1*  (si;^r)s  /    f    T  V)  ii-f(ff-)  d^r  dr. 

C3ee  (23.5).  ] 

Evidently,  when  we  extend  the  summation  convention  also  to 
*y  ,  /3  ,  it  follows  from  (26.11)  that  (26.10)  can  be 
written  as 


(26.13)  /a,   l2(s)=  K   (si)  ff'Af  (<r,q(<r))  dgvI2(si;  «r ,  *  ) . 


As  we  observed  in  #  23,  the  functions  A    (s) 
defined  by  (23.2)  are,  formally,  like  those  designated 
/<-  J  (s)  in  (21.4),  in  that  only  a  single  quadrature  is 
involved.   3y  comparing  (23,2)  with  (24.2),  which  is 
equivalent  to  (21.4),  we  see  that  (23.2)  is  a  special  case 
of  (24.2)  with  L  V  (r)  in  place  of  IT  "  (r)  and  with 
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0 
If  (si)   replacing  Af   (t,   q.{*)    ).       Thus,   If  we  Refine 

[  see    (24,5)  ] 


V:  rr»i'nC<r,rJ 
(26.14)      *VlIl(si;  o-,  «-)■    /  -  L*\(v)    dv   , 

it  follows  by  the  argument  which  established  the  equiva- 
lence of  (24.3)  and  (24.4)  that,  for  each  ordered  pair  of 
the  indices  (  °<  ,  V)  the  corresponding  term  in  the  right 
member  of  (23.2)  is  the  double  Stieltjes  integral 


(26.15)    rc  (si)  ff   ^'"hr.qfcndg*111  (si;<r,tr). 

Furthermore,  extension  of  the  summation  convention  to 
c/  ,  «/  enables  us,  therefore,  to  write  (23.2)  in  the  form 


(26.16)/AIJ1(s)»  KUl)ffAf   (•r,q(r))dg>'111  (si;<T,<*) 


IVi 

Finally,  we  note  that  the  functions  /&*- ^     (s) 

III 
defined  in  (25.2)  are  formally  like  the  functions yU.  ^    (s) 

in  (24.2).   Given  v>  ,  for  each  j  ( j  -  1,  .  .  .  ,  n) ,  the 

function  P  .  .  (t)    in  (25.2)  corresponds  to  n' ,*  (•*)  in 

«  i  Y  «  ff 

(24.2),  the  function  d  <f  j  (*,q(ir)  )  to  Af    (t,  qM),  and 
the  constant  VC.(sl)  to  the  function  A  cf  (r,  q(*«)).   Conse- 
quently, the  argument  by  which  we  showed  that  (24.3)  is 
expressible  as  the  double  stieltjes  integral  (24.4)  enables 
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as  to  conclude  that ,  given  >> ,  for  each  ordered  pair  of 
the  indices  *,  Vt   and  for  each  j ,  the  sum  of  the  corres- 
ponding terms  in  the  right  member  of  (25.2)  is  the  double 


Stieltjes  integral 


(86.17)  Kw/V  foyi*>VV>  J  CZlCh  '  ''  ^ 
in  whiJh  tl^e  integrator  function  is  [see  (24.5)] 

(26.18)  *f      V*;*;*»     /  rp*    coi/ 

Summing  on  the  index  j  In  (26.17)  gives  us,  as  the  sum  of 
those  terms  in  the  right  member  of  (25.2)  and  an  ordered 
pair  of  indioes  (*,  V),  the  expression 

<26.19)  A"v  ^,  ^cvJ/^;rM.w>j^.r;^ 

The  remarks  following  (25.2)   apply  equally  well  to   (26.19).. 
Summing  on  *,V    in  (26.19)   enables  us  to  write  (25.2)   as 


(26.20)^     tut  K/^JJ"  Af  *.  (r,$U))  i^.  ICii^Oj 


V  t 


in  which  the  right  member  is  the  sum  of  nm*"  double  Stieltjes 
integrals,  n  such  integrals  for  eaoh  ordered  pair  of  the 
indioes  «f,  V. 


1X9 

37.      The   vu  notions  /^jjjjiL^I^  re  §L°ii^i^-  t  ?  6 , . Fa  i  r 


Atf*  >  Atf  ^  °f„I)ist^ban£fej;  ans. 

To  comoute  the  second  differential  effect  on  a 
trajectory  element  -  say,  on  range  -  due  to  a   disturbance 
(A     q.o,Af  ),   we  have  seen   that  one  or  more  of  the  func- 
tions x  i  defined  by  (16 .1)  must  be  calculated.   Jn 
(£0.15)  we  have  seen,  furthermore,  that  each  second,  varia- 
tion  3:  i   is  expressible  in  terms  of  the  n  functions 
^  \t    „    .  .  ,  f   n   appearing  in  the   ith  row  of  the  de- 
terminant  (20.8)   and  the  n  functions 
/CM  (s),  J>«  1.  '.•••'.  n]   defined  by  (20.13).   In  if 
21-26  we  have  seen  how  each  term  in  any  one  of  the  func- 
tions A>(s)  is  expressible  as  a  double  Stieltjes  inte- 
gral.   Evidently  we  must  add  these  integrals  to  obtai 
for  ']    given  ->   f    the  corresponding/^.    It  is  fortunate 
that  all  these  integrals  have  the  same  region  of  integra- 
tion— the  square  designated   J  in  f   22  and  3  in  9  23. 
Moreover,  many  of  these  integrals  nave  the  same  integrands 
but  different  integrator  functions;  we  foi     at  we  could 
combine  those     -rais  having  this  property  into  a  single 
integral.     t  this  combining  process  is  only  partially 
finished;  it  is  our  purpose  in  the  present  section  to  com- 
plete it. 

In  order  to  help  us  fix  in  mind  what  needs  to 
be  done,  let  us  suppose,  for  example,  that  we  wish  to 


120 


obtain  the   second  differential  effect  on  range  of  the  com- 
bined presence  of  F  ren<?e-wlnd     wx(y)   and  a  non-standard 
temperature 0U) .        Let  us  further  suppose  that,    among  the 
fourteen  disturbance  functions  jAq>    ,....,   &<{      $      all  are 
identically  zero  except  those  which  represent     wx(y)    and 

0  (y)  ;   say  these   are    h<f  l"'t    &<{  9 ',   respectively.        For 
simplicity,  we   shall  assume   also   that   the  initial  condi- 
tions are   standard,    so  that    A  q0  -  0,    for  i   »  1,..«,6. 
Then,  for  eaoh  value   of  v>  (  \>    »   1,...,6),   the   corresponding 
function  /^^  will   be   the    sum  of  ell  those   terms  in  the 
right  member  of    (21.1)   which  contain   &f    ,   or    ^  f  9,   or 
both  Af®  and  ^f9,    or  the   derivatives    &fyt  ACf*,    in 
their  integrands.       Such  terms  will  be  those  whose  inte- 
grands  contain  the  products  AfQ6f9t    (Af8)2,    (4<f9)2   , 

4({>0   Af  J,    ^tf^f9,    A(f8Ap8,   and    A  <f>  %  <p  9. 

To  these  products  there  correspond  the  ordered  pairs 
(8,9),  (8,8),  (9,8),  (9,9)  of  indices  <*  ,  (I     on  the  dis- 
turbance functions*   For  example,  the  terms  which  corres«- 
pond  to  the  ordered  pair  18,9)  consist  of  those  containing 

a   8      9 
the  product  »*f  y  ALf>  .     and  half  of  those  containing  the 

product  A<f  Q   A</>  9;  the  other  half  of  the  latter  terms, 
together  with  all  the  terms  containing  the  product 
^f  y  **f        correspond  to  the  ordered  pair  (9,8)  ♦   The 
terms  containing  (4f8)2  and  A  if  8  A^f  8  are  those  which 
correspond  to  the  ordered  pair  (8,8),  etc.   Clearly  we 
shall  obtain  all  the  terms  which  occur  in  the  function  A4, ^ 


J^Cj  ju 


if  we  add  t<        1    '<  Lie ct Ions  of   terras  which 

correspond  to  the  :d  pairs  of  Indices  (8,9), 

(8,3),  (9,8),  ana"  (9,  y        spectivel:  s  can  be  most 

easily  accomplished  :  (  find  an  ex  ression  for  the 

sum  of  all  the  terms  right  member  of  (2Q,13)  which 

correspond  to  a  given  oJ  °  and  a  given  ordered  pair 

(  <<  ,  /3  }  of  indice  l\  — -b-      .  r.ions,  ini 

conditions  assumed  st  d, 

i  we  seek  is,  thus  the  sum  of  (22.21), 
(23.16),  (84.4)  ana  (  i.8).   [Since  the  initial  conditi 
are  normal*  the       ties  K"  (si),  IfC..  (s\) ,    defined  by 
(20o12)  and  occurring    the  functions  considered  in  $  2ov 
are  all  zero.   Hence    e  of  the  functions  /^J-^lsj,  etc., 
discussed  in  f  26  occ    in  the  expression  we  desire.j 
oince  the  integrand  in  each  of  the  expressions  (22.21), 

(at)  (?> 

(£3.16),  and  (25.4)  is  6  f       (r^(r))  A<p  (<T,q(<r)   and  the 
region  of  integration  in  each  is  the  square  &,    the  sum  of 
these  three  expressions  is 

(27.1)  J  )  A(f   (r,cjr))  Af      (er,q(c-))dg  (M)  (.}  (s.si;*-,  •*  ) 
in  which  the  integrator  function  g    ~j  (s,si;  f, t  )  is 

(27.2)  g   ,  *s,si;<r,r)=g  vI   +  g»IX  +   S*111. 

The  integrands  of  'che  n   integrals  in  the  right  member  of 
(25-8)  are  dependent  upon  j  and  are  conceivably  diff. 
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for  different  values  of  j.   Thus,  in  general,  we  are  not 
able  to  combine  the  integral  in  (27.1)  with  those  in  (25.8), 
to  obtain  in  the  form  of  one  double  Stieltjes  integral  the 
sum  of  all  the  terms  in  the  right  member  of  (21.1)  which 
correspond  to  an  ordered  pair  of  indices.   But,  as  we 
observed  in  the  remark  following  (25.9),  for  trajectories 
along  which  it  is  possible  to  use  altitude  y  as  inde- 
pendent variable  and  when  the  disturbances  depend  only 
upon  altitude,  then  the  integrals  (25e8)  are  identically 
zero. 

Consequently,  in  the  case  of  the  illustrative 
example  considered  earlier  in  this  section  (involving  the 
two  disturbance  functions  Af 8,  A  <f>  9  dependent  on   y 
alone  along  a  one-branched  trajectory,  with  standard 
initial  conditions),  for  each  »>  (  )>  *  1,  .  .  .  ,  6)  the 
function  /*-^   (s)  is  the    sura  of  the  four  double  Stieltjes 
integrals  obtained  from  (27.1)  by  giving  (  <*  ,  (3  )     the 
pairs  of  values  (8,8),  (8,9),  (9,8),  and  (9,9).    Of 
course,  the  two  integrals  corresponding  to  the  second  and 
third  of  these  pairs  of  indices  are  identical. 

If  initial  conditions  are  not  standard,  then  the 
integrals  considered  in  }    26  must  also  be  computed,  and 
added  to  any  obtained  by  the  use  of  (27.1). 


V  ■  »  • 
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Chapter  VIX 

•        •         ••-.*•■       " 
J^cSHANijj'5  TWO^y.  OH  7.TE  QriliGT)  DIFFERENTIAL. 

"*   S*<*     '' '  •  '        .  •  • 

AN  ANALOGOUS  linaQRi ±    CN  TILS  SECOND  DIFFERENTIAL. 

28 .   One-b ranched  Tra jectorlen.   Uniqueness  of  solutions 
of  D iffere n t i al  ttquetions  of  Motion. 

Our  theorem  (15,10)  on  the  existence  of  the  seco.  ' 
differential  or  the  mapping  defined  by  (7.1),  (7.2)  is  a  re- 
sult of  the  sa.'ue  type  es  that  of  Bliss,  (11.3),  on  the  first 
differential,    iv.   si  . larity  arises  from  the  fact  that,  in 
both  theorems,  the  nor  of  the  disturbance  includes  the  de- 
rivative of  the  disturbance  functions  with  respect  to  the 
variable   y  when  the  disturbance  is  one  (such  as  wind,  ab- 
normal temperature,  or  abnormal  density)  which  is  assumed 
dependent  only  on  altitude  r   B        As  pointed  out  in  the 
paragraph  preceding  (11,3),  :..c3hane's  main  theorem  [(11,7) 
of  Qwcs]  3  ^s  8n  improvement  on  (11,3)  in  the  sense  that  it 
establishes*  the  existence  of?  the  first  differential  when 
Nly»  from  which,  the  objectionable  derivative  is  missing, 'is 
used  as  norm.  .  Our  'chief  objective  in' this  chapter  will  be 
to  establish  a  theorem  on  the  second  differential  analogous 
to  that  of-  ;,.cshane  on  the  first  differential, -and  thus  to 
sharpen  our j conclusion  (15.10)  in  a  similar  way  to  that  in 
which  I.cShane's  theorem  improves  on  that  of  Bliss, *   But  tto 
accomplish , this  objective  we  have 'to  introduce  a  new  norm, 
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which  can  be  larger  than  N$»y  s"d  whieh  may  not  reduce  to 
No  even  when  the  disturbance  depend*  ealy  on  altitude, 
yor  this  new  no-wi  includes  the  vaiues  of  fat  f Ira jj  partial 
derivatives  of  the  disturbance  tn  a  neighborhood  of  the 
aumittlt;  and  the  derivative  of,  say,  range~wind  wx(y)  with 
respect  to  y  ia  not  necessarily  zero  near  the  sucurdt. 

In  order  to  provide  the  background  necessary  for 
a tat lag  licShene's  the  :rem,  and  also  to  assemble  for  the 
reader1*  convenience  certain  material  from  [lies]  which  is 
essential  to  the  proof  of  cur  theorem  on  the  second  dif- 
ferential, we  shall  again  quote  several  passages  Fronj  that 
paper  in  thie  chapter .   We  begin  with  certain  remarks  on 
one -branched  trajectories,  leading  to  e  t  iecreiu  on  the  con- 
tinuity and  uniqueness  of  the  solutions  of  tha  differential 
equations  of  (notion  (  USi.   The  excerpts  whiah  we  quote 
in  this  section  ore  taken  from  $f7,  8,  9,  and  10  Qf  the 
manuscript  of  [kcs}a 

"In  |  1  (of  [iCos];  i  £  of  the  present  pap«r)  it 
was  pointed  out  that  for  t^e  purposes  of  ballistic*  it  -s 
unobjectionable  to  include,  in  the  definition  of  t&e  •'  '.» 
the  derivatives  of  the  disturbances  with  respect  to  t.u 
variables  vx,  vy,  vZt   x,  rad  ^\   it  is  tha  inclusion  ;Of  da* 
rivatives  with  respect  to  y  thct  we  desire  to  avoid*   En 
the  definition  of  N8  tha  derivatives  oi  f     wJt:  rsapeot  t^ 
the  qj  ore  present,  but  the  derivative  with  resyeet  to  a 
is  absent.   Thus  if  v,*e  identify  c  with  y  the  situation 
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is  what  we  are  seeking.   Since  s  ■  y,  the  norm  Ns  will  be 
written  es  Ny.   Suppose  then  that  we  are  interested  only  . 
in  the  subset  of  trajectories  with  Vy>  0.   On  these  we 
can  introduce  y  as  independent  variable  and  apply  the- 
orems" (7.8),  (10.1),  and  (11.3).   "In  fact,  if  we  are 
interested  only  in  the  effects  of  changes  of  initial  condi- 
tions, wind,  density,  and  temperature,  the  latter  being  re- 
garded as  a  function  of  y  alone,  4?  is  identically  zero 
and  A<f   is  independent  of  qJ  ,  so  Ny  m   Nq#   The  norm  Ny 
doe?  not  then  merely  lack  the  objectionable  property  of  N8 
when  s  is  t  or  x  or  anything  but  y;  it  has  coincided  with 
the  very  reasonable  norm  No.   The  same  situation  prevails 
if  we  regard  only   the  trajectories  with  vy  *  0;  these  in- 
clude dive-bombing  trajectories,  and  level-bombing  trajec- 
tories if  we  cut  off  the  short  time  at  the  top  during  which 
the  bomb  is  sheltered  in  the  bomb-bay. 

"However,  for  many  trajectories,  including  all 
ordinary  ground-to-ground  artillery  trajectories,  it  is  im- 
possible to  use  y  as  independent  variable,  because  the 
trajectory  has  an  ascending  and  a  descending  branch.   We 
can  use  the  time  t  as  independent  variable.   But  if  ^  , 
©,  wj,  wz  are  merely  continuous  functions  of  y,  the  stand- 
ard theorems  on  differential  equations  will  not  even  guar- 
antee the  uniqueness  of  the  solution.   Our  first  task, 
then,  is  to  show  that  even  when  ^  ,  etc.,  are  merely  contin- 
uous in  y,  the  equations  of  motion  are  uniquely  solvable.* 
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Ivi03hane's  proof  of  the  continuity  and  uniqueness 
of  the  solutions  of   the   system   (6.3)  makes   strong  use  q?  i» 
lemma,    (8.4)   of    [ttcs] ,   concerning  the  propagation  of  dis- 
turbances a  Ion/?  the   descending  branch  of  a   trajectory. 
This  sane    lemma  will  also  he  employed  in  our  proof  of  t)  . 
principal  theorem  of  this  chapter  on  the   second  differen- 
tial,       iie   state  it   in    -he  next  excerpt,  fr.oin  fislcsji   after 
first   quoting  some   introductory  passages. 

"Let  us  suppose-   that  we  desire   to  solve'*   (6»3) 
"with   certain  initial  position  and  velocity,  the  density  ?   s 
temperature  9  and  wind   components  wx,   wz  being  continuous 
except  perhaps   at   a  finite  number  of  values  of  y,   where   they 
may  have   simple   jump  discontinuities,   and  having  bounded 
partial  derivatives .with ' respect  to  x   and   z  continuous  ex- 
cept et   finitely  many   values  of  y.        If  vv0  >"0,    up  to  the 
summit  we   can  introduce     y     as   Independent  variable   in'* 
(5.3)    ''end  find    that  the   resulting  equations,  are  uniquely 
solvable.        Likewise   if  tv:o  trajectories  have   coincident  x, 
t,    z,    t-v,   vy,    vz  et   some  level  y  on  the   descending  branch, 
these   descending  branches  are  identical.       however,   these  , 
standard  theorems   do   not   at  once   inform  us   that  a   given 
ascending  branch  might  not  be   followed  by  two   distinct  de- 
cending  branches.        3y  continuity,    the  values  of  x,    y,    z, 
vx»    vy  "  Of    vz  are  uniquely  deter  nine d  at  the   summit  because 
the   ascending  branch  is  unique.        ,/e  must  show  that  these 
summital  values  determine   a   unique  descending  branch. 
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"Let  us  first  transfer:. 1  equations"  (6.3).    "sup- 
pose that  we  nre  given  the  descending  branch  of  a  trajec- 
tory, corresponding  to  tine?  t  in  an  interval  ftQjtg],  and 

defined  by  functions  x(t) Yz(t).   This  we  shall 

call  the  'original  trajectory'.   V/e  understand  that  t0 
corresponds  to  the  summit.  Let  x(t),  .  .  .  ,  Vz(t), 
t 3  ■  t  ■  t4,  be  another  arc  of  trajectory  (the  'second  tra- 
jectory1) with  v"y  <  0  and  with  y  between  y(to)  and  y(t2). 
To  each  time  t  between  £3  and  T4  there  corresponds  a 
time   s   between  t0  and  tg  for  which 

(28.1)         y(s)  =  y(t). 

If  in  particular     S3,    s4     correspond  to  ^3,   T±     respect- 
ively,   equation"    (28.1)    ^establishes  a   one-to-one   corres- 
pondence between  s  and    t,   which  we   can  write   in  either  of 
the   forms 

(2392i)  s  -    s(t)         (t3  -   t  =    t4    ), 

(£8.2ii)  t   =   t(s)         (33  »  s  =    S4    ). 

This  new  variable   s  we  introduce  es  independent  variable 
on  the  second  trajectory." 

There  follow  in  [L:cS]  the  details  of  the  trans- 
formation of  (6.3).   We  do  not  give  them  here  as  they  ere 
not  necessary  in  what  follows.   The  result  is  a  system  of 
equations  in  which  s   is  independent  variable,  from  which 
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y  is  absent i  find  which  is  satisfied  by  the  ".-second  trajec- 
tory",   [j./hen  y  is  nodded,  It  can  be  obtained  from1  (:.8,,1).J 
By  introducing  the  nynsbQls  q}t    .  .  ,  ,  q'J     for  the  w  l-. 
ables  vy,  Vy,  v„,,  k,  7>»    t  in  that  order,  this  system  can 
be  written  in  a  fom  analogous  to  (7,1),  and  the  "original 
trajectory"  satisfies  e  pinilar  system.   In  terras  of  the 
new  independent  variable  s  and  the  new  dependent  variables 
qi  defined  above  we  can  now  state  '    ian<  "   lemms  pr  the 
propagation  of *  disturbances  along  the  .  seen&ing  branch. 

(28.3)  Lemma  £  (3,4)  of  [,:c3]  J:  '  Corresponding  tc  the  de- 
scending branch   aj.  -  cj/Ms),  (to  *  s  =  tg) ,  of  the 
original  trajectory  there  is  u  constant  I-Io  with  the 
following  property.    I?  qi(s)   and  q.2(s), 
(S3  =  s  =*  S4) ,  are  descending  arcs  of  trajectories 
with  tQ  ■<•  S3  *  S4  ■  tg,  both  being  computed  with 
the  same  functions  F*,  ay,  ¥>  **  but  with  possibly 
different  initial  values,  and  if 

(28.31)    q^(s)/«l1(s)  -1/2,  q^(s)/ql(s)    =  1/2, 

(S3  ■  s  «  s4) t 

then 

(28.311)    lq2(a)-o1(si|  =  K^  |q2(s3)-q1 (s3) | , 

(S5  -  s  -  S4). 

TA'cShane  •  used   the   lemma    (23.3)    to   show  the   urn  r.     - 
nees  o-"  the   solutions   of    (fi.3)  J    i.e.,   that  two  trajectories 
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satisfying  the  system  (6.3)  and  having  the  same  initial 

values  Vy0  ■  0,  vxo,  vzo,  *o»  zo*   t*Oi  Yo  aro  identical. 
Together  with  (7.3)  it  was  also  used  to  prove  the  continuity 
of  the  solutions  of  (6,3)  as  functions  of  the  disturbances 
(  fr  10  of  [lcs] ) .   Lore  specifically,  the  latter  property 
was  established  in  the  sense  described  in  the  following 
passage  from  f  10  of  [lcs]  : 

"Now  we  assume : 

(28.4)  The  function  is  is  a  continuously  differentiable 
function  of  x,  y,  z,  vx,  vy,  vz>  wx,  w2,  ^  ,  6 

regarded  as  its  independent  variables. 

(28.5)  The  functions  wx,  wz,  ^  ,  6,  ax,  ay,  az,   and 
their  first-  order  partial  derivatives  with  re- 
spect to  all  variables  except  y,  are  bounded  and 
are  continuous  except  perhaps  at  a  finite  number 
of  values  of  y,  where  they  may  have  simple  jump 
discontinuities. 

(28.6)  "Let  (28.4)  and  (28.5)  be  satisfied.   Let  the 
original  initial  values,  and  the  wind,  density, 
temperature,  and  added  accelerations  be  such  that 
the  solution  of  (6.3)  lies  in  Ro  for  ti  ■  t  ■  t2. 
Then  for  all  other  initial  values,  etc.,  satisfy- 
ing (28.4)  and  (23,5)  and  lying  sufficiently  close 
to  the  original  values  in  the  sense  of  the  norm  No, 
equations  (6.3)  will  again  have  a  unique  solution, 
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and  this  will  lie  in  R0  for  t^  -  t  «  tg,  and  as  i;0 
tends  to  zero  the  change  in  the  solutions  (with  the 
usual  norm  (7.6) )  will  also  tend  to  zero.   In  this 
theorem  we  may  replace  equations  (6.3)  by  any 
equivalent  system." 

29.  KcShane's  Theorem  on  the  First  Differential, 

In  this  section,  we  shall  state  McShane's  theorem 
[  (11.7)  of  [KcsJ  ]  on  the  existence  of  the  first  differen- 
tial when  the  norm  of  the  disturbance  is  Niy  (Nis  when 
s  «  y) ;  wo  shall  also  describe  in  detail  the  new  coordinate 
system  introduced  by  McShane  in  the  proof  of  his  theorem, 
for  ve  shall  use  it  in  the  proof  of  our  theorem  on  the  sec- 
ond differential  in  £  31.   Quotations  occurring  below  are 
from  the  original  manuscript  of  [lies]  • 

We  shall  assume  that  we  have  given  a  two-branched 
trajectory  T  satisfying  the  system  (6.3).   We  suppose  that 
T  has  a  summit  at  the  point  which  corresponds  to  a  value 
t0  of  t  between  ti  and  t2.   As  an  additional  hypothesis 
we  also  assume : 

(29.1)  "On  a  neighborhood  of  the  summita]  values 

(x(t0),  y(t0),  .  .  .  ,  v2(t0))   all  eleven  func- 
tions cp  are  continuous  end  continuously  differ- 
entiable ;  and  the  F*  are  everywhere  twice  con- 
tinuously dii'ferentiable. 
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•In  praotics  this  hypothesis  is  quota  JiCirraleaa. 
The  'normal'  trajectory,  from  which  disturbances  r^  mees-. 
ured,  is  always  computed  with  constant  wx  end  w<j  [ordta- 
arily  both  are  identically  zero)  anc"  with  ^>  and  9  de- 
pendent on  y  alone,  and  either  analytic  or  at  worat  com- 
posed of  two  analytic  arcs,  meeting  (with  discontinuous  de- 
rivatives) at  the  lowest  altitude  of  the  stratosphere." 

Although  T  is  a  two "branched  trajectory,  we  use 
the  norm  appropriate  for  one- branched  trajectories.   This 
fact  cannot  be  made  too  emphatic,  for  the  crux  of  tic  Shane '  s 
theorem  is  that  the  norm  Nly  can  be  used  in  the  proof  of 
the  existence  of  a  first  differential  on  a  two- branched 
trajectory;  we  are  not  limited  to  the  use  of  the  less  de- 
sirable norm  N]_t,  despite  our  inability  to  use  y  as  in- 
dependent variable  along  a  two-branched  trajectory.    (As 
an  "extra  dividend",  of  course,  the  norm  Niy  reduces  to 
the  more  reasonable  norm  No  when  the  disturbance  depends 
only  on  altitude.) 

Suppose  then  that  we  have  an  original  trajectory 
T  lying  in  R  end  satisfying  (29.1).   If  Ax0,  A  y0,  a  z0l 
A  vx0»  L   vyo,  ^vz0,  AtD,  4  f  ,  4$>    A  wx,  A  wz  are  dis- 
turbances with  Sufficiently  small  norm  No,  then  the  equa- 
tions of  motion  with  initial  conditions  xo  =xo  +  A  xo» 
etc,,  with  density,  etc.,  ^  =  f  +  a^>  ,  9  =  9  +  a  6  , 
wx  =  wx  +  A  wx,  wz  =  wz  +  A  wz,  will  have  solutions  like- 
wise defined  on  [ti,  t£\  .    ;,e  now  choose  a  value  t3  between 
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t0  and  tg,  and  between  y(ts)  end  yUg)  change  to  y  as  in- 
dependent variable.   V/e  find  that  ailon^  the  new  trajectory 
f  the  veriables  x,  etc.,  differ  from  the  original  vplues 
x,  etc.,  by  the  quantities 

4i(y)  -  x(y)  -  x(y),      A  z(y)  »  z(y)  -  z(y)  , 
At;y)  =  T(y)    -   t(y),      Avx(y)  «  vx(y)  -  vx(y), 
Avy(y)  =  vy(y)  -  vy(y),    Avz(y)  »  vz(y)  -  vz(y). 

ive  now   state  L,, Share's  theorem,  on  the   existence 
of   a  first  differential. 

[23.2]    Theorem    (^cshane)  ;      Let  T  be   an  original  trajectory 
ou    which    (29,1)    is  satisfied.       Let   the    summit   cor- 
respond to  t0l    and  lot  t3  be   any  number  such  that 
t0-^t3-<t2   .  Then  to  each   set  of   disturbances 

Ax0,    •    •     .    ,   A  wz     there    corresponds  a   set  of  six 
functions 

Y?  i(y)  =  Y?   i(y;  AXQ)   .  .   ,  ,  AWZ) 

(i  -  1,  .  .  .  ,  6),  y(t2)  -  y  »  yitg),   with  the 
following  properties.   They  are  continuous  on 
[y("t2)>  y(t3)"J.    /ith  the  usual  norm  (7.6),  they 
are  linear  functions  of  the  disturbances,  and  are 
continuous  not  only  when  the  disturbances  are'  meas- 
ured by  Niy,  but  even  when  they  are  measured  by 
N0.    if  we  assign  vy,  vx,  vz,  x,  z,  t  the  alter- 
native symbols  q1,  .  .  .  ,  q*3,  respectively,  *J>  (y) 


is  the  dlt 

S  Ct  3 

0>      • 

{  •  y  **  :    ) 

te   sy--  s    in 

proving    (29.2)   is  described  in  the  'following  pai 
quoted  x  he  manuscript  of    (} 

f«  ^ 
•      chosen   in   #6)    "are   translated   horizontally  so 
s   on   the   y- 
y-axis  so  that  at  the   summit   the   velocity  vector- is- 
(vx(t0)*   0,0) ,   with  vx{t0)  ">0.. 
as  b   the   trajectory  dc 

except  at   the   summit.     .  This   is   com 
sential,   hut   it  hardly  seems  worth  compile 

merelj  :lude  differential 

of  "boomerang®. 

"Now  we  construct   i  in 

(x,  -y,    zf    Vx,    vy,    vzJ-spfjoe* 
unchanged;  the   first  two  are   ir- 
es follows.        .   -  choose  a   posil  ding  t 
"clearance"')    and    in   the   yz-plan 
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line  at  altitude  y(t0)-fc.   This  line  is  tha  upper  boundary 
of  the  halft-plane 

(29,3)   W;   x  *  0,   y  »  y(t0)  ~»t. 

Let  V   be  a  point  not  on  tit    and  let  (x,  y,  z,  vx>  Vy,  vz)  be 
its  coordinates,   './e  define  p  to  be  the  distance  from  i  to 
;/,   If  y  =  y(to)  -  lc  we  pass  through  P  a  horizontal  half- 
plane  bounded  by  J.   If  y>y(t0)  ~*t  we  pass  through  P  a 
half-plane  bounded  by  the  upper  boundary  of  W.   In  either 
case  the  half-plane  through  i-  intersects  the  original  tra- 
jectory in  exactly  one  point;  this  is  rather  obvious,  and 
will  be  proved  incidentally  in  the  sentences  following 
(29.8).   To  this  point  corresponds  a  value  of  the  time  on 
the  original  trajectory.   Thi3  time  we  denote  by  the  symbol 
s.   Thus  if  y  =  y(to)  ~  hr  and  x  <   0,  to  ?  there  corres- 
ponds one  velue  of  s,  and  this  s  is  uniquely  -Jeter  lined  by 
y,  end  is  less  than  t0.    If  y  =  y(t0)  -  k  and  x  >  0,  to  I 
there  corresponds  one  value  of  s,  and  this  s  is  uniquely 
determined  by  y,  and  is  Greater  than  t0.    ..here 
y  >y(t0)  -  fC  ,  the  bwo  'variables  s,p  ere  uniquely  deter- 
mined by  (x,y;  and  conversely;  hut  y  alone  does  not  deter- 
mine s f    nor  does  s  alone  determine  y. 

•There  is  a  tirae  at  which  y(t)  on  the  escendin*? 
branch  of  the  original  trajectory  is  equal  to  y(t0)  -  K  ; 
this  tine  we  denote  by  sa.   There  is  a  time  at  which  y{t) 
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on  the  descending  branch  of  the  original  trajectory  ij 
equal  to  y(t0)  -  K   ;  this  time  we  denote  by  s<i. 

"Now  let  T  be  another  trajectory  (defined  by 
functions  x(s),  y(s),  etc.)  having  a  single  point  in  common 
with  the  yz-plane,  at  which  point  vx  >  0  and  y  >•  y(to)-  ffc'« 
On  the  ascending  branch  up  to  height  y(t0)  -  Vc  t    a   is 
uniquely  determined  by  y  and  conversely,  and  x  ■  -p.   Let 
t"(s)  be  the  time  at  which  y  »  y(s).   Then 

dy(T(s))/ds  »  dy(s)/ds, 


or 


Hence 


vv  dT/ds  «  vv(s ) . 


(29.4) 


dT/ds  =  vv(s 


V. 


(29.5) 


The  equations  of  motion  become 

dVy/ds  =  -  EUyVy/Vy  -   L  9y)Vy/Vyp 

dvx/ds  =  (-  Eux  +  ax)vy(  >Ay, 
dp/ds  =  [  vxvy  (s)/Vy]sgn^s~t;c 
dvz/ds  =  (-  Eu2  +  az)vy(s)/vy, 
dz/ds  =  vzvy(s)/vy, 

V   dt/dS  =  Vy(s)/Vy. 

In  the  r  ight  members  the  functions  E  are  to  be  computed  as 
functions  of  x,  etc.  along  f,  with  5c,  y  replaced  by  -p  and 
y(a)    respectively.        A  similar    argument  applies  to  the   part 


136 


of  the   descending  branch  wi th  y  *  y(t0)    -   vc  ;   the   only  dif- 
ference  is   that  x  is   to  be  replaced  by  +p  instead  of  -p, 
Which  also  accounts  for  the  factor  sgn    (s-t0)    in  the  third 
of   the    equations    (29.5). 

"There  remains  the  portion  of  t  with  y  >  y(t0)-*C, 
As  an  auxiliary,   we  use   polar  coordinates    (p,  "£)   to  replace 
(x,y),       The  angle  0  is  uniquely  determined  by  s,    and 

(29.6)  cos  0(e)   -  x(s)/     [x(s)2  ♦  {y(s)-y(t0)   *  l^p]  * 
sin  0(8)    -  y(s)/    [x(s)2  ♦  {y(s)-y(t0)   *te}2]  i/% 

At   tine   t  at  which  y   >y(t0)    -  Ire  $   the  rectangular  co- 
ordinates   (x,y)    transform  into    (p,0~),   where 

(29.7)  ©~{t)   -  0"(s(t)). 

By  elementary  computations  we  find   that 

dp/dt  ■  vx  cos  0  +   Vy  sin  0", 

(29.3)  ,  •'  > 

dO/d^  =  -    (vx/p)    sin  5  +    (vv/p)    cos  0\ 

By  applying  the  seme  argument   to  T,   we  find   that  the   second 
of  equations    (29.3)    is   valid  if    the  dashes  are   omitted  and 
T  replaced  by  s.       '.Ye  now   add  the  assumption    (retained 
throughout  the   rest  of  the  proof)    that   rc    is   small  enough 
so  that  Vx  remains  positive  between  Sq  and   s<3 ;   leter  we 
shall   reduce    *C    still  further.        Then    (vx/p)    sin  0   is  posi- 
tive  for  sa  <  s  <  sa.        For  s  <  t0  we  know   that  Vy  is  nega- 
tive,   while  for   s  >  to    it   is  positive.        At   s   =  t0, 
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cos  Q    (s)    conges  from  negative   to  positivo.        30 
(Vy/p)    cos   9   is  negative  except  at   to.        Thus  dQ/ds  is 
negative  for  Se  «  s   »  sdi    and   so   s  is  a  single-valued 
function  of  9  for  -  *T  »  9  »  0. 

"Furthermore,   by  making  T'  sufficiently  close   to  T 
in   the  usual  norm    (that   is,    by  making  the  norm  No  of  the 
disturbances  small),   we   can  make  the  right  members   of    (29,8) 
differ  from  the   corresponding  expressions  for  T  by  an 
amount   less  than   half  the  minimum  of    (vx/l. >}  sin9-  (vy/p)  cos9, 
so  that  the   right  member  of   the   second  equation  is  surely 
negative.       Then  9    (t)   has  a   single-valued  differential, 
inverse  £(?).        By   (29.7), 

(29,9)      dt/ds  -    Lu  9(b)    cl-l/[d  ?(T)/dt] 

-vx(s)p   sin   i){s)   +   vy(s)p"  cos  9(s) 
-vxp(s)    sin  ?♦  Vyp(s)    cos  © 

The   equations   of  motion  for     sa  *  s     '"  sa  now   become 
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(29.10 


dvy/ds  *    (-Euy  -   g  +    ay)dt/ 
dvx/ds   =    (-Eux  +   ax) dt/ds, 
dp/ds   »    [  vx  cos   9(s)    +  ?y  sin  9(s)   ]   dt/ds, 
dv2/ds  -    (-Eu2  +   a2)   dlT/ds, 
dz'ds   =  vz  d^/ds, 
dt/ds   =  dt/ds, 
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where  in  the  right;  members  dtVds  is  to  he  replaced  by 
its  value  from  (29.9),  y  by  p  sin  ©(s)  ♦  y(t0)  -  K   , 
£  by  p  qos  ©(s),  and  the  sine  and  cosine  of  Q(s)  are  in 
turn  to  be  replaced  by  their  values  from  (29. 6).* 

30,  Another  New  Nora \> 

Theorem  (89.2)  leads  us  to  conjecture  that  the 
existence  of  the  second  differential  of  the  mapping  (6.1) 
could  be  established  with  the  use  of  the  norm  N2y.   put 
an  example  which  will  be  given  in  $  32  shows  that  this  is 
not  possible. 

We  shall,  however,  demonstrate  the  existence  of 
the  second  differential  relative  to  a  new  norm— presently 
to  be  defined  and  to  be  designated  as  Ngy.   This  new  norm 
is  not  as  objectionable  as  the  norm  N2tJ  f°r  }t  does  notj, 
as  does  N2t»  include,  the  absolute  values  of  the  first  and 
second  partial  derivatives  of  the  disturbance  functions 
with  respect  to  altitude  y  at  every  point  of  R0.   On  the 
other  hand,  it  is  not  as  desirable  a  norm  as  A/gyf  because 
it  does  include  (whereas  Ngy  does  not  include)  the  absolute 
values  of  the  first — but  not  the  second— derivatives  Of 
the  disturbance  functions  with  respect  to  altitude  at 
points  of  Rp  which  lie  in  an  arbitrarily  small  neighborhood 
of  the  summit. 
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The  npw  norm  we  define  as  follow*: 

(90.1)  Definition:     Let  k  be  any  positive  number,   end  let 
Uk  denote  the  k-neighborhood    [in  Euclidean  7-epace] 
of  the   summit    (x©,    y0,    z0,   t0,    y^,   vy0,    yZq).       J"0** 
<fc|ty  disturbance 
<j\  s   (Ax0»    A  y0»    A  z0t    Atp,    AvXQ,   Avyo,-Avzo; 

4p  ,   £  #  ,   A  wx,    4  wz) 

k 
let  NgY    (A  )    be   the   larger  of  the   two  numbers; 

(i)  A/ay  (£); 

(ii)    sup    [(Afi)2  ♦    (A  #   .)2 

for  all  o^i  of  R0  in  Ufc,  in  which  qi  is 
replaced  by  x,  y,  z,  t,  vx,  vy,  v2  in  turn. 


31.   Existence  of  the  Second  Differential  When  the  Norm 

7 — * ™~~ ■ 

Is   N3y. 


We  now  state   our  theorem    m  the   second  differen- 
tial   comparable   to   that  of  fccShane         "he  first  differen- 
tial,   (29.2). 

(31.1)   yheorem;     Let  T  be  an  original  trajectory  on  which 
(£9,1)    is   satisfied.       Let   the  summit     f  T  corres- 
pond to  t  =   t0,    and   let  t3  be   any  number  such  that 
to  ■*  t3  <  tg.       Then  to  each  disturbance 
A  -(Axo^yoi^Zo^to^v^o^Vyo^Vao  ;4f,A9,4wx,Awz) 
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"there   correspond  two   sets   of  six  functions, 

1?  My)   -  v?i(y  iZ) 

and 

J  My)  -J  My;*,*  >, 

(i  -  1,    .   .     .    ,    6),    y(t2)   ^  y  ^  y(t3),   with  the   fol- 
lowing properties.        The    functions    f?    (y)   are  the 
same   as    those   so  designated   in    (29.2)    and  hence  have 
the   properties   stated  in   the    conclusion  of  that 
theorem.        With  the    usual  norm    (7.6),   the    y     (y) 
are   quadratic  functions  of   the  disturbance   2  and   are 
continuous   on   the    interval    fyCtg),    y(tg)J      when   A    is 
measured  by  the  norm  Niy(A).        If  we   assign  to  vy> 
vx»    vz»    *»    z,    t  the    alternative   symbols    q   ,    ...    ,    q6, 
respectively,    then    "J"^)    *s  the   second  differential 
of      q     in  the   sense   that   to  each  positive   €  there 
corresponds   a   positive     <?(  €  )    such  that,  whenever 
A    belongs   to   (j>yy  and  has   norm  N2y     •<£(€  ),    the   in- 
equality 

(31.2)       I  Aq(y)    -    I?  (y)    -     (1/2)   ^  (y)  |    -  <£  Nly  N^yf 

[y(t2)  -  y  -  y(t3)], 

is  satisfied. 

The  quadratic  character  and  continuity  property 
of  the  functions  3- My)  are  consequences  of  (16.1). 
In  order  to  establish  (31.2),  we  require  first  a  somewhat 
more  precise  estimate  for  the  constant  Kg  in  (11.31) • 
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In  his  proof  of  (11.3)  Mcshane  defines 
e^s)  ^Aqi(s)  -  ^(s)  [  not  the  same  functions  «3 
those  in  (20.1)  ]  and  then  deduces  the  inequality 

(31.3)  d|e(s)l/ds  -  A(s)  |e(s)  |+K0N0U  q0t  A  <f)Nj.B(4  qp,A<T)» 

valid  whenever  |e(s)|  >  0,  in  which  K  is  constant  and  A(s) 
is  the  least  upper  bound  of  Fypc(£f(s,  q(s  ) )  )  t-f*    j  (s,  q(s  ))u^u«J 

for  unit  vectors  u.   He  then  states  that  lemma  (9.3), 
applied  to  (31.3)  shows  the  existence  of  a  constant  Kg  such 
that   |e(s)|  -  KoN0  (  A   q0,  &*P)   Nis(  A  q0,-  H°) ,  which  is 
equivalent  to  (11.31 ).   The  details  of  this  application 
of  (9.3)  are  omitted  in  [ifcs]  .   Upon  supplying  them.,  v, 
find  that 

/  «*•  sx 

(31.4)  |e(s)  |*  exp  j   A(t)dt ..  <    le  (si)  |*  [kNoNis  [expJ*A(*)drJdt}f 

in  which  we  have  written  N0      .  c  for  the  first  and  sec- 
ond norms,  respectively,  occurs      the  right  member  of 
(31.3).   But  ei(si)  -  0,  so  (31.4)  reduces  to 

s 

(31.5)  Uq(s)   -  i?(s)|    -  K[  J{exp  J"  A(fK      }   dt]   N0Nls. 


On  the  original  trajectory,    s  and  t  ar«  identical. 
So  the  equations    (29.10)   hold  when  the   dashes   are  omitted 
and  dt/ds  is   replaced  by  unity.       Let  k  be   a  positive  number, 
and  let  Uk  be  the  k-neighborhood  of  the  summit.       By  the 
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hypothesis  (29.1),  returning  to  the  F ( <f  (s»q)  notation  in 
(29.10),  for  all  q  in  a  neighborhood  U  of  the  summital 
values  q(t0)  the  functions  <f  qi  and  Frof(  <f  (a,q) )  remain 
bounded.   Clearly  we  lose  no  generality  by  assuming  UtfUfc. 
Then  there  exists  a  positive  number  K0   so  small  that  on 
the  corresponding  interval  sa0  *  s  =  s^o   "the  quantities 
q(s)  on  the  original  trajectory  remain  in  U,   Furthermore, 
we  suppose  that  the  norm  N0  of  the  disturbance  is  small 
enough  so  that  the  quantities  q(s)  on  the  disturbed  tra- 
jectory also  remain  in  U.   Let  A  be  an  upper  bound  for 
A(s)  on  U.   Then  for  each  positive  tt  «  KT0  and  the  cor- 
responding s0  and  sd  determined  by  MT  ,  it  follows  from 
(31  o 5)  that  there  exists  a  constant  K<±  determined  by  A, 
K,  and  fC0  such  that 

(31.6)  |Aq(sd)  -  f?(sd)  [  -  K4  Ufl  -  8a)  N0  Nis. 

How  let  €   be   any  positive  number.       Let  KQ  be 
the   positive  number  whose  existence  is    assured    by    (£8.3). 
We  now   choose   a  positive  number    ta    an  all   enough  to  satisfy 
If    =    lt0  and   K  =  y(t0)    -  y(tg),   as  well  as  the  additional 
requirement 

(31.7)  sd   -    se  "£    €/4  KoK4. 

(This  is  possible,  since  the  left  member  of  (31.7)  tends 
to  zero  with  K   . )   If  (4qo»A(P  )  is  a  disturbance,  we  de- 
fine ^£f>(s>0j  to  be  4A?(s,q(s))  for   sa  ^  s  ^  sd  and 
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to  bo  ^(Pdjq)  otherwise.   £et  A  z  f  *  Acf-   A  x<f   ,  so 
th$t 

(31.8)  (4q0>*f  >  "  (A<*o»  A  I?)    ♦  (  °»  *  a*5 '  )• 

On  the  original  trajectory  T  the  disturbances  A<f  and 
Aitf   are  identical.   Hence   [Cf.  (8.8)] 

(31.9)  r)(s;Qo,(f';Aq0,  *  if  )  ,  ^(s;q0,9  J  Aq0,*^  ) 
and 

(31.10)  J  (s;  Ai<f  ,  AX(^  )  ,  3-(s;  ^,  Acp  ) 

for  t3  £~  s  -  t2.   From  the  definition  of  A  ±(f  f   We 
obviously  have 

(31.11)  N0(.  "o,  *lf  )  -  N0(  Mo»  ^f)  . 

From  (31.8) ?  (31.11),  .    he  triangle  inequality  it  fol- 
lows that 

(31.18)    N0(  0,  A2tf)  ~  £  L     :  )fAf  )  . 

On  the  interval  sa  ■  s  «~  *£  the  functions  A  ^  (£ 
are  dependent  on  s  alone,  so  tney  c     bute  nothing  in 
(ll.Bii)  or  in  (12.7).   Thus 

(31.13)    sup  1^1^^(5,0)1    (ti^s  =  t2) 

-   sup    l-A^f  qj    (s,q)|       (ti^s   ^"sa    sad  Sfl^s  ^~  t%) , 
and 
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(3i.i4)    eup  I  A  x  q  qjqk(s,q) I      (ti  *■  a  -  t2) 

-  sup  I  /if  qjqk  (s,q)   I 

(tx  "  s   *   sa,   and     Bfl  -   s  ■  tg) . 

But  on  the  intervals  last  indicated  the  variables  q«J  are  all 
the  variables  Vy,  .  .  .  ,  t  except  y,   Hence  the  quantities 
in  (31.13)  and  (31.14)  cannot  exceed  the  supremuni  of  the 
partials  of  the  A  (f   with  respect  to  all  the  variables  ex- 
cept y.   From  this  and  (31.11)  we  see  that 

(31.15)    Nls(Aq0,  Arf)   t  Nly(  Aqot  Af  ) 

and 

(31.15)    N2S(  A<l0,*lrf  )  ^N2y(*qo»  *?) 

Besides   the   original  trajectory  T  and  the  dis- 
turbed trajectory  T  we   shall  also  consider  a   trajectory 
Tis    Q.  «  qi(s)    computed  with  the  same   initial  values 
q0  ♦   A  qo  as  those   of  f  but  with  functions    (f  ♦    A  l  ^     in- 
stead  of    lf+    A  «p  .        Thus  we  may   consider  that  T  is  ob- 
tained from  Ti  by  applying  the   disturbance    (   0   ,    A  2  *P  ) • 

On  the  original   trajectory     T     we  have 
ql    (s)   ■  Vy(s)   <  0,    (sa  ■  s  -  t2)»       Hence  on  this   interval 
-ql  has  a  positive   lower  bound,  which  we  denote  by     inf    (-q^)» 
By  applying  the   continuity  theorem   (28.6)    to  the  equations 
(29.5),    (29.10),   we   see   that  there  exists   a  positive     &i 
such  that  if 
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<-  *c 


(51.17)  N0(Aq0,  Af)  <   £lB 
then 

(31.18)  qi(s)   arxi    q(s)   are  in  U  for      aa  ■"  s  ■  sd, 
and 

(31.19)  lq1(s)-q1(s)  |    and    | "q1  ( s ) - <J< s )  |      are   less  than 

inff-q1)/^     for   sd  -   s  *  t2. 

Henceforth  we  assume   that    (31.17)    it    satisfied. 

The   trajectories  Ti  and  T  are   identical  for 
ti  *   s   »   s6      since  T  is  obtained   from  Ti  by  disturban. 
(°»    A  Z^   )   v      cl    ^re   zero   for     s<sa.        If   (31.17)    holds, 
so  does    (31.18),    and  equently,    by    (31.6),   with   q^ts) 

replacing  q(s),   we  have 

(31.20)  ki(sd)    -   q(sd)    -     y?2    (sd)| 

-  K^4    (sd-sa)   No(0,    A  2<        Nis(0,    A8cp)t 

in  which      f?8(s)    5  £  (s;    q0,  f  ;     C,    A  2  <p  ) . 

Moreover   in  McShane's  deduction  csl ,   #  6)    of 

(31.3),   when    AF  =  0  the  factor  Nis    (  Aq0,A^)    in  the 
right  member  of  that   inequality  is  necessary  only  because 
of  the   presence   of  the  factor 

(31.21)  J        &fQy,    (s,    q  *     t   A^)    dt 
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in  the   right  member  of  his   equation   for  <!(£*  (s)  )/d«:.        In 
our   application  of   (31.6)   above      qi(s)    replaced  q(s);    so 
in    (31.21)    the   integrand  is    A2tfqj(s»    Ql  ♦    (<1  -   41   )    )> 
which  for     0  ■    ~  -  1   and  sa  ■  s   ■   sd     does  not   exceed  a 
constant  multiple   of     N2y    (Aqo,^^  )>    and  is  identically- 
zero  both  when   si   ■   s  <  sa   and  when     sa  <  s  ■   82»        Conse- 
quently  we  may  replace     N]_g(0,  ^2^)    bv  N2y    ^A<lo»AcC    )    in 
(31.20).        But  then,    from    (31.12)    and   the    fact  that 

I?  2^s^    *  °   sillce      ^i^f   =   ^  ^P   on  the   original  trajectory, 
we  see    that    (31.20)    can  be  replaced  by 

(31.22)  |q1(sd)-q(sd)|^  2K4(sd-sa)    ^fAq^  Af  )N2y(  Aq0,^)  , 

using   (31.11)   and    (31.13).        It  then  follows   from   (31.7) 
and    (31.22)    that 

(31.23)  ki(sd)-q(sd)|    -    (ff/2K    )   N0(Aq0,A<p)    N2y  (  A  q0,  A<p)  , 

If  (31.17)  holds,  then  so  does  (31.19),  and  we  can 
apply  lemma  (28.3)  with  q  in  place  of  q2  and  sd  in  place 
of  S3.   From  this  and  (31.23)  we  find  that  the  inequality 

(31.24)  |q1(s)-q(s)i  =  (  G/Z)    N0(Aq0,4<f)  N2y(Aq0,^) 

is  valid  for  s<i  *  s  ■  t2,  end,  in  particular,  for 

1 3   *  s  ■  t2* 

By    (15.10)    there   exists   a  number  K3  such  that 
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(31.25)   |qi(s)-q(s)-  |?(s)-(l/2)  J  (s  )  | 

■K3Nis(M0»4lf)  [N28^<lo.Al^)  Y> 


t3  »  s  »  t2,  in  which  •?  (s )  is  either  member  of  (31.9)  and 
TJ"  (s)  is  either  member  of  (31,10) • 

We  now  define 

(31.26)  &{e)    =  min    [  S  lf     €/2K3J. 
Then,    if 

(31.27)  N2|(  ^  q0,  A<p  )   <    &  (€)    , 

it  follows,  by  (31.16),  (31.26),  (31.25),  and  (31.15), 
that 

(31.28)  kiU)-q(s)  -  f?(s)-(l/2)  y  (s)  | 

-  (<?/2)  Nxy(^q0,Af)  N2y(Aq0,^), 

t3  ^~s  ^t2.   But  if  (31.27)  is  valid,  so  is  (31.17),  by 
virtue  of  (31.26)  and  the  inequality  No  =  Niy  ■  N2y« 
Hence  (31.24)  also  holds,  and  from  (31.28)  and  (31.24)  it 
follows  that 

(31.29)  |q(s}-q(s)-  1J(s)-(l/2  J(s)  |-  €  Niy(Aq0,^f)N2y(Aq0,4<p) 

whenever  (31.27)  is  satisfied.  If  we  now  replace  s  by 
«(y)>  y(t2)  *  y  »  y(t3),  we  have  (31.2).  This  completes 
the  proof  of  our  theorem. 
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32.  A  Disturbance  Which  Has  Mo  Second  Differential 
Effect  Relative  to  the  Norm  Ngy. 

In  this  section  we  cite  an  example  to  show  the 
impossibility  of  proving  the  existence  of  a  second  differ- 
ential of  the  mapping  (6.3)  relative  to  the  norm  N^y. 

We  shall  assume  an  original  trajectory  such  that 
air-density  (»  and  air-temperature  0  are  constant,  2  is 
positive  except  when  air-speed  u  is  zero,  wx  *  wz  »  0, 
and  with  initial  values  x<>"  yo  ■  Zo  ■  vXo  ■  vZo  ■  0, 
v  0  >  0.   Let  Y  be  the  summital  height  on  this  original 
trajectory.   Let  £   be  a  disturbance  consisting  only  of 
the  following  abnormal  conditions: 

AYo  a  r,   for  each  real  r; 
o  ,    y  -  Y 


Awx  -  \      y  -  Y,   Y  «  y  -  Y  «■  2r 
2r,     Y  +  2r  *  y 

For  4  y0  >  o,  if  sa  and  sd  are  the  times  of 
passing  altitude  Y  on  the  ascending  ord  descending  branches 
of  the  disturbed  trajectory,  we  have 

(32.1)   sd  -  sa  -  (8  *yJs)1/Z 

+  (infinitesimal  of  higher  order)  , 

Near  the  summit  on  the  disturbed  trajectory  vx  is  ne«r0, 
ux  is  nearly  -r,   E  is  approximately  equal  to  its  value 
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]$o  at  the  summit,  and  Vx  is  therefore  nearly  E0r.   Hence 
at  3d  we  have 

(32.2)  vxlsd)  *  EQ  r  (sd-sa) 

+    (infinitesimal  of   higher  order). 

It  follows  from  theorem    (29.2)   that,    at  y  -  0  on  the  de- 
scending branch   of  the   disturbed   trajectory,    x  is  nearly 
proportional  to     vx{sa).        So  for    4  y0  >  0     we   hove,    com- 
bining   (32.1)    and    (32.2) 

(32.3)  AX  -  (constant).  r3/2, 

the   constant  factor  being  non-zero.   But  since 
N2y  (A)  -  |2r|,  (32.3)  is  equivalent  to 

(32.4)  AX  -      bant)  [N2y  (a)  ] 

the   new  constant  factor      being  non-zero. 

In  case  A  yp  <  0,  however,  the  summital  height 
on  the  disturbed  trajectory  is  Y   AyOJ  because  A  wx  =  0 
for  y  ■  Y.   It  follows  that  in  this,  case  ax  *  0. 

Hence  dX  =  dx  (  K |  y  -  0;  s  I  and 
d2X  s  d2x  (  £ |  y  =  0}  =0.   But  then  we  s2  -  from  (32,4) 
that 

(32.5)  |  A  X  -  dX  -  |  d2x|  =  K  CNgy(X)  33/£ 

in  which  K  is  a  non-zero  constant.        The  ratio  of  the  left 
member  of    (32.5)    to      CN2y   (  &  )12     therefore  does  not  tend 
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to   zero  with  Ngy(A  ),    so  the  second  differential  relative 
to  this  norm  does  not  exist# 
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